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Summary 


A  novel  actuator  for  spacecraft  attitude  control  with  liquid  flywheel  is 
presented.  The  main  characteristic  of  this  new  concept  of  reaction  wheel  is 
that  a  conductive  liquid  rather  than  a  solid  mass  is  accelerated  to  change 
the  angular  momentum  of  the  equipment  and,  as  a  consequence,  to  provide 
torque  to  the  spacecraft.  The  conductive  liquid  is  accelerated  using  a 
conductive  Magneto- Hydro- Dynamic  (MHD)  pump.  Two  different  configu¬ 
rations  of  the  device  have  been  studied  on  the  basis  of  the  optimization  of 
the  dimensionless  moment  of  inertia  and  the  minimization  of  the  viscous 
shear. 

A  2-dimensional  Finite  Difference  Hybrid  Model  (FDHM)  has  been 
developed  on  the  basis  of  the  MHD  set  of  equations  under  the  hypothesis  of 
low  Magnetic  Reynolds.  The  model  solves  numerically  the  time  dependent 
axially  symmetric  problem  of  an  electrically  conductive  liquid  rotating  in 
a  torus  with  rectangular  cross  section  due  to  the  interaction  of  a  radial 
magnetic  field  and  an  axial  electric  field.  The  electric  side  of  the  problem 
has  been  solved  by  means  of  the  node  method  applied  to  a  network  of  elec¬ 
tric  resistances  and  voltage  generators  representing  the  back  electromotive 
voltage  induced  by  the  spinning  liquid  through  the  magnetic  field.  The 
fluid- dynamics  side  of  the  problem  has  been  solved  using  a  Crank- Nicols on 
method  over  a  non  uniform  and  collocated  grid.  The  grid  generator  has 
been  written  to  be  sensitive  to  the  Hartmann  number  of  the  problem. 
Several  simulations  has  been  made,  over  600,  in  order  to  test  the  FDHM 
and  in  order  to  outline  the  performances  of  the  device.  The  model  has 
been  tested  numerically  checking  the  accuracy  of  the  computed  dynamic 
and  electric  quantities. 

A  Lumped  Parameter  Model  (LPM)  has  been  developed  using  two  sta¬ 
tionary  analytical  solutions  of  the  MHD  set  of  equations.  The  LPM  has 
been  compared  with  the  results  obtained  with  the  FDHM,  showing  good 
agreement  in  the  estimation  of  the  coefficients  of  the  model  in  both  the 
cases  of  high  and  low  electrically  conductive  liquid.  Finally  a  parametric 
performance  analysis  has  been  conducted  by  means  of  the  best  set  of  formu¬ 
las  to  compute  the  LPM  coefficients  in  order  to  outline  the  performances 
of  the  device  for  several  geometric  dimension  of  the  torus  containing  the 
conductive  liquid  and  for  several  values  of  electric  conductivity  of  the  liquid. 
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Chapter 

Introduction 


The  ability  to  control  spacecraft  attitude  and  the  reliability  of  the  attitude 
control  system  are  key  components  of  space  missions.  Also,  future  Earth  obser¬ 
vation  missions,  as  well  as  astronomical  observation  missions,  require  satellites 
which  are  more  agile  and  with  high  pointing  stability  over  the  higher  lifetime 
possible. 

A  number  of  techniques  have  been  developed  and  rehned  to  underpin  attitude 
control  and  re-pointing,  but  basically  since  the  beginning  of  the  space  explo¬ 
ration  era  all  techniques  used  for  high  performance  attitude  control  exploit  the 
use  of  reaction  and,  in  some  cases,  momentum  wheels. 

The  technology  improvement  during  this  time  allowed  reaction  wheels  to  im¬ 
prove  their  performance  in  terms  of  control  resolution  and  reliability  and  effi¬ 
ciency  of  the  electronic  part  but  there  are  some  drawbacks  related  to  standard 
reaction  wheel  design  that  were  not  improved  over  time.  In  particular  we  can 
mention  some  of  the  more  important  drawbacks  like: 

•  Standard  reaction  wheels  use  moving  mechanical  parts  which  reduce  the 
overall  reliability  of  the  equipment  and  increase  its  overall  complexity. 

•  Jitter  and  other  unbalances  associated  with  moving  masses  lead  to  un¬ 
wanted  high  frequency  disturbances  on  the  satellite  itself  which  could  be 
critical  for  demanding  high  resolution  observation  missions. 

•  The  presence  of  electric  motors  in  standard  reaction  wheels  induces  a 
ripple  disturbance  on  the  torque  provided. 

•  The  torque  provided  by  standard  reaction  wheels  is  typically  non  linear 
in  the  neighbourhood  of  the  zero  angular  velocity. 

In  the  last  years  some  articles  have  been  published  regarding  the  possibilities 
to  use  active  fluid  loops  as  actuator  for  spacecraft  attitude  control,  [1]  [2] .  This 
kind  of  actuators  uses  hydraulic  pumps  to  accelerate  a  low  viscosity  liquid  in¬ 
side  a  torus  with  circular  cross  section.  These  devices  have  the  advantage  to 
reduce  the  mechanical  unbalances  typical  of  classical  reaction  wheels  with  solid 
flywheel  but  they  do  not  increase  the  reliability  of  the  actuator  because  of  the 
presence  of  a  spinning  electric  engine  in  the  hydraulic  pump. 

A  novel  actuator  for  spacecraft  attitude  control  is  proposed.  The  main  char¬ 
acteristic  of  this  new  concept  of  a  reaction  wheel  is  that  a  conductive  liquid 
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Cap.  1  Introduction 


rather  than  a  solid  mass  is  accelerated  to  change  the  angular  momentum  of 
the  equipment  and,  as  a  consequence,  to  provide  the  torque  to  the  spacecraft. 
The  conductive  liquid  inside  an  annulus  with  rectangular  cross  section  is  ac¬ 
celerated  using  the  Lorentz  force  originated  from  the  interaction  between  an 
electric  current  flowing  through  the  fluid  and  a  magnetic  held.  In  this  way 
the  electric  engine  can  be  removed  with  the  advantage  to  remove  any  ripple 
disturbance  and  to  improve  the  overall  reliability  of  the  equipment. 

In  the  past  and  in  particular  in  the  last  years  the  interest  in  this  kind  actuators 
increased,  [3]  [4]  [5]  [6]  [7]  [8]  [9]. 

The  lack  of  accurate  mathematical  models  describing  the  behaviour  of  the  de¬ 
vice  brought  this  study  to  start  from  the  basis  of  the  magneto-hydro-dynamics 
problem  to  obtain  a  numeric  model  able  to  model  the  real  physics  of  the  prob¬ 
lem.  The  model  is  then  used  to  better  characterize  the  actuator  and  to  have 
information  about  its  performances  and  limitations. 
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2 


Configurations:  Distributed  Magnetic 
Field 


In  this  chapter,  a  preliminary  analysis  of  the  optimal  configuration  for  a  re¬ 
action  wheel  with  distributed  conduction  Magneto- Hydro- Dynamics  pumps 
is  presented. 


2. 1  Working  principle .  5 

2.2  Configurations  of  the  device .  6 

2.3  Dimensionless  Moment  of  Inertia .  7 

2.4  Comparison  between  configurations .  8 
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Cap.  2  Configurations:  Distributed  Magnetic  Field 


2.1  Working  principle 

The  main  feature  of  this  new  concept  of  reaction  wheel  is  that  an  electrically 
conductive  fluid  rather  than  a  solid  mass  is  accelerated  to  change  the  angular 
momentum  and  to  provide  a  controlling  torque  to  the  spacecraft.  The  Lorentz 
force  is  used  to  accelerate  the  conductive  liquid  [10].  This  driving  force  orig¬ 
inates  from  the  interaction  between  a  magnetic  held  and  an  electric  current 
howing  through  the  huid,  Figure  2.1. 

The  externally  applied  magnetic  held  is  provided  by  permanent  magnets  be¬ 
cause  the  use  of  magnetic  coils  would  increase  the  total  power  consumption  of 
the  actuator  and  this  is  not  convenient  for  a  device  whose  use  is  addressed  to 
space  application,  where  power  resources  are  usually  limited.  The  magnetic 
held  in  this  section  is  considered  distributed  to  all  the  volume  where  the  liquid 
is  placed. 


Figure  2.1:  Section  of  an  annulus  with  applied  electric  field  and  radial  mag¬ 
netic  field  B  (r). 

The  spin  of  the  liquid  can  be  controlled  in  intensity  and  direction  by  means  of 
the  electric  current  sent  to  device.  The  inversion  of  the  direction  of  the  electric 
current  switches  the  direction  of  the  Lorentz  force  acting  as  body  force  on  the 
huid. 
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The  liquid  is  placed  in  an  annulus  with  a  rectangular  cross  section  and  it  ro¬ 
tates  along  the  azimuthal  direction.  The  Lorentz  force  acts  in  the  direction  in 
which  the  electrically  conductive  fluid  is  free  to  flow. 


2.2  Configurations  of  the  device 

The  maximum  Lorentz  force  obtainable  from  the  interaction  between  an  elec¬ 
tric  current  and  a  magnetic  field  appears  when  the  two  fields  are  mutually 
perpendicular  [11].  This,  together  with  the  advantage  to  have  the  conductive 
liquid  placed  in  an  axially-symmetric  container  leads  to  two  possible  arrange¬ 
ments  of  magnetic  field  and  electric  current  and  two  different  configurations 
for  the  MHD  actuator,  Figure  (2.2): 

•  Axial  magnetic  field  configuration-,  this  conhguration  takes  into  account 
an  electric  current  along  the  radial  direction  r  and  a  magnetic  field  along 
the  axial  direction  2;  [5]  [12]  [4]. 

•  Radial  magnetic  field  configuration:  this  conhguration  switches  the  hrst 
conhguration  considering  an  electric  current  along  the  axial  direction  z 
and  a  magnetic  held  along  the  radial  direction  r  [11]  [3]. 


Figure  2.2:  First  (top)  and  second  configuration  (bottom)  for  electric  field  E 
and  magnetic  field  B 
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Both  the  conhgurations  provide  a  force  in  the  requested  direction.  The  advan¬ 
tage  to  use  one  of  two  conhgurations  has  been  evaluated  on  the  basis  of  the 
maximization  of  the  dimensionless  moment  of  inertia. 


2.3  Dimensionless  Moment  of  Inertia 


The  dimensionless  moment  of  inertia  gives  a  numerical  indication  of  how 
efficiently  a  mass  is  used  to  generate  angular  momentum.  To  obtain  the  value 
of  Iz  for  a  hollow  cylinder  we  need  to  divide  the  dimensional  moment  of  inertia 
by  the  mass  and  the  squared  external  radius  of  the  cylinder.  We  obtain: 


h 


m  (fg  -1-  rf ) 


2 


(2.3.1) 


where  is  the  outer  radius  of  the  hollow  cylinder,  is  the  inner  radius  of  the 
hollow  cylinder,  is  the  ratio  Ti/rg  and  m  is  the  mass  of  the  liquid. 


We  can  see  that  Iz  is  independent  from  the  height  of  the  cylinder.  It  assumes 
its  maximum  value  when  the  inner  radius  approaches  the  outer  radius,  i.e.  for 
thin  cylinders. 
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Dimensionless  Moment  of  Inertia  -  Iz 


Figure  2.4:  Dimensionless  Moment  of  Inertia 


2.4  Comparison  between  configurations 

We  start  pointing  out  that  in  both  the  configurations  the  dimension  parallel  to 
the  magnetic  held  is  constrained  in  order  to  have  a  homogeneous  magnetic  held 
in  the  air  gap  of  the  magnetic  circuit.  It  means  that  the  height  h  of  the  device 
in  the  hrst  conhguration  and  the  radial  length  of  the  cross  section,  L  =  rg  —  rj, 
in  the  second  one  cannot  be  greater  than  a  certain  value. 

In  Figure  (2.5)  it  is  evident  that  where  the  dimensionless  moment  of  inertia 
has  its  maximum  we  have  the  minimum  of  the  dimensional  moment  of  inertia. 
For  the  axial  magnetic  held  conhguration,  if  the  device  has  an  almost  optimal 
dimensionless  moment  of  inertia  it  is  not  possible  to  have  a  great  dimensional 
moment  of  inertia  due  to  the  limitation  of  the  height  h  of  the  device. 
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X 10  ’  Moment  of  Inertia  -  F  —  [kgrrd] 


Dimensionless  Moment  of  Inertia  -  F 


Figure  2.5:  Dimensional  and  Dimensionless  Moment  of  Inertia  -  I'’*  Config¬ 
uration 


For  the  radial  magnetic  field  configuration  the  device  can  be  characterized  by 
an  almost  optimal  dimensionless  moment  of  inertia  and  it  it  is  possible  to 
obtain  the  desired  dimensional  moment  of  inertia  modifying  the  height  of  the 
device,  not  constrained  in  this  configuration. 


X 10  ’  Moment  of  Inertia  -  F  —  [kgm^]  Dimensionless  Moment  of  Inertia  -  F 


Figure  2.6:  Dimensional  and  Dimensionless  Moment  of  Inertia  -  2^*  Config¬ 
uration 


We  could  think  to  stack  several  annuli  according  to  the  axial  magnetic  field 
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configuration  in  order  to  obtain  the  desired  moment  of  inertia  while  optimiz¬ 
ing  the  dimensionless  moment  of  inertia,  Figure  (2.7),  as  already  studied  in 
[5].  This  approach,  anyway,  would  increase  the  contact  surface  between  the 
liquid  and  cavity,  and  then  the  viscous  shear,  more  than  how  it  would  increase 
choosing  the  second  conhguration. 


Figure  2.7:  MHD  Fluid  Loop  -  1®*  Configuration 


We  can  conclude  that  the  second  conhguration,  characterized  by  a  radial  mag¬ 
netic  held  and  an  axial  electric  held,  maximizes  the  dimensionless  moment  of 
inertia  while  optimizing  the  contact  surface  of  the  conductive  liquid.  It  rep¬ 
resents  a  more  easily  scalable  conhguration  in  which  the  required  moment  of 
inertia  of  the  liquid  hywheel  can  be  met  just  varying  the  height  of  device. 

The  device,  according  to  the  selected  conhguration,  has  electric  plates  on  the 
top  and  on  the  bottom  sides  and  two  concentric  cylindrical  magnets  generating 
a  radial  magnetic  held.  The  conductive  huid  is  placed  in  the  air  gap  of  the 
hollow  cylinder.  Figure  (2.8).  The  spin  direction  of  the  liquid  can  be  switched 
by  changing  the  direction  of  the  electric  held  between  the  two  plates. 


Figure  2.8:  MHD  Fluid  Loop  -  2”*^  Configuration 
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Chapter 


Assumptions:  Axially-symmetric  MHD 
problem  and  Low  Re^  approximation. 


In  this  chapter,  the  MHD  set  of  equations  describing  the  problem  are 
presented  in  details.  The  equations  are  set  to  solve  the  2-dimensional 
axially-symmetric  problem  of  an  electrically  conductive  liquid  rotating  in 
a  torus  with  rectangular  cross  section  under  the  interaction  of  a  radial 
magnetic  field  and  an  axial  electric  field. 


3.1  Complete  MHD  set  of  equations .  12 

3.1.1  The  Low  Magnetic  Reynolds  Approximation .  13 

3.2  Axially  symmetric  MHD  equations .  15 

3.2.1  Magnetic  Field .  16 

3.2.2  Navier- Stokes  equations  and  continuity  equation  for  the 

velocity .  17 

3.2.3  Electric  current  density  and  Back  electromotive  force  .  18 

3.2.4  Boundary  and  Initial  Conditions .  20 
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3.1  Complete  MHD  set  of  equations 

The  physics  of  the  Magneto-Hydro-Dynamics  (MHD)  phenomena  is  generally 
described  by  a  set  of  equations  taking  into  account  the  mutual  interactions 
among  the  velocity  field,  the  electric  held  and  the  magnetic  held,  [13]  [14], 
The  complete  MHD  set  of  equations  includes  Maxwell’s  equations,  (3.1.1) 
(3.1.2)  (3.1.3)  and  (3.1.4),  Ohm’s  law  applied  to  a  moving  conductive  medium 
(3.1.5),  together  with  the  huid-dynamics  equation  for  a  viscous  huid  on  which 
acts  an  electromagnetic  body  force  (3.1.6)  and  the  continuity  equation  (3.1.7). 


^  X.  dB 

(3.1.1) 

V  -B  =  0 

(3.1.2) 

V  X  B  =  47r/imJ 

(3.1.3) 

V- J  =  0 

(3.1.4) 

J  =  a(E  +  U  X  B) 

(3.1.5) 

(9U 

p—  +  p(U  ■  V)U  =  -VPn  +  z/pAU  +  J  X  B 

(3.1.6) 

<1 

o 

(3.1.7) 

where  E  is  the  electric  held,  B  is  the  magnetic  held,  J  is  the  electric  current 
density,  U  is  the  velocity  of  the  huid  and  /i,  p,  pm  and  a  are  respectively 
the  dynamic  viscosity,  density,  magnetic  permeability  and  electric  conductiv¬ 
ity  of  the  working  liquid.  In  order  to  simplify  the  set  three  assumptions  are 
considered: 

•  The  how  is  considered  laminar 

•  The  problem  is  considered  axially-symmetric  and  it  is  solved  over  the 
r  —  z  domain. 
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•  The  hypothesis  of  Low  Magnetic  Reynolds,  Low-Re^.,  is  applied. 

The  hypothesis  of  axially-symmetric  problem  is  a  direct  consequence  of  the 
distribution  of  magnetic  held,  considered  to  be  radial  all  over  the  annulus,  and 
consequence  of  the  presence  of  the  electric  armatures  extending  all  over  the 
upper  and  lower  surface  of  the  actuator.  The  hypothesis  of  laminar  how  is 
consistent  considering  a  limited  maximum  velocity  of  the  how  and  the  typical 
stabilizing  ehect  that  Lorentz  force  induces  in  MHD  hows.  The  hypothesis  of 
Low-Rcm  will  be  discussed  in  detail  in  the  next  section. 

3.1.1  The  Low  Magnetic  Reynolds  Approximation 

The  Magnetic  Reynolds  number,  Rcm,  gives  information  about  how  the  mag¬ 
netic  held  is  inhuenced  by  the  velocity  held  and  the  electric  current.  We  look  at 
the  ehect  of  a  prescribed  magnetic  held  on  the  how.  To  ensure  that  B  remains 
unahected  by  U  we  must  restrict  ourselves  to  low  magnetic  Reynolds  numbers: 

Rm  =  =  fim<j\Vo\L  <  1  (3.1.8) 

where  Uo  is  the  mean  velocity  of  the  huid.  However,  this  is  not  overly  restric¬ 
tive,  at  least  not  in  the  case  of  liquid-metal  MHD.  For  example  for  case  we  are 
going  to  analize  u  ~  \rv? j s  ,  L  ~  0.005  and  internal  friction  keeps  |Uo|  to  a 
level  of  around  O.Olm/s  ~  Im/s.  This  gives  Rm  ~  10“^  -i- 10“^. 

The  essence  of  the  low  Magnetic  Reynolds  Rm  approximation  is  that  the  mag¬ 
netic  held  associated  with  induced  currents,  |J|  ~  a  |Uo|  X  |B  I,  is  negligible  by 
comparison  with  the  imposed  magnetic  held. 

There  are  three  distinct  cases  which  commonly  arise. 

1  The  imposed  magnetic  held  is  static,  the  how  is  induced  by  some  external 
agent,  and  friction  keeps  u  to  a  modest  level  in  the  sense  that  |Uo|  rj/ 1 

2  The  imposed  magnetic  held  travels  or  rotates  uniformly.  This  induces 
a  how  which,  due  to  friction  in  the  huid,  is  somewhat  slower  than  the 
speed  of  the  held. 

3  The  imposed  magnetic  held  oscillates  extremely  rapidly,  in  the  sense 
that  the  skin-depth  5  =  {2/ is  much  less  than  the  characteristic 


13 


DISTRIBUTION  A.  Approved  for  public  release:  distribution  unlimited. 


Cap.  3  Assumptions:  Axially-symmetric  MHD  problem  and  Low  Rem  approximation. 


length  of  the  problem,  /  being  the  held  frequency.  The  magnetic  held  is 
then  excluded  from  the  interior  of  the  conductor  and  inertia  or  friction 
in  the  huid  ensures  that  |Uo  I  «/^ 


Categories  (1)  —  (3)  cover  the  majority  of  hows  in  engineering  applications. 
Typical  examples  are  the  magnetic  damping  of  jets,  vortices  or  turbulence  (1), 
magnetic  stirring  using  a  rotating  magnetic  held  (2)  and  magnetic  levitation 
(3).  In  the  cases  we  are  going  to  study  the  applied  magnetic  held  will  be  steady 
and  as  said,  the  Magnetic  Reynolds  is  Rm  1.  We  can  then  assert  that  all 
the  case  we  are  going  to  simulate  belong  to  the  hrst  case. 

We  now  discuss  the  simplihcations  which  result  in  the  governing  equations 
when  Rm  is  low  and  the  imposed  magnetic  held  is  steady.  Let  Eq,  Jo  and  Bq 
represent  the  helds  which  would  exist  in  a  given  situation  if  U  =  0,  and  let  e,  j 
and  b  be  the  inhnitesimal  perturbations  in  E,  J  and  B  which  occur  due  to  the 
presence  of  a  vanishingly  small  velocity  held.  These  quantities  are  governed 
by: 


V  X  Eo  =  0  (3.1.9) 

Jo  =  aEo  (3.1.10) 

~  dh  ,  , 

Vxe  =  -—  (3.1.11) 

j  =  a  (e  +  U  X  Bo)  (3.1.12) 

where  we  have  neglected  the  second-order  term  U  x  b  in  (3.1.12).  Now  Fara¬ 
day’s  equation  gives  e  ~  Ub  and  so  the  perturbation  in  the  electric  held  may 
also  be  neglected  in  (3.1.11).  Ohm’s  law  now  becomes: 

J  =  (J  (Eo  -|-  U  X  Bq)  (3.1.13) 

However,  Eo  is  irrotational  and  so  may  be  written  as  — V0  ,  where  0  is  an 
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electrostatic  potential.  Our  final  version  of  Ohm’s  law  is  therefore: 

J  =  a(-V0  +  U  X  Bo)  (3.1.14) 

while  the  leading-order  term  in  the  Lorentz  force  per  unit  volume  is: 

F  =  J  X  Bo  (3.1.15) 

Eqs.  (3.1.14)  and  (3.1.15)  are  all  that  we  need  to  evaluate  the  Lorentz  force  in 
low-i?m  MHD. 

There  is  no  need  to  calculate  b  since  it  does  not  appear  in  the  Lorentz  force. 
Moreover,  J  is  uniquely  determined  by  (3.1.14)  since: 

V-J  =  0  (3.1.16) 

V  X  J  =  aV  X  (U  X  Bq)  (3.1.17) 

and  a  vector  held  is  unambiguously  determined  if  its  divergence  and  curl  are 
known. 

From  now  on  we  shall  drop  the  subscript  on  Bq  and  Eq,  on  the  understanding 
that  B  represents  the  imposed,  steady  magnetic  held  and  E  is  the  externally 
applied  electric  held  [15]. 

3.2  Axially  symmetric  MHD  equations 

The  set  of  equations  on  which  the  numerical  model  has  been  built  is: 


V  •  B  =  0 

(3.2.1) 

V  X  E  =  0 

(3.2.2) 

V- J  =  0 

(3.2.3) 
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J  =  a(E  +  UxB)  (3.2.4) 

(9U 

p—  +  p(U  ■  V)U  =  -VPn  +  i^pAV  +  J  X  B  (3.2.5) 

at 

V-U  =  0  (3.2.6) 

where  Eq.  (3.2.1)  states  the  conservation  of  the  magnetic  field  flux  density, 
Eqs.(3.2.2)-(3.2.4)  model  the  electric  side  of  the  problem  while  Eqs. (3.2.5)- 
(3.2.6)  the  fluid-dynamics  side,  [16]. 

The  shape  of  the  previous  set  of  equations  in  the  cylindrical  coordinates  r  and 
X  is  given  in  detailed  way  in  the  following  paragraphs. 

3.2.1  Magnetic  Field 


The  Eq. (3.2.1)  represents  the  conservation  of  the  magnetic  flux  density: 


r  or  oz 


(3.2.7) 


The  low-i?em  approximation  implies  that  the  induced  magnetic  field  generated 
by  the  electric  current  flowing  through  the  liquid  is  neglectable  compared  to 
the  applied  magnetic  held.  In  this  way  the  magnetic  held  B  can  be  considered 
steady  over  the  time  and  it  can  be  solved  independently  from  Eqs.  (3.2.2)- 
(3.2.6). 

As  we  assume  the  device  to  be  characterized  by  a  pure  radial  magnetic  held 
with  an  assigned  value  Bq  at  the  inner  radius  r*,  the  Eq. (3.2.6)  can  be  reduced 
to: 


V.B=id^=0 

r  or 


(3.2.8) 


and  then: 
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rBr  =  const  =  riB^  ^  B  =  — — ^  (3.2.9) 

r 

From  now  we  will  not  specify  anymore  the  direction  of  magnetic  field. 

3.2.2  Navier-Stokes  equations  and  continuity  equation 
for  the  velocity 

In  this  paragraph  Eqs. (3.2.5)  and  (3.2.6)  will  be  simplihed  using  the  hypothesis 
of  laminar  flow  and  axial  symmetry. 

The  Eq.  (3.2.6)  is  the  continuity  equation  for  the  velocity  and  it  can  be  rewritten 
as: 


V-U 


ld{rUr)  ^  dU, 
r  dr  dz 


(3.2.10) 


in  which  the  azimuthal  component  has  been  neglected  due  to  the  hypothesis 
of  axial  symmetry. 

Limiting  the  analysis  to  the  the  laminar  case,  we  can  consider  the  vorticity 
equation  along  the  direction  9: 


V  X  U 


dUr 

dz 


d^ 

dr 


0 


(3.2.11) 


and  introducing  the  potential  of  the  velocity  <h  we  can  say  that  Eq.(3.2.11)  is 
satished  everywhere  in  the  plane  r  —  zif: 


( 


< 


\Jr 

u. 


(9$ 

dr 

d^ 

dz 


(3.2.12) 


Substituing  Eq.(3.2.12)  into  Eq.(3.2.10)  we  obtain; 


^  I  d  [  dcp 


H - =  V^<h 

dz^ 


(3.2.13) 


17 


DISTRIBUTION  A.  Approved  for  public  release:  distribution  unlimited. 


Cap.  3  Assumptions:  Axially-symmetric  MHD  problem  and  Low  Rem  approximation. 


Eq.(3.2.13)  is  the  Laplace  equation  in  the  space  r  —  z  and  it  is  well-known  that 
if  there  is  no  source  term  and  all  the  boundary  conditions  are  Neumann  ho¬ 
mogeneous  boundaries  conditions  the  solution  is  constant  all  over  the  domain. 
We  can  then  say  that  the  radial  and  axial  components  of  the  velocity  are  zero 
under  the  hypothesis  of  laminar  flow. 

Thanks  to  this  conclusion,  Eq.  (3.2.5)  can  be  reduced  to  the  only  azimuthal 
component: 


dUe 

'  dt 


= 


r  dr 


dr 


+ 


d^Uo  Uf> 


dz"^ 


J zBj, 


(3.2.14) 


in  which  the  hypothesis  of  laminar  flow  has  been  applied  and  pressure  gradient 
along  the  azimuthal  coordinate  has  been  neglected  [17]  [18]. 


3.2.3  Electric  current  density  and  Back  electromotive 
force 

The  electric  side  of  the  problem  is  classically  solved  reducing  the  Eqs.  (3.2.2)- 
(3.2.4)  to  a  single  eqnation  in  the  electrostatic  potential.  The  solntion  of  this 
eqnation  gives  information  abont  the  electric  cnrrent  necessary  to  integrate 
the  Navier-Stokes  eqnation.  This  model  will  be  solved  with  a  classical  Finite 
difference  algorithm  and  for  this  reason  will  be  called  FDM,  Finite  Difference 
Model. 

The  Eq. (3.2.2),  also  called  Faraday  Law,  states  that  a  variable  magnetic  field 
generates  a  not  null  curl  of  the  electric  field.  In  the  case  in  study  the  mag¬ 
netic  held  is  constant  over  the  time  and  the  electric  held  will  be  considered 
irrotational.  Moreover  the  low-Rcm  approximation  implies  that  the  indnced 
magnetic  held  generated  by  the  electric  current  howing  through  the  liquid  is 
neglectable  compared  to  the  applied  magnetic  held.  We  then  have 


V  X  E  =  0 


dEr  dE, 


dz  dr 


=  0 


(3.2.15) 


In  order  to  satisfy  the  (3.2.15)  we  introdnce  the  potential  of  the  electric  held 
(f)  dehned  as: 
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( 


< 


Ej. 

E, 


d4> 

dr 

d(j) 

dz 


(3.2.16) 


In  this  way  Eq.(3.2.15)  is  solved  in  all  the  r  —  z  domain. 

The  Eq.  (3.2.3)  describes  the  conservation  of  electric  current;  as  for  the  Eq.  (3.2.2) 
it  can  be  expressed  in  cylindrical  coordinates  as: 


V- J 


Jr  dJr  dJz 
El  ^ - 1  - 1 

r  dr  dz 


(3.2.17) 


Finally  Eq.  (3.2.4)  also  called  Ohm’s  law  for  moving  conductors  links  the  veloc¬ 
ity  held  with  the  magnetic  held,  electric  current  and  electric  held.  We  rewrite 
its  components  along  r  and  and  using  Eq.(3.2.16)  we  have: 


Jr  =  cr  {Er  -|-  (U  X  B)  ■  f )  = 

7? 

=  aEr  =  -CT— 

dr 


(3.2.18) 


Jz  =  O'  {Ez  (U  X  B)  ■  5)  = 

=  a{Ez-UeB)  =  -a(^^+UeB 


(3.2.19) 


It  is  convenient  to  assemble  Eqs.(3.2.16),  (3.2.17),  (3.2.18)  and  (3.2.19)  in  order 
to  have  all  the  electric  quantities  function  of  the  only  electrostatic  potential. 
We  compute  the  derivatives  with  respect  of  r  and  x  of  the  homologous  compo¬ 
nents  of  the  electric  current  density: 


dJr  dJ(j) 

dr  dr"^ 


dJz  ( dUg\ 


(3.2.20) 

(3.2.21) 
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Substituting  these  results  in  Eq.(3.2.17)  we  obtain  the  equation  in  the  electro¬ 
static  potential: 


1(90  (9^0  (9^0  ^ 

0^  +  ^  +  ^  =  “^  ) 

An  alternative  method  to  investigate  the  electric  current  distribution  is  to 
apply  the  node- voltage  analysis  to  an  electric  network  consistent  with  the  mesh- 
grid  used  for  solving  the  Navier-Stokes  equation.  Thanks  to  low-i?em,  we  can 
say  that  the  electric  network  will  be  only  made  by  electric  resistances  and 
voltage  generators  representing  the  back  electromotive  voltage  induced  by  the 
spinning  liquid  through  the  magnetic  held.  This  approach  will  be  described 
in  the  discussion  of  the  numerical  method  used  to  solve  the  problem  with  the 
name  FDHM,  Finite  Differences  Hybrid  Model  because  the  Finite  differences 
approach  will  be  used  only  for  the  huid-dynamics  part  of  the  problem  while  the 
algebraic  node- voltage  analysis  will  be  used  for  the  electric  part  of  the  problem. 

3.2.4  Boundary  and  Initial  Conditions 

The  complete  set  of  equations  for  the  low-i?^  MHD  has  been  reduced  to  a 
system  of  two  PDF  equations,  Fq.(3.2.14)  and  Fq.(3.2.22),  in  the  unknowns 
Ug,  velocity  of  the  huid,  and  0,  electrostatic  potential.  A  further  relation  is 
needed  to  link  the  two  PDE  equations;  the  relation  chosen,  Eq.(3.2.19),  pro¬ 
vides  the  conversion  of  the  electrostatic  potential  to  the  2;-component  of  the 
electric  current  density  needed  to  compute  the  azimuthal  velocity  held. 

In  order  to  obtain  a  solntion  for  the  problem  we  have  to  set  the  right  bonndary 
conditions  for  all  the  nnknown  taken  into  account. 

The  liqnid  is  enclosed  in  a  torus  with  rectangular  cross  section.  Neglecting  the 
effect  of  the  capillarity  in  the  zones  next  to  the  corners  and  the  effect  that  a 
not  wetting  liquid  could  have  on  the  velocity  held,  we  set  the  no-slip  condition 
on  all  the  boundaries  of  the  cross  section.  It  means  that  the  azimnthal  velocity 
of  the  huid  at  the  boundaries  is  set  to  the  same  velocity  of  the  torus,  and  con¬ 
sidering  the  torus  at  rest  and  bounded  to  a  hxed  reference  system  the  velocity 
of  the  hnid  at  the  bonndaries  is  set  to  zero. 
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Ue  {vi,  z,  t)  =  Ug  (re,  t)  =  Ug  (r,  0,  t)  =  Ug  (r,  h,t)  =  0 


(3.2.23) 


The  initial  condition  for  the  velocity  is: 


Ug  {r,z,0)  =  0 


(3.2.24) 


d(f)  {vi ,  z) 
dr 


H  r,  h)  =  f{t) 


=  0 


r,  0)  =  0 


d(f)  {re ,  z) 
dr 


b) 


u,(n,z] 


Ue  (r,  h)  =  0 


0  Ue{re,z)=0 


f/»(r.0)  =  0  c) 


Figure  3.1:  Boundary  Condition  applied  to  the  Finite  Difference  model. 

a)  Reference  system  and  dimension  of  the  cross  section 

b)  Electric  boundary  conditions 

c)  Fluid- dynamics  boundary  conditions 


On  the  electric  side  we  shall  impose  boundary  conditions  allowing  the  electric 
current  to  flow  through  the  liquid.  We  choose  an  external  electric  field  along  z 
and  then  we  consider  a  net  flux  of  electric  current  only  through  the  boundaries 
parallel  to  the  r  axis  while  insulating  the  boundaries  parallel  to  the  axis. 
The  external  electric  field  is  generated  by  time  varying  electrostatic  potential 
at  the  upper  boundary,  z  =  h  and  we  consider  a  null,  constant  over  the  time. 
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reference  electrostatic  potential  at  the  lower  boundary,  z  =  0. 


o 

o 

(3.2.25) 

0  {r,  h)  =  f  it) 

(3.2.26) 

The  condition  of  isolated  wall  for  the  boundaries  parallel  to  the  axis  can 
be  assigned  with  the  use  of  Eq.(3.2.18),  putting  the  radial  component  of  the 
velocity  equal  to  zero: 

T  -  -n 

Jr  ^  O  ^ 

or 

(3.2.27) 

and  obtaining: 

d(j){ri,z)  _ 
dr 

(3.2.28) 

d(j)  (re,  ^)  _  g 
dr 

(3.2.29) 
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Chapter 


Methods  and  Procedures:  Axially- 
symmetric  Numeric  Models 


In  this  chapter  a  classical  2-dimensional  Finite  Differences  model  (FDM) 
will  be  described  and  applied.  A  second,  alternative  method,  called  Finite 
Difference  Hybrid  Model  (FDHM),  is  developed  on  the  base  of  the  MHD  set 
equation.  The  FDHM  is  then  tested  in  order  to  highlight  its  conservative 
capabilities. 
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4.1  FDM:  Finite  Differences  Model 


The  FDM  has  been  solved  numerically  using  a  consistent  and  conservative 
scheme  over  a  collocated  not  uniform  mesh  [19]  [20]  applied  to  the  following 
set  of  equations: 


P- 


dUo 


dt 


=  P 


Id  f  dUe\  d^Ue  Ue 

dz"^  ^2 


r  dr  \  dr 


+  Jz^ 


1  dcf)  d'^(j)  dUe 

r  dr'^  dz'^  dz 

j  _  _  ^ 

dr  ^ 

or 

Jz  =  -a{^  +  UeB 


dz 


Bl'r'l  = 


TiB, 


i^O 


(4.1.1) 

(4.1.2) 

(4.1.3) 

(4.1.4) 

(4.1.5) 


In  this  collocated  grid  the  velocity,  Ug,  and  electrostatic  potential,  0,  are  cal¬ 
culated  at  the  center  cell  while  the  radial  and  axial  electric  current  densities 
Jr  and  Jz  are  evaluated  on  the  cell  faces  respectively  parallel  to  the  2;  and  r 
axes.  Figure  4.1  . 


Figure  4.1:  Arrangment  of  the  physical  quantities  in  a  collocated  grid. 

The  algorithm  starts  from  the  calculation  of  the  electrostatic  potential  at  the 
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center  cell  considering  a  null  initial  condition  for  the  velocity  field  Ug]  once  the 
electrostatic  potential  is  known  in  all  the  domain,  it  is  possible  to  calculate 
the  current  densities  by  means  of  Eqs.(4.1.3)-(4.1.4).  The  calculated  electric 
current  densities  refer  to  the  cell  faces  and  in  order  to  proceed  to  the  calculation 
of  velocity  held  they  need  to  be  interpolated  to  the  center  cell  in  order  to  have 
the  Lorentz  force  at  the  center  cell.  The  mono-dimensional  and  time-dependant 
Navier-Stokes  equation  is  then  solved  used  the  Crank-Nicolson  method  over  a 
not  uniform  grid  over  the  space  r  —  z  and  a  uniform  time  sequence. 

In  the  following  the  discrete  formulas  for  Eqs.(4.1.1)-(4.1.4)  will  be  given  in 
details. 

4.1.1  Conservative  and  Consistent  Laplace  and  diver¬ 
gence  operator  in  cylindrical  coordinates 

In  the  collocated  grid,  the  conservation  of  the  electric  current  in  differential 
form  refers  to  the  cell  faces: 


V- J 


Jr  dJr  dJz 

- ^  - ^ 

r  dr  dz 


(4.1.6) 


The  terms  of  Eq.(4.1.6)  can  be  expressed  in  the  discrete  form  as  follows: 


Jr  _  ^  i’^r)k+l/2,l  +  i’^r)  k-ll2,l 
r  Tk  2 

dJr  _  i'Jr)k+l/2,l  ~  i'^r)k-l/2,l 
dr  rk+i/2  -  rk-i/2 

dJz  _  iJz)k^l+l/2  ~  ('^2)fc,«-l/2 
dz  ^1+1/2  —  Zl-i/2 


(4.1.7) 

(4.1.8) 

(4.1.9) 


The  contribution  of  the  electrostatic  potential  to  the  electric  current  J®  can  be 
separated  from  the  one  generated  by  the  electric  induction 


J  = -V0  +  U  X  B  =  r  + J"  (4.1.10) 
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Figure  4.2:  Finite  Differences  Model  stencil. 


The  components  of  the  conductive  part  of  electric  current  density  at  the  cell 
faces  can  be  related  with  the  values  of  the  electrostatic  potential  at  center  cell. 


('^r)fc+l/2,Z  ~ 

i'^z)k+l/2,l  ~ 
iJz)k-l/2,l  ~ 


0fc+l,Z 

-  4’k,l 

Tk+l 

-  Tk 

4>k,l  - 

4>k-l,l 

Tk  - 

Tk-l 

0A:,Z+1 

-  4>k,l 

^Z+1 

-  Zi 

4>k,l  — 

4>k,l-l 

zi  - 

Zl-1 

(4.1.11) 


(4.1.12) 


Replacing  Eqs.(4.1.11)-(4.1.12)  into  Eqs.(4.1.7)-(4.1.9),  for  the  Eq.(4.1.7)  we 
have: 


r 


Oil  0fc+l,z  +  012  4>k,l  +  Ol3  4>k-l,l 


(4.1.13) 
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1 


ail  — 


2rfc  \  rfc+i  -  Tfc 

1/1  1 


^12  — 


2rfc  \rfc+i  -  Tfc  Tfc  -  Tfc-i 
1  /  1 

Ol3  — 


2rfc 


(4.1.14) 


for  the  Eq.(4.1.8): 


dJ^ 

dr 


—  0,21  4>k+l,l  +  ®22  <Pk,l  +  ^23  <Pk-l,l 


(4.1.15) 


021  — 


0,22  —  — 


r fc+1/2  —  rk-1/2/  k+1  ~ 

I  \  f  I  1 


,^^fc+l/2  -  rfc-1/2/  V  Tk+l  -  Tk  Tj  -  Tj-i 

a..  ^  ^ - )  (^ 

V^fc+1/2  - ’^fc-1/2/  \rk-rk-i 


(4.1.16) 


for  the  Eq.(4.1.9): 


dz 


—  O3I  4>k,l+l  +  O32  4>k,l  +  033  (j)k,l-l 


(4.1.17) 


031  — 


032  —  — 


Zl+1/2  —  Zl-i/2  J  V^«+l  ~ 

1  \  /  1  1 


^i+1/2  —  Zl-112 


O33  —  — 


Zi+i  -  Zi  Zi-  Zi-I 

1 


Zl+1/2  —  Zi-i/2  J  \Zl  —  Zi-i 


(4.1.18) 


Considering  that  the  velocity  U  is  purely  azimuthal  while  the  magnetic  held 
B  is  purely  radial,  the  inductive  part  of  the  Eq.(4.1.10)  can  be  written  as: 
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f  9  z 

f  9  z 

J“  =  U  X  B  = 

Ur  Ug  U, 

= 

0 

0 

Br  Bg  Bz 

5,.  0  0 

{BrUg)z  (4.1.19) 


The  divergence  of  J“  has  only  the  component  parallel  to  the  axis  and  using 
the  values  of  and  Ug  at  the  cell  faces  we  have: 


V  - J"  = 


dJ^  {J") 


dz 


k, 1+1/2 


-  W 


^  ik, 1-1/2 


Zl+1/2  —  Zl-i/2 

~  i^r)k  i^9)k, 1+1/2 


+  i^r)k  i^e)k, 1+1/2 


(4.1.20) 


=  -C 


Zl+l/2  —  Zl-i/2 


Finally  the  equation  in  the  electrostatic  potential  at  the  cell  center  can  be 
expressed  as: 


(Oii  +  a2l)  0fc+l  z+  (ai2  +  022  +  O32)  (t>k,l  +  (Ol3  +  ^23)  0fc-l,Z+  ,  ^  , 

(4.1.21) 

®31  0fc,Z+l  +  O33  4>k,l-l  =  C 

4.1.2  Interpolation  Operator  for  not  equispaced  grid 

The  Navier-Stokes  equation  gives  values  of  Ug  at  the  center  cell  but  in  order  to 
compute  the  inductive  part  of  the  Eq.(4.1.20)  values  of  Ug  at  the  cell  faces  are 
needed.  On  a  collocated  grid  system,  an  interpolation  operator  Ac^./  can  be 
dehned,  which  represents  the  interpolation  from  the  cell  center  to  cell  faces: 

U/  =  Ac^./  Uc  (4.1.22) 

where  U  denotes  the  velocity  vector  {Ur,Ug,Uz).  The  detailed  formula  of  the 
interpolation  operator  on  collocated  grid  is: 
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Ufc+i/2,«  —  Ac_s./r+  Uc  —  +  (1  ~  (4.1.23) 

Ufc,z+i/2  =  Ac_s./z+  Uc  =  /3+Ufc^i  +  (1  ~  l3+)  Ufc,i+i  (4.1.24) 

The  subscript  c  — )■  /r+  and  c  — )■  fz-\-  denotes  the  interpolation  from  the  cell 

center  normal  to  the  r  and  2:  axes;  the  symbol  +  denotes  that  the  cell  faces  are 

located  to  the  right  and  top  of  the  cell  center. 

The  coefficients  of  the  interpolation  operator  Ac_s./  are  dehned  on  the  base  of 
the  grid  arrangement: 

^  rfc+i  -  rk+1/2  ^  zi+i  -  zi+i/2  (4.1.25) 

rfc+i  -  Tk  zi+i  -  zi 

An  equivalent  operator  can  be  dehned  for  the  interpolation  from  the  cell  faces 
to  center  cell  and  it  can  be  used  to  obtain  the  electric  current  density  at  center 
cell  once  known  it’s  value  at  the  cell  faces.  Due  to  the  position  of  the  center 
cell,  considered  to  be  in  the  middle  between  the  two  faces  of  the  cell,  we  have: 

Jc  =  Af^cJf  (4.1.26) 


where  Aj_^c  has  the  meaning  of: 


(Jr)  li,t  —  {Jr) Jr 


(Jr) 


J  k,l 


w  k,l  ~  {Jz) fz 


('^^)fc,i+l/2  +  1-1/2 


k,l 


(4.1.27) 

(4.1.28) 


4.1.3  Grid  Generator 

The  mathematical  model  is  solved  on  a  not  equispaced  grid  in  order  to  study 
with  accuracy  the  boundary  effects  due  to  the  coupling  between  the  electric 
and  huid-dynamic  problems.  It  is  possible  to  dehne  two  different  characteristic 
lengths  for  the  boundary  layers,  [21]  [22]  [23]: 
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•  for  the  boundaries  normal  to  the  magnetic  held,  parallel  to  the  axis, 
the  boundary  layers  are  called  Hartmann  layers  and  their  lengths  are: 

SHrtm  =  (4.1.29) 

•  for  the  boundaries  parallel  to  the  magnetic  held,  i.e.  parallel  to  the  r 
axis,  the  boundary  layers  are  called  Sherclih  layers  and  their  lengths  are: 

where  L  is  the  diherence  between  the  inner  and  outer  radius  of  the  cavity 
containing  the  liquid. 

With  the  purpose  to  make  several  simulations  for  diherent  values  of  magnetic 
held  B,  electric  conductivity  of  the  huid  a,  length  of  the  magnetic  gap  L  and 
height  of  liquid  torus  h  it  has  been  necessary  to  understand  how  to  modify  the 
computational  grid  in  function  of  the  dimension  of  Surtm  and  5sh- 
The  grid  generator  create  four  diherent  kinds  of  grid.  For  the  r  axis: 

•  when  hurtm  <  ~  h)  /6  the  grid  on  the  r  axis  is  characterized  by  an 

equispaced  mesh  in  the  zone: 


Vi  <r  <ri  +  25Hrtm  re-25Hrtm<r  <re  (4.1.31) 

while  in  the  zone  r*  +  25Hrtm  <  r  <  r*  +  25Hrtm  the  grid  is  logarithmic 

•  if  Surtm  >  ij'e  —  Ti)  /Q  the  grid  is  logarithmic  from  Vi  to  rg. 

For  the  z  axis: 

•  when  Ssh  <  h/6  the  grid  on  the  axis  is  characterized  by  an  equispaced 
mesh  in  the  zone: 


0  <  z  <  25sh  h  —  25sh  <  z  <  h  (4.1.32) 

while  in  the  zone  25sh  <  ^  <  h  +  25sh  the  grid  is  logarithmic 
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•  if  6sh  >  h/6  the  grid  is  logarithmic  from  0  to  h. 


Figure  4.3:  Grid  for  dsh  <  h/6  and  durtm  <  {r^  -  r,)  /6. 


Figure  4.4:  Grid  for  Ssh  >  h/6  and  dHrtm  >  (re  —  Ti)  /6. 


4.1.4  Navier-Stokes  Equation 

The  Navier-Stokes  equation  Eq.(4.1.1): 


P 


dUe 

dt 


P 


1  dUe  d'^Ue 
- -  H - -  + 

T' 


d^Ue 

dz^ 


+  Jz^ 


(4.1.33) 
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is  solved  using  a  Crank-Nicolson  scheme  on  a  five  points  stencil,  (Appendix  A). 
The  finite  difference  scheme  has  been  calculated  for  a  not  equispaced  mesh.  The 
time  derivative  has  been  evaluated  with  a  second  order  finite  difference  scheme 
with  constant  time  step.  The  shape  of  the  discretized  version  of  the  derivatives 
present  in  Eq.(4.1.33)  is: 


IdUn 


r  dr 


bi,i  Uk+i,i  +  &i,2  Uk,i  +  &i,3  Uk-i,i  (4.1.34) 


fell  =  — 


fel2  —  - 


rk,i  \rk+i,i  —  ^k,i, 

1/1  1 

- + - 


\  fk+l,l~'^k,l  f'k,l  —  f'k-1,1 
rk,i  \rk,i  -  rk-i,i 


(4.1.35) 


d^Ue 

dr‘i 


fe2,l  Uk+l,l  +  fe2,2  Uk,l  +  fe2,3  Uk-l,l  (4.1.36) 


fe21  — 


fe22  — 


k+l,l  ~  ^ k,l 

1  1 


rk+i,i  ~  f'k,i  Tk^i  —  rk-i,i 

fe23  =  - 

'f'kd  —  f'k-ld 


(4.1.37) 


d^Ue 

dz'^ 


fes,!  Uk,l+1  +  fe3,2  Uk,l  +  fe3,3  Uk,l-1  (4.1.38) 


fe31  — 


fe32  — 


^A:,Z+1  —  Zk^l 

1  1 


^k^l+1  ^k,l  ^k^l—1 

Zk,l  —  Zk,l-1 


(4.1.39) 
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Ue  Uk,i  _  ,  jj 

~  ~^4,2  Uk,l 

^  n,i 


(4.1.40) 


dUe  ^  t/r;'  -  Ul, 


dt 


-)■ 


At 


(4.1.41) 


The  right  side  of  the  discrete  version  of  Eq.(4.1.33)  has  then  the  following 
shape: 


F  {u,  J^,  B,  r,z)  =  n  [(6i,i  +  62,1)  Uk+i,i 

+  (^1,2  +  ^2,2  +  ^3,2  +  ^4,2)  Uk^l 

+  (&i,3  +  M  (4-1-42) 

+  &3,1  Uk,l+1  +  &3,3  +  {Jz)k^l  Bk,l 


Applying  the  Crank-Nicolson  algorithm  we  obtain: 


_  'll^  1 

=  Yp  Z,B,T,z)  +  F'(u,.].,B,r,z)]  (4.1.43) 

The  (4.1.43)  leads  to  a  tridiagonal  system  that  can  be  easily  solved  with  tridi¬ 
agonal  matrix  algorithms. 

The  computation  of  the  velocity  Ug  extends  to  the  all  the  center  cell  nodes. 
The  no-slip  condition  is  implemented  adding  node  on  the  boundaries,  where 
Ug  is  known  and  equal  to  zero,  in  {k  —  1/2,  /)  and  {k  -|-  1/2,  /)  for  the  bound¬ 
aries  parallel  to  the  axis  and  nodes  in  {k,l  —  1/2)  and  {k,l  -\- 1/2)  for  the 
boundaries  parallel  to  the  r. 
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4.1.5  Discrete  boundary  conditions 


In  order  to  solve  the  electrodynamic  problem  (3.2.22)  the  boundary  conditions 
(3.2.25) (3.2.26)  and  (3.2.28)-(3.2.29)  shall  be  applied. 

The  Dirichlet  conditions  on  the  top  and  bottom  side  of  the  cross  section  are 
simple  to  express.  They  are  known  values  of  the  electrostatic  potential  and 
they  appear  inside  the  vector  of  known  terms  of  the  algebraic  system. 

The  Neumann  conditions  on  the  right  and  left  side  of  the  cross  section  needs  to 
be  studied.  The  classical  way  to  implement  the  conditions  on  the  derivative  is 
using  ghost  nodes  on  the  grid  coupling  this  conditions  with  the  general  equation 
for  the  electrostatic  potential. 

The  ghost  nodes  are  virtual  nodes  external  to  the  computational  grid  for  which 
is  not  supposed  to  know  the  value  of  the  unknown. 

For  the  left  side  boundary  we  have  that  the  Eq. (3.2.28)  can  be  expressed  as: 


dej)  {ri,z)  _  g  4>k+i,i  ~  <Pk,i  <Pk,i  —  0fc-i,z 

dr  rfc+i  -Tk  Tk-  Tk-i 


(4.1.44) 


On  this  side  of  the  cross  section,  the  node  {k  —  1,  /)  is  a  ghost  node  and  we  can 
highlight  the  electrostatic  potential  on  this  node  in  order  to  obtain: 


(j>k-l,l 


4>k+l,l 


( Tk  -rk-i\ 
\rk+i  -Tk) 


+  (j>k,l 


r-k  -  rk-i\ 

rk+i  -  Tk) 


(4.1.45) 


Introducing  Eq.(4.1.45)  into  Eq. (3.2.22)  we  obtain  the  Neumann  condition  ap¬ 
plied  using  central  hnite  differences  scheme. 

In  the  same  way  for  the  right  boundary  of  the  cross  section,  we  have: 

(90  (fe,  Z)  _  g  0A:+1,Z  "  0A:,Z  4’k,l  —  0fc-l,Z  _  g  ^ 

dr  rk+i  -rk  rk  -  rk-i 


which  means: 


4>k+l,l 


4>k,l 


rk+1  -rk\ 
rk  -  rk-ij 


+  0fc-l,Z 


/ rfc+i  -rk\ 
\rk-  rk-ij 


(4.1.47) 
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With  the  same  procedure  we  obtain  the  Neumann  condition  on  the  right  bound¬ 
ary  of  the  cross  section. 

4.1.6  Voltage  Drive  and  Current  Drive 

The  model  can  simulate  the  reaction  of  the  conductive  liquid  to  a  certain  kind 
of  law  f{t)  for  the  electrostatic  potential  assigned  between  the  lower  and  the 
upper  boundaries  of  the  cross  section. 

It  does  not  seem  to  be  able  to  simulate  an  assigned  trend  g(t)  of  the  electric 
current;  it  is  not  possible  to  assign  as  boundary  condition  a  given  law  for  the 
electric  current  flowing  through  the  conductive  liquid. 

In  order  to  make  it  possible  we  have  to  consider  the  Eq.(3.1.14).  Being  the 
external  electric  held  along  we  can  replace  the  divergence  with  the  total 
derivative  of  the  electrostatic  potential  with  respect  to  2;.  We  have: 

J  =  -a^  +  aUxBo  (4.1.48) 

oz 

Integrating  over  the  surface  in  the  plane  r  —  6  we  obtain  the  total  current 
howing  in  the  conductive  liquid: 

1  = -a  f  ^dS +  a  [  V  xBodS  (4.1.49) 

Js  dz  Js 

In  order  to  simplify  the  derivative  of  the  electrostatic  potential  with  respect  of 
2;  component  we  proceed  integrating  along  2;: 


dS  dz  -\-  a 


V  X  Bn  dS  dz 


(4.1.50) 


Due  to  the  conservation  of  electric  current  and  being  the  electrostatic  potential 
just  function  of  the  2:  we  have: 

O  ^  ph  p 

I  =  — -A0  -  -  /  V  X  Bo  dS  dz  (4.1.51) 

^  ^  Jo  Js 


35 


DISTRIBUTION  A.  Approved  for  public  release:  distribution  unlimited. 


Cap.  4  Methods  and  Procedures:  Axially- symmetric  Numeric  Models 


What  we  have  found  is  the  relation  between  the  total  electric  current  and  the 
voltage  applied  to  the  upper  boundary  taking  into  account  the  back  electro¬ 
motive  force.  We  can  say  that  we  can  numerically  impose  a  certain  electric 
current  of  known  law  I {t)  through  the  following  relation: 

(j){r,h,t)  =  — ~  4  /  [  U  X  30  dS  dz  (4.1.52) 

Jo  Js 

where  with  ^(r,  h,  t)  we  intend  the  electrostatic  potential  at  the  upper  boundary 
being  zero  its  value  at  the  lower  boundary. 

4.1.7  Iterative  scheme 

The  solving  algorithm  uses  in  sequence  all  the  relations  presented  in  previous 
section.  In  order  to  have  a  clearer  vision  of  the  algorithm  a  flux  diagram  is 
presented,  Figure  4.5. 
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Figure  4.5:  Flux  Diagram  of  the  FDM  model 


4.1.8  Conservation  problems 

The  FDM  shows  problem  of  conservation  of  electric  cnrrent.  Fignre  (4.6)  shows 
the  total  cnrrent  flowing  throngh  the  cross  section  along  the  axial  direction 
It  is  evident  that  the  cnrrent  shonld  be  constant  for  each  valne  of  ^  becanse 
of  the  conservation  of  the  electric  cnrrent.  The  algorithm  shows  conservation 
errors  of  electric  cnrrent  near  the  armatures  of  the  device.  This  is  not  accept¬ 
able  because  this  conservation  errors  lead  to  an  underestimation  of  the  Lorentz 
force  acting  on  the  conductive  liquid. 

We  then  try  to  solve  this  problem  with  an  alternative  method.  This  method 
will  be  discussed  in  the  next  sections. 
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FDM  —  Axial  current  —  ^{z) 
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Figure  4.6:  FDM  -  Total  axial  electric  current 


The  radial  component  of  the  electric  current,  Figure(4.7)  has  no  problem  and 
it  shows  a  perfect  antisymmetric  trend  with  null  mean. 


FDM  —  Radial  current  —  Fiz) 


Figure  4.7:  FDM  -  Total  radial  electric  current 
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4.2  FDHM:  Finite  Differences  Hybrid  Model 

The  Finite  Difference  Hybrid  Model  is  an  alternative  algorithm  to  solve  the 
MHD  equations.  The  difference  with  the  FDM  is  in  the  way  the  electric  part 
of  the  model  is  solved.  Instead  of  solving  the  electrostatic  potential  equation,  an 
electric  network  is  created  in  order  to  obtain  the  distribution  of  the  electrostatic 
potential  and  of  the  electric  current  inside  the  cross  section  of  the  device. 
The  electric  network  is  then  solved  applying  the  nodal  analysis.  This  model 
has  been  inspired  by  [24],  anyway  this  model  solves  numerically  the  Navier- 
Stokes  equations  instead  of  imposing  a  hxed  velocity  prohle  and  it  is  set  up 
to  compute  the  transient  between  the  state  of  null  velocity  to  the  state  of 
completely  developed  MHD  laminar  flow.  The  same  grid  generator  of  the 
FDM  is  used  in  the  FDHM.  The  cell  centers  become  here  the  nodes  of  electric 
network  while  the  cell  faces  are  the  branches. 

4.2.1  Equivalent  electric  network 

The  equivalent  electric  network  is  made  of  radial  and  axial  resistances  and 
voltage  generators. 

The  radial  resistances  link  the  nodes  in  the  radial  direction.  Applying  the 
conservation  of  the  electric  current  along  the  radius  on  a  sequence  of  cylindrical 
surface,  a  law  for  the  density  of  electric  current  can  be  found: 


Jr)  ^ 

dr  ^  =  (4.2.1) 

Jr  (^l)  Jq 


Introducing  first  the  laws  linking  the  electric  current  to  the  electric  current 
density  and  the  voltage  to  the  electric  held: 


I 


Jr  (r)  dS 


2n 


JrVdzdO  =  271  Jovl 


0  JO 


(4.2.2) 


A0 


Erdr 


'rk 


1 

a 


^’■'=+1  Jo  Vi 


dr 


V-fc 


Jq  r-j 

a 


In 


(4.2.3) 
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© 


Figure  4.8:  Equivalent  electric  network. 

a)  Nodes  and  cells 

b)  Electric  radial  resistances 

c)  Electric  axial  resistances  and  back  electromotive  voltage 
b)  Assembled  electric  network 


and  then  the  Ohm  law,  we  obtain: 


In  f 

r.  _  ^  _  J_  VO+1/ 

"  I  27rcT  I 

The  axial  resistances  link  the  nodes  in  direction  parallel  to 

A0  1  I 

I  a  TT  (r2  -  rl^,) 


(4.2.4) 


(4.2.5) 
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The  voltage  generators  model  the  back  electromotive  voltage  generated  by  the 
conductive  liquid  spinning  in  the  magnetic  held.  Their  values  can  be  computed 
integrated  the  cross  product  between  velocity  and  magnetic  held  over  the  length 
parallel  to  ;2.  Using  the  Ohm  Law  for  moving  conductor  we  have: 

pSzi 

^(pind  =  u  X  Bdz  =  -UqB  (r)  5zi  (4.2.6) 

Jo 


4.2.2  Boundary  Conditions 

The  FDHM  does  not  need  specihc  equations  to  impose  the  boundary  condi¬ 
tions.  The  boundary  conditions  on  the  side  of  the  cross  section  parallel  to  the 
X  axis  are  imposed  not  taking  into  account  any  radial  resistances  that  could 
bring  electric  current  in  the  direction  of  the  boundary.  So,  the  nodes  adjacent 
to  vertical  side  of  the  cross  section  have  just  three  branches:  two  axial  branches 
and  one  radial  branch  towards  the  inside  of  the  cross  section. 

The  boundary  conditions  on  the  side  of  the  cross  section  parallel  to  the  r  axis 
are  applied  imposing  the  same  value  of  the  electrostatic  potential  to  all  the 
nodes  on  the  boundaries. 


4.2.3  Voltage  and  Current  drive 

In  the  FDHM  the  voltage  and  current  drive  are  realized  with  two  different 
conditions  on  the  nodes  on  the  upper  and  lower  boundaries. 

In  order  to  simulate  the  voltage  drive  the  voltage  on  the  boundaries  is  assumed 
to  be  known.  The  lower  armature  is  set  to  a  reference  voltage  while  the  nodes 
on  the  the  upper  armature  are  set  to  same  value,  in  general,  function  of  the 
time. 

In  the  current  drive  it  is  necessary  to  control  the  current  flowing  through  the 
conductive  liquid  and  to  have  an  equipotential  surface  on  each  armature.  In 
order  to  satisfy  these  two  constrains  the  nodes  on  each  armature  have  been 
collapsed  to  only  one  node  where  to  impose  a  net  electric  current. 
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4.2.4  Iterative  scheme 

The  algorithm  to  solve  the  FDHM  is  presented  in  the  following  flnx  diagram. 


Figure  4.9:  Flux  Diagram  of  the  FDHM  model 


4.2.5  FDHM  performances 

The  FDHM  appears  to  solve  the  conservation  problem  of  the  FDM  algorithm, 
Fignres  (4. 10)  (4. 11).  The  total  axial  electric  current  is  now  constant  all  over 
the  cross  section  and  equal  to  the  supply  value. 
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FDHAl  —  Axial  current  —  I^{z) 


z  —  [mm] 


Figure  4.10:  FDHM  -  Total  axial  electric  current 

The  radial  component  of  the  electric  current  shows,  as  in  FDM,  an  antisym¬ 
metric  trend  with  null  mean. 


FDHM  —  Radial  current  —  I^{z) 


Figure  4.11:  FDHM  -  Total  radial  electric  current 


The  relative  errors  of  the  electric  and  fluid-dynamic  part  can  be  computed  and 
they  are  respectively  of  the  order  of  10“^°  and  10“®. 


43 


DISTRIBUTION  A.  Approved  for  public  release:  distribution  unlimited. 


Cap.  4  Methods  and  Procedures:  Axially-symmetric  Numeric  Models 


The  error  of  the  electric  problem  is  computed  as  the  standard  deviation  of 
the  difference  between  the  commanded  electric  current  and  computed  electric 
current  flowing  through  the  cross  section. 

Err,i  =  a[F{z)]  (4.2.7) 


Figure  4.12:  FDHM  -  Error  of  the  electric  side  of  the  MHD  problem 


FDHM  —  DynamicsRelativeErvor 


Figure  4.13:  FDHM  -  Error  of  the  dynamic  side  of  the  MHD  problem 
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The  error  of  the  dynamic  problem  is  computed  as  the  difference  between  the 
absolute  value  of  the  torque  computed  with  the  numerical  velocity  distribution 
and  the  difference  between  the  viscous  shear  stress  moment  and  the  moment 
of  Lorentz  force: 


Err-dyn  = 


r|  -  \Mlo,  - 
iri 


(4.2.8) 


4.3  Conclusions 

In  this  chapter  two  different  numerical  methods  have  been  analysed  in  order 
to  simulate  the  behaviour  of  an  electrically  conductive  liquid  spinning  in  an 
axial  symmetric  torus  under  the  interaction  of  a  magnetic  held  and  an  electric 
current.  The  hrst  method,  called  FDM,  shows  problems  of  conservation  of 
the  electric  current  due  to  the  numeric  method  adopted.  For  this  reason  an 
alternative  method,  called  FDHM,  has  been  set  up  and  studied.  The  FDHM 
shows  great  accuracy  in  the  estimation  of  both  the  dynamic  and  electric  side 
of  the  MHD  problem. 
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Results  and  Discussion:  FDHM 


In  this  chapter  several  simulations  based  on  the  FDHM  are  presented  in 
order  to  outline  the  performance  of  the  device  in  study  together  with  the 
performances  of  the  numeric  model  itself. 
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5.1  Electric  Current  and  Velocity  Fields 

The  FDHM  method  developed  and  described  in  the  previous  section  gives  the 
possibility  to  visualize  the  distribution  of  the  electric  current  flowing  through 
the  fluid  and  of  the  the  velocity  of  the  fluid  in  the  cross  section  of  the  device. 
The  axial  electric  current,  Figure(5.1),  appears  to  be  forced  to  the  boundaries 
of  the  cross  section  parallel  to  the  ^  axis.  The  reason  of  this  behaviour  is  the 
back-electromotive  force,  function  of  the  velocity  held,  that  pulls  the  electric 
current  line  towards  the  zones  where  the  velocity  of  the  huid  is  lower. 


FDHM  —  Axial  current  density  —  Jz{r^  z) 


z  [mm]  j.  _  [uim] 

Figure  5.1:  FDHM  -  Axial  electric  current  density 


For  the  same  reason  the  radial  electric  current,  Figure(5.2),  shows  an  anti¬ 
symmetric  behaviour  in  order  to  allow  the  electric  charges  to  move  from  the 
equipotential  armatures  to  the  boundaries  parallel  to  the  z  axis  where  the  ve¬ 
locity  of  the  huid  and  so  the  back-electromotive  force  is  lower.  The  velocity  held 
computed  by  the  FDHM,  Figure(5.3),  has  a  typical  shape  diherent  from  the 
classical  hydrodynamic  hux.  It  shows  a  strong  gradient  next  to  the  boundaries 
parallel  to  the  ^  axis.  The  reason  of  this  resides  in  the  strong  concentration  of 
the  axial  component  of  the  electric  current  in  those  zones.  The  core  zone  of 
the  hux  shows  a  slight  linear  trend  in  the  radial  direction,  increasing  from  the 
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FDHM  —  Radial  current  density  —  Jr(r,  z) 


Figure  5.2:  FDHM  -  Radial  electric  current  density 

inner  radius  to  the  outer  radius.  This  typical  trend  is  caused  by  the  radial 
magnetic  distribution  that  increase  in  the  direction  from  the  outer  radius  to 
the  inner  radius.  These  two  opposite  trends  of  the  velocity  and  of  the  magnetic 
held  generate  an  almost  constant  back  electromotive  force  acting  on  the  electric 
current. 


FDHM  —  Linear  velocity  —  Ue{r^  z) 


z  —  [mm] 


r  —  [mm] 


Figure  5.3:  FDHM  -  Linear  velocity 
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Figure(5.4)  shows  that  the  angular  velocity  in  the  core  zone  of  the  flux  is 
constant;  this  helps  to  understand  that  the  azimuthal  velocity  of  the  fluid  in 
the  core  region  varies  linearly  with  the  radius.  This  is  consistent  considering 
that  the  magnetic  held  varies  linearly  with  the  inverse  of  the  radius. 

FDHM  —  Angular  velocity  —  uJz{r^  z) 


Figure  5.4:  FDHM  -  Angular  velocity 

5.2  Torque- Current  characteristic 

Due  to  the  presence  of  a  liquid  spinning  the  behaviour  of  this  actuator  appears 
to  be  different  from  the  classical  reaction  wheels. 


X  10"’ 


Time  —  [s] 


Figure  5.5:  FDHM  -  Torque 
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In  particular,  the  viscous  shear  leads  to  a  different  behaviour  of  the  Torque- 
Current  characteristic:  while,  for  the  classical  reaction  wheel,  the  torque  pro¬ 
vided  is  constant  when  the  electric  current  fed  to  the  device  is  constant,  here, 
the  torque  is  constant  for  a  linear  time- varying  electric  current.  Figure  (5.5). 
This  comes  directly  from  the  nature  of  the  model  itself.  It  models  a  dumped 
system  and  the  angular  velocity  of  the  liquid  can  be  increased  only  providing 
an  increasing  torque  acting  on  the  liquid. 


5.3  Current  and  Voltage  drive 

The  FDHM  is  able  to  simulate  both  the  current  and  voltage  drive  Figures 

(5.6)-(6.7). 


X  10"^ 


Figure  5.6:  FDHM  -  Current  Drive 


X  io~' 


Figure  5.7:  FDHM  -  Voltage  Drive 
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The  device  appears  to  be  characterized  by  non-linear  effect  in  the  neighbour¬ 
hood  of  the  null  velocity.  It  can  be  seen  by  the  trend  of  the  electrostatic 
potential  for  the  current  drive  and  of  the  electric  current  in  the  case  of  voltage 
drive.  The  reason  can  be  found  in  the  modihcation  of  the  velocity  prohle  before 
than  the  fully  developed  MHD  flow  is  reached. 


5.4  Angular  Momentum  and  Power  Consump¬ 
tion 

The  FDHM  can  give  information  about  power  and  angular  momentum,  Figures 
(5.8)  -  (5.9). 

The  power  consumption  of  this  devices  is  typically  very  low  due  to  the  high 
electric  conductivity  of  the  working  fluid.  The  electrical  resistance  of  the  fluid 
at  rest  is  of  the  order  of  10“®  D;  while  the  liquid  increases  its  spin  velocity 
the  back  electromotive  force  increases  pulling  the  electric  current  to  flow  in 
narrows  layers  called  Hartmann  layers  parallel  to  the  z  axis. 


Time  —  [s] 


Figure  5.8:  FDHM  -  Power 
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X  10  ’ 


Figure  5.9:  FDHM  -  Angular  Momentum  F 

This  leads  to  increasing  values  of  the  electrical  resistance  with  increasing  power 
consumption. 

If  the  device  is  supplied  with  a  linearly  increasing  electric  current,  the  angular 
momentum  has  a  linear  trend  with  respect  of  the  time,  Figure  (5.9).  This  is  a 
direct  consequence  of  the  Torque- Current  characteristic  of  this  kind  of  device 
for  which  the  torque  provided  is  constant  if  the  electric  current  is  linear  over 
the  time. 


5.5  Simulations 

In  order  to  test  the  accuracy  of  the  FDHM  and  in  order  to  understand  the 
performances  of  the  device,  several  simulations  have  been  done. 

The  simulations  have  been  computed  varying  the  intensity  of  the  magnetic  field 
Hq,  the  electric  conductivity  of  the  liquid  a,  the  inner  radius  r*  and  the  height 
h  of  the  the  torus  containing  the  liquid.  Other  parameter  have  been  considered 
fixed,  such  as  the  outer  radius  rg,  the  viscosity  of  the  liquid  fi  and  density  p. 
The  following  table  summarizes  the  grid  of  values  of  the  simulations. 

The  results  are  provided  in  terms  of  physical  quantities  that  will  be  described 
in  the  following  subsections.  The  results  are  ordered  according  to  electrical 
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\kg/m^'\ 

1000 

1000 

1 

[Pa  s\ 

10-3 

10-3 

1 

[S/m] 

10® 

10 

6 

[T] 

0.8 

0.2 

4 

[mm] 

49 

45 

5 

[mm] 

50 

50 

1 

[mm] 

10 

30 

3 

Table  5.1:  Values  grid  for  the  FDHM  simulation 


conductivity  of  the  the  fluid  simulated.  In  order  to  highlight  the  capabilities  of 
the  FDHM,  the  numeric  error  at  the  end  of  the  each  simulation  of  the  dynamic 
and  electric  side  of  the  problem  are  presented. 

The  results  presented  in  the  next  sections  refers  only  to  liquid  with  electric 
conductivity  a  =  10  S/m,  in  5.6,  and  a  =  10®  S/m  in  5.7.  For  all  the  inter¬ 
mediate  values  of  electric  conductivity  the  results  are  given  in  the  Appendices: 
a  =  10^  S/m  in  Appendix  B,  a  =  10^  S/m  in  Appendix  C,  a  =  10"^  S/m  in 
Appendix  D,  a  =  10®  S/m  in  Appendix  E.  The  results  are  discussed  in  5.8. 
For  each  electrical  conductivity  both  the  current  drive  and  the  voltage  drive 
have  been  simulated.  For  the  current  drive  an  input  of  a  =  lA/ s  has  been 
given,  for  the  voltage  drive  b  =  lt)~^V/ s. 

The  fluid  has  been  considered  at  rest  for  t  =  Os  and  the  simulations  have  been 
stopped  when  the  torque  reached  the  plateau. 


5.5.1  Angular  Momentum  and  Torque 

The  angular  momentum  F  during  the  simulation  is  computed  numerically. 
Each  cell  is  considered  as  an  inhnitesimal  hollow  cylinder  spinning  at  the  center 
cell  angular  velocity.  The  angular  momentum  is  then  computed  as  follows: 


K 


^”  =  Y  ^  (^'+1/2  -  S-1/2)  (r^+1/2  -  rt  -1/2) 


(5.5.1) 


k=l  l=l 


In  order  to  increase  the  results  readability,  we  decided  not  to  plot  the  increas¬ 
ing  angular  momentum  with  the  time,  but  the  specihc  value  of  the  angular 
momentum  F  for  an  electric  current  of  /  =  lA.  It  has  been  computed  divid¬ 
ing  the  angular  momentum  at  the  end  of  the  simulation  by  the  total  electric 
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current  flowing  through  the  conductive  liquid  at  the  end  of  the  simulation. 


r  = 


T 

¥ 


(5.5.2) 


EOS 


The  torque  provided  is  computed  with  a  first  order  hnite  difference  (same  as 
Crank-Nicolson  method)  of  the  total  angular  moment: 


pn  _ 


where  dt  is  the  integration  time  step. 


rn  _  r 


n—1 


dt 


(5.5.3) 


5.5.2  Characteristic  Time 

The  characteristic  time  Tc  has  been  computed  in  the  post-processing  phase, 
hnding  the  time  when  the  provided  torque  is: 

f  it)  =  (l  -  e"^)  V^ax  Tc  =  t  (5.5.4) 

The  characteristic  time  is  computed  by  interpolation  in  the  case  the  previ¬ 
ous  relationship  is  not  satished  exactly  in  some  of  the  time  integration  points. 


5.5.3  Power 


Because  of  the  number  of  simulations  and  due  to  the  fact  that  the  power  is 
constantly  increasing  with  the  time  because  of  the  increasing  supplied  electric 
current,  the  power  plotted  in  the  results  is  a  specihc  value  of  the  power.  We 
can  say  it  is  an  equivalent  electric  resistance. 

The  specihc  power  P  is  computed  taking  into  account  the  total  electric  power 
at  the  end  of  the  simulation,  as  the  product  between  the  voltage  between  the 
armatures  and  the  total  electric  current  howing  through  the  conductive  liquid, 
divided  by  the  the  square  of  the  total  electric  current: 


P  = 


'A(j)  P' 

'Af 

[  \ 

EOS 

P 

EOS 


So,  P  is  an  electric  power  over  square  electric  current. 


(5.5.5) 
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5.5.4  FDHM  Errors 

Together  with  the  physical  quantities  just  described,  in  the  results  can  be  read 
the  numeric  error  of  the  solution  of  the  electric  and  dynamic  side  of  the  problem 
for  each  simulation.  The  meaning  of  the  errors  is  the  same  described  in  the 
equations  Eq.(4.2.7)  and  Eq.(4.2.8). 
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5.6  Results:  FDHM  for  a  =  10  S/m 

5.6.1  Current  Drive 


^  jq-4  h  =  10mm 


L  —  [mm] 


a  =  10  S/m  —  50  mm 


^  j^-4  h  =  20mm 


L  —  [mm] 


X 10'^  ^  “  30mm 


L  —  [mm\ 


- So  =  0.2T 

- Bo  =  0.4T 

-  -  -  So  =  0.6  T 
- So  =  0.8T 


Figure  5.10:  Current  Drive  -  Torque  for  cr  =  10  S/m  -  FDHM 


h  =  10mm 


L  —  [mm] 


a  =  10  S/m  Toat  =  50  mm 
h  =  20mm 


L  —  [mm] 


h  =  30mm 


L  —  [mm] 


- So  =  0.2T 

- Bo  =  0.4T 

- Bo  =  0.6  T 

- So  =  0.8T 


Figure  5.11:  Current  Drive  -  Charaeteristic  time  for  cr  =  10  S/m  -  FDHM 


(7  =  10  S/m  =  50  mm 


h  =  10mm 


L  —  [mm] 


h  =  20mm 


L  —  [mm] 


h  =  30mm 


L  —  [mm] 


- So  =  0.2T 


- So  =  0.4T 


So  =  0.6  T 


- So  =  0.8T 


Figure  5.12:  Current  Drive  -  Specific  Power  for  cr  =  10  S/m  -  FDHM 
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O'  =  10  S/m  =  50  mm 


L  —  [mm] 


^  j^-4  h  =  20mm 


L  —  [mm\ 


X  ^  “  30mm 


L  —  [mm] 


- So  =  0.2T 


- So  =  0.4T 


- So  =  0.6T 


- So  =  0.8T 


Figure  5.13:  Current  Drive  -  Specific  Angular  Momentum  for  a  =  10  S/m  - 
FDHM 


a  =  10  S/m  Tcyu^t  —  50  mm 


^jq-6  h  =  10mm 


L  —  [mm] 


^  jq-6  h  =  20mm 


L  —  [mm] 


X 10'^  ^  “  30mm 


L  —  [mm] 


Figure  5.14:  Current  Drive  -  Relative  Dynamic  Error  for  a  =  10 
FDHM 


- So  =  0.2T 


- So  =  0.4T 


-  -  So  =  0.6  T 


- So  =  0.8T 


S/m  - 


O'  =  10  S/m  rout  =  50  mm 
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L  —  [mm] 


^  jq-12  h  =  30mm 


L  —  [mm] 


- So  =  0.2T 


- So  =  0.4T 


So  =  0.6  T 


- So  =  0.8T 


Figure  5.15:  Current  Drive  -  Standard  deviation  of  the  electric  current  for 
a  =  10  S/m  -  FDHM 
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5.6.2  Voltage  Drive 


(7  =  10  S/m  =  50  mm 


X 10’^  ^  “  10mm 


L  —  [mm] 


^  jq-7  h  =  20mm 


L  —  [mm] 


^  jq-7  h  =  30mm 


L  —  [mm] 


Figure  5.16:  Voltage  Drive  -  Torque  for  a  =  10  S/m  -  FDHM 
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Figure  5.17:  Voltage  Drive  -  Characteristic  time  for  ct  =  10  S/m  -  FDHM 


h  =  10mm 
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Figure  5.18:  Voltage  Drive  -  Specific  Power  for  a  =  10  S/m  -  FDHM 
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O'  =  10  S/m  =  50  mm 


Figure  5.19:  Voltage  Drive  -  Specific  Angular  Momentum  for  cr  =  10  S/m  - 
FDHM 


a  =  10  S/m  —  50  mm 


Figure  5.20:  Voltage  Drive  -  Relative  Dynamic  Error  for  cr  =  10  S/m  - 
FDHM 


(7  =  10  S/m  Taut  =  50  mm 


Figure  5.21:  Voltage  Drive  -  Standard  deviation  of  the  electric  current  for 
a  =  10  S/m  -  FDHM 
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5.7  Results:  FDHM  for  cr  =  10®  S/m 


5.7.1  Current  Drive 


<7  =  10^’5/m  Tout  =  50mm 


^  jq-5  h  =  10mm 


L  —  [mm] 


- So  =  0.2T 


- So  =  0.4T 


-  So  =  0.6  T 


- So  =  0.8T 


Figure  5.22:  Current  Drive  -  Torque  for  a  =  10®  S/m  -  FDHM 
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Figure  5.23:  Current  Drive  -  Characteristic  time  for  a  =  10®  S/m  -  FDHM 
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-  -  So  =  0.6  T 


- So  =  0.8T 


Figure  5.24:  Current  Drive  -  Specific  Power  for  a  =  10®  S/m  -  FDHM 
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O'  =  r^yu^t  —  50mm 


X 10’^  ^  “  10mm 


L  —  [mm] 


^  jq-5  h  =  20mm 


L  —  [mm] 


X 10'^  ^  “  30mm 


L  —  [mm] 


- So  =  0.2T 


- So  =  0.4T 


So  =  0.6  T 


- So  =  0.8T 


Figure  5.25:  Current  Drive  -  Specific  Angular  Momentum  for  a  —  10®  S/m  - 
FDHM 


o  =  10^*5/m  Ffynt  —  50mm 


Figure  5.26:  Current  Drive  -  Relative  Dynamic  Error  for  a  =  10® 
FDHM 
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Figure  5.27:  Current  Drive  -  Standard  deviation  of  the  electric  current  for 
a  =  10®  S/m  -  FDHM 
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5.7.2  Voltage  Drive 


a  =  lO^S/m  raat  =  SOmm 


Figure  5.28:  Voltage  Drive  -  Torque  for  a  =  10®  S/m  -  FDHM 
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Figure  5.29:  Voltage  Drive  -  Characteristic  time  for  a  =  10®  S/m  -  FDHM 
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Figure  5.30:  Voltage  Drive  -  Specific  Power  for  a  =  10®  S/m  -  FDHM 
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O'  =  r^yu^t  —  50mm 


X 10’^  ^  “  10mm 
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Figure  5.31:  Voltage  Drive  -  Specific  Angular  Momentum  for  a  —  10®  S/m  - 
FDHM 
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Figure  5.32:  Voltage  Drive  -  Relative  Dynamic  Error  for  a  =  10®  S/m  - 
FDHM 
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Figure  5.33:  Voltage  Drive  -  Standard  deviation  of  the  electric  current  for 
a  =  10®  S/m  -  FDHM 
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5.8  Discussion 

In  this  section  we  will  discuss  the  plots  present  in  the  previous  sections  in  or¬ 
der  to  highlight  and  justify  the  behaviour  of  the  device  and  its  sensitivity  to 
the  geometrical  dimensions  of  the  torus,  the  physical  properties  of  the  working 
fluid  and  the  intensity  of  the  magnetic  field  in  both  the  case  of  current  and 
voltage  drive. 

The  torque  provided  by  the  device  is  highly  dependent  by  the  intensity  of  the 
magnetic  field  applied.  In  particular  for  increasing  intensity  of  the  magnetic 
field  the  torque  increases.  It  is  possible  to  see  anyway  two  different  trends 
according  to  the  electrical  conductivity  of  the  working  fluid. 

For  low  electrical  conductive  fluid,  Figure  (5.10),  the  flow  is  almost  biparabolic, 
and  then  quite  different  from  a  MHD  flow  at  high  Ha.  The  viscous  shear  is 
mainly  linked  with  the  overall  dimensions  of  the  cross  section  of  the  torus. 
So  the  torque  provided  increases  for  increasing  height  h  and  width  L  of  the 
cross  section  because  the  viscous  shear  decreases.  Being  the  viscous  shear  in¬ 
dependent,  in  this  case,  by  the  magnetic  field,  the  provided  torque  increases  for 
increasing  Bo  because  the  Lorentz  force  is  proportional  to  the  magnetic  field 
applied. 

For  high  electrical  conductive  fluid.  Figure  (5.22),  the  viscous  shear  is  inversely 
proportional  to  the  thickness  of  the  Hartmann  and  Shercliff  layers,  which  are 
proportional  to  I/Bq,  and  to  the  overall  dimensions  of  the  cross  section.  The 
torque  provided  increases  for  increasing  height  h  and  width  L  of  the  cross  sec¬ 
tion.  It  increases  for  increasing  intensity  of  the  magnetic  field  but  less  than 
the  case  of  low  electrically  conductive  fluid.  The  reason  resides  in  the  fact  that 
now,  the  viscous  shear  is  proportional  to  the  intensity  of  the  magnetic  field. 
The  torque  provided  in  case  of  voltage  drive  is  affected  by  the  total  electric 
current  flowing  through  the  liquid.  Being  the  electric  current  influenced  by  the 
electric  resistance  of  the  fluid  we  can  see  in  Figure  (5.16)  and  Figure  (5.28), 
that  the  torque  is  not  proportional  to  the  height  h  of  the  device.  Furthermore 
for  high  electrically  conductive  liquids.  Figure  (5.28),  the  torque  increases  for 
decreasing  intensity  of  the  magnetic  field;  it  is  because  the  total  electric  cur¬ 
rent  increases  when  the  Hartmann  layers  become  thicker  and  this  has  a  major 
effects  on  the  resulting  Lorentz  force  than  the  intensity  of  the  magnetic  field. 
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The  same  cannot  be  said  for  low  conductive  fluid,  Figure  (5.10)  and  Figure 
(5.16),  where  the  torque  has  the  same  trend,  just  scaled  due  to  the  different 
drive. 

The  characteristic  time  Tc  gives  informations  about  the  response  of  the  system 
to  a  voltage  or  current  input  linear  with  the  time,  according  to  the  specihc 
drive. 

For  low  conductive  fluid.  Figure  (5.11)  and  Figure  (5.17),  the  characteristic 
time  is  the  same  for  both  current  and  voltage  drive.  This  is  a  proof  of  the  fact 
that  the  flow  is  almost  equivalent  to  classical  Poiseuille  flow  and  so  not  influ¬ 
enced  by  the  magnetic  held.  In  this  this  case,  the  characteristic  time  is  mainly 
inhuenced  by  the  moment  of  inertia  of  the  spinning  liquid,  that  increases  for 
increasing  height  h  and  width  L  of  the  cross  section. 

For  high  conductive  liquids.  Figure  (5.23)  and  Figure  (5.29),  on  the  contrary 
we  have  some  differences.  The  system  appears  to  have  a  faster  response  in  the 
case  of  voltage  drive  than  in  the  case  of  current  drive.  The  reason  is  in  the  way 
the  device  manages  the  electric  current.  For  current  drive,  the  device  has  to 
distributed  all  the  electric  current  given  by  the  power  supply  to  the  fluid.  This 
leads  the  system  to  be  slower  in  the  answer.  On  the  contrary  in  the  voltage 
drive  the  system  chooses  by  itself  the  current  according  to  the  voltage  applied 
and  to  the  shape  of  the  velocity  prohle.  In  Figure  (5.6)  and  Figure  (5.7)  we 
can  see  that  in  the  case  of  voltage  drive  the  electric  current  increases  faster  in 
the  hrst  phase;  this  gives  more  Lorentz  force  to  the  liquid  that  tends  faster  to 
a  fully  developed  laminar  MHD  flow  than  the  case  of  current  drive.  For  high 
electrically  conductive  fluids  we  can  see  that  the  characteristic  time  decreases 
for  increasing  intensity  of  the  magnetic  held,  as  the  increasing  viscous  shear 
contributes  to  reach  faster  a  fully  developed  laminar  MHD  how. 

The  specihc  power  P  has  the  meaning  of  equivalent  electric  resistance  and  it 
is  computed  using  the  voltage  applied  across  the  device  and  the  total  electric 
current. 

For  low  conductive  huid.  Figure  (5.12)  and  Figure  (5.18),  there  is  no  diherence 
between  voltage  and  current  drive  and  the  specihc  power  P  can  be  considered 
as  the  classical  electrical  resistance  of  the  huid  seen  as  simple  electric  conduc¬ 
tor. 

For  high  conductive  huid.  Figure  (5.24)  and  Figure  (5.30),  there  are  some  little 
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differences  and  these  are  due  to  the  different  answer  of  the  electric  side  of  the 
system  to  the  voltage  or  current  drive.  In  Figure  (5.6)  and  Figure  (5.7)  we  can 
see  that  the  voltage  in  the  case  of  current  drive  and  the  current  in  the  case  of 
voltage  drive  have  two  different  knees  with  opposite  concavity.  This  influences 
the  total  electric  power  and  as  consequence  the  specihc  power,  chosen  as  ref¬ 
erence  quantity. 

The  specihc  angular  momentum  F  is  the  angular  momentum  divided  by  the 
electric  current  howing  through  the  huid.  It  has  the  meaning  of  how  much 
angular  momentum  we  can  get  from  the  device  if  the  total  current  is  Xj A. 

In  the  case  of  current  drive  in  both  the  case  of  high  and  low  conductive  huid, 
Figure  (5.25)  and  Figure  (5.13),  we  have  that  the  specihc  angular  momentum 
is  numerically  the  same  as  the  torque.  The  reason  is  the  particular  value  of 
the  current  slope  in  time,  I  =  1  A,  that  makes  the  specihc  angular  momentum 
and  torque  numerically  equivalent.  In  the  case  of  voltage  drive.  Figure  (5.19) 
and  Figure  (5.31),  the  behaviour  of  F  changes  and  it  is  due  to  the  not  constant 
electric  current  over  the  several  simulation  made.  We  will  clarify  better  the 
behaviour  of  the  angular  momentum  in  the  second  phase  of  the  project  by 
means  of  a  lumped  parameter  model. 

Finally  in  the  results  can  be  found  the  computational  errors  made  in  the  so¬ 
lution  of  the  electric.  Figures  (5.15)-(5.27)  and  Figures  (5.21)-(5.33),  and  the 
errors  of  the  dynamic  part  of  the  problem.  Figures  (5.14)-(5.26)  and  Figures 
(5.20)-(5.32).  The  relative  error  in  the  solution  of  the  dynamic  part  is  of  the 
order  of  10“®  while  the  absolute  error  associated  to  the  electric  part  is  always 
less  than  10“^°.  They  are  all  sufficiently  low  to  consider  the  simulations  trust¬ 
worthy  and  it  is  also  a  proof  of  the  accuracy  of  the  FDHM  over  a  wide  range 
of  input  values. 
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Methods  and  Procedures:  Lumped  Pa¬ 
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In  this  chapter,  a  Lumped  Parameter  Model  similar  to  the  model  used  for 
classical  reaction  wheels.  The  coefficients  of  the  model  are  computed  on  the 
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6.1  Lumped  Parameter  Model 

The  lumped  parameter  model  of  a  classical  reaction  wheel  with  rigid  flywheel 
is  well-known  and  it  take  into  account  the  electric  and  dynamic  phenomena 
happening  while  the  flywheel  spins.  It  has  the  following  shape: 

{hdJw  =  KtI 

(6.1.1) 

^(j)  KyiOyj  T  R^ql 


The  purpose  of  this  chapter  is  to  derive  a  lumped  parameter  model  for  the 
device  in  study.  The  model  for  the  MHD  reaction  wheel  is  supposed  to  be  very 
similar  to  the  classical  reaction  wheel  one.  Some  contributions  need  anyway 
to  be  added  in  order  to  take  into  account  the  viscous  phenomena  that  are  ne- 
glectable  in  the  classical  reaction  wheel. 

The  Lumped  Parameter  Model  (LPM)  for  MHD  reaction  wheel  has  the  follow¬ 
ing  shape: 

{I zCijjUJm  KjjigUJiu  -f-  Kjl 

(6.1.2) 


In  the  model  the  angular  velocity  of  the  fluid  has  the  meaning  of  mean 
relative  angular  velocity  between  the  fluid  and  the  case  of  the  reaction  wheel. 
The  model  is  similar  to  the  model  of  a  classical  reaction  wheel  with  an  extra¬ 
term  taking  into  account  the  viscous  torque,  proportional  to  the  relative  angu¬ 
lar  velocity  This  term  leads  to  a  slight  different  behaviour  of  the  Torque- 
Current  characteristic:  while,  for  the  classical  reaction  wheel,  the  torque  pro¬ 
vided  is  constant  when  the  electric  current  fed  to  the  device  is  constant,  here, 
the  torque  is  constant  if  the  electric  current  varies  linearly  over  the  time.  This 
come  directly  from  the  nature  of  the  device  itself.  It  is  modeled  as  a  dumped 
system  and  the  angular  velocity  of  the  liquid  can  be  increased  only  providing 
an  increasing  torque  acting  on  the  liquid. 

The  first  equation  of  (6.1.2)  models  the  dynamic  part  of  the  problem  stating 
that  the  provided  torque  is  the  resultant  between  the  moment  of  the  viscous 
shear  and  moment  of  the  Lorentz  force;  the  second  equation  of  (6.1.2)  models 
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the  electric  part  of  the  problem  stating  that  the  applied  voltage  is  equal  to 
sum  of  the  resistive  voltage  losses  and  the  induced  voltage  generated  by  the 
electrically  conductive  liquid  spinning  through  the  magnetic  held. 

In  the  following  two  sections  the  coefficients  of  the  dynamic  and  electric  parts 
of  the  model  will  be  derived  [11]. 


6.2  Lumped  Parameter  Dynamic  Model 

In  order  to  derive  the  dynamic  part  of  the  Lumped  Parameter  model  we  con¬ 
sider  the  Navier  Stokes  equation  through  an  integral  balance  of  moments  along 
the  direction. 


pri  X  ^^6dV  =  J  rr  x  [J^B  (r)]  OdV 

,  ,  .  .  d  ru> 

+  /  prr  X 

Js±r 


dr  \  r 


OdS  +  /  prr  x 

Js±z 


m 

dz 


eds 

(6.2.1) 


A  stationary  solution  needs  to  be  derived  and  used  to  obtain  quantitative 
evaluation  of  the  viscous  shear  at  the  boundaries  together  with  an  equivalent 
moment  of  inertia  of  the  spinning  liquid  and  the  back  electromotive  force  gen¬ 
erated  by  the  spinning  liquid  inside  the  superimposed  magnetic  held. 

To  compute  the  moment  of  the  viscous  shear,  we  consider  the  last  term  of 
right  side  of  Eq.  (6.2.1)  and  we  obtain  for  the  boundaries  of  the  rectangular 
cross  section: 


[  pr^^(^]dS  +  [  pr^dS  (6.2.2) 

Js±r  or  \  r  /  Js±z  oz 

where  S'  T  r  and  S  P  z  are  respectively  the  surfaces  of  the  torus  parallel  to  the 
and  r  axes.  For  the  terms  on  the  right  of  the  Eq.  (6.2.2)  a  detailed  expression 
can  be  given: 
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=27r/irf 


du^  (r,  z) 


dr 


dz 


r=ri 


27ifiri 


dujz  {r,  z) 


dr 


dz 


(6.2.3) 


dUg 

ur^—db 

dz 


=27ifa  /  r' 


.du^  {r,z) 


'n 


dz 


dr 


z=0 


+  2nfi  /  r 


^duz  ir,z) 


dz 


dr 


z=h 


(6.2.4) 


The  moment  generated  by  the  viscous  shear  can  then  be  written  as: 


■^vis  ^vis^w 


(6.2.5) 


where: 


2ti jarf  duz  (r,  z) 


K, 


u. 


w  Jo 


dr 


2TTfj,rl  duz  (r,  z) 


dz  T 


r=ri 


u. 


w  Jo 


dr 


2nfi  3  dojz  (r,  z) 


CU-ji 


dz 


dr  +  ^ 


2=0 


CU-ji 


dz 


dz 

(6.2.6) 

dr 

z=h 


The  term  Uz  {r,  z)  acquires  the  meaning  of  dimensionless  angular  velocity 
normalized  with  respect  of  the  mean  angular  velocity. 

The  torque  generated  by  the  Lorentz  Force  can  be  derived  from  the  hrst  term 
of  the  right  side  of  the  (6.2.1).  Introducing  the  shape  of  the  radial  magnetic 
held  we  have: 


Mior  =  r  X  9  {JzBr)  dV  =  Bf^rihl  =  Kjl 


(6.2.7) 


where: 
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Kj  =  B^rih  (6.2.8) 

The  term  on  the  left  side  of  the  (6.2.1)  represents  the  resulting  torque  acting  on 
the  conductive  fluid,  i.e.  the  time  derivative  of  the  angular  moment.  In  order 
to  have  this  term  function  of  the  only  time  derivative  of  the  average  relavite 
angular  velocity  of  the  fluid  we  start  considering  the  angular  moment  in  the 
reference  frame  of  the  spacecraft: 

r  =  /  pr‘^ujz  (r)  zdV  =  (6.2.9) 

Jv 


where: 


lo  (r,  z)  drdz 

^ ^ -  (6.2.10) 

UJuj  Iz 

and  Iz  is  the  moment  of  inertia  of  a  rigid  hollow  cylinder  with  the  same  dimen¬ 
sions  of  the  cavity  containing  the  liquid.  In  this  way  the  time-derivative  of  the 
angular  moment  can  be  considered  as: 

t  =  Uy,cJz  (6.2.11) 

due  to  invariance  with  respect  of  the  time  of  the  coefficient  c^. 

6.3  Lumped  Parameter  Electric  Model 

A  formula  for  the  lumped  parameter  electric  model  of  the  proposed  MHD 
reaction  wheel  can  be  obtained  by  integrating  over  the  volume  the  microscopic 
Ohm’s  Law  for  moving  conductor  and  introducing  the  relative  velocity  between 
the  liquid  and  the  magnetic  held  considered  stiffiy  linked  to  the  boundaries: 

[  EzdzdS  =  [  —dzdS+  [  ruz{r)  BrdzdS  (6.3.1) 

Jv  Jv  ^  Jv 
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Considering  the  electrostatic  potential  0  uniform  over  the  surface  and  applying 
the  conservation  of  the  electric  current  we  obtain: 


h  ruz  (r,  z)  BMzdS 
=  /—  +  - 


(6.3.2) 


where  the  coefficient  of  the  electric  current  is  the  equivalent  electric  resistance 
of  the  conductive  liquid: 

Req  =  ^  (6.3.3) 

The  second  term  on  the  right  side  of  the  (6.3.1)  is  the  induced  voltage  that  can 
be  written  as: 

Acpind  =  KyUJw  (6.3.4) 


where: 


Ky  =  2TiBQri 


ruz  (r,  z)  drdz 


(jJyjS 


(6.3.5) 


6.4  Analytical  solution  for  the  angular  and  lin¬ 
ear  velocity 

The  principal  problem  to  be  solved  to  derive  the  LPM  is  to  hud  an  analytical 
solution  for  the  angular  and  linear  velocity  of  the  conducting  liquid.  This  step 
is  necessary  to  compute  the  LPM  coefficients  which  appear  to  be  function  of 
this  distributions  and  because  the  LPM  is  referred  to  the  average  angular  of 
the  fluid. 

In  this  paragraph  two  different  solutions  will  be  developed  and  then  compared 
with  the  numerical  solution  provided  by  the  FDHM  model. 


6.4.1  Monodimensional  solution  under  the  hypothesis 
of  Low-Magnetic  Reynolds 

For  convenience  we  proceed  writing  the  whole  set  of  equation  in  dimensionless 
variables.  We  obtain: 


V'  ■  b  =  0 


(6.4.1) 
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v'-j  =  0 

(6.4.2) 

j  =  e  +  u  X  b 

(6.4.3) 

V'  X  e  =  0 

(6.4.4) 

du  1  ,  Ha'^ . 

(6.4.5) 

r 

s  =  — 

(6.4.6) 

Te 

V'  ■  u  =  0 

(6.4.7) 

in  which  the  following  substitutions  have  been  done: 

u 

(6.4.8) 

Bo 

(6.4.9) 

E 

Bo  Uo 

(6.4.10) 

3 

a  Bo  Uo 

(6.4.11) 

tUo 
^  L 

(6.4.12) 

V'  =  L  V 

(6.4.13) 

The  coefficients  in  Eq.  (6.4.5)  are  the  Hartmann  Number,  Ha  defined 

as  BoL^Ja/ p, 

and  the  Reynolds  number,  Re,  defined  as  pLUo/fi.  In  Eqs.(6.4.8)-(6.4.13)  Uq 
has  the  meaning  of  average  velocity  of  the  fluid,  Bq  is  the  magnetic  field  density 
at  the  outer  radius  and  re  is  the  outer  radius  of  the  liquid  flywheel. 


The  main  simplifying  assumption  adopted  is  to  consider  infinite  the  length 
of  the  annulus  along  the  axis  of  symmetry.  It  implies  the  neglection  of  the 
boundary  effects  next  to  the  two  electric  armatures  and  it  implies  moreover 
that  all  the  derivative  along  are  considered  neglectible.  The  superimposed 
magnetic  field  is  assumed  to  be  dependent  just  from  the  radial  coordinate. 
This  assumption  together  with  the  uniformity  of  the  imposed  electric  field  e 
makes  the  problem  axially  symmetric,  i.e.  all  the  derivative  with  respect  of  the 
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azimuthal  direction  can  be  neglected. 


Figure  6.1:  Section  of  an  infinite  annulus  with  applied  electric  field  E  and 
radial  magnetic  field  B  (r) . 

We  start  by  analyzing  the  continuity  equations  of  the  magnetic  flux  density 
b,  current  density  j  and  velocity  u. 

It  is  possible  to  simplify  the  continuity  equations  for  the  dimensionless  mag¬ 
netic  flux  density  b,  Eq.(6.4.1),  the  dimensionless  electric  current  density  j, 
Eq.(6.4.2),  and  the  dimensionless  velocity  of  the  fluid  u,  Eq.(6.4.7).  Introduc¬ 
ing  the  expression  of  the  divergence  in  cylindrical  coordinates  and  simplifying 
we  obtain: 


Vb 


s  ds 


{sbr)  -|- 


1  dbg 

s  do 


A 

ds 


(sbr)  =  0 


(6.4.14) 
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1  (9  ,  ,  ,  1  djg  djz  d  /  ■  \ 


s  ds ' 


VU  =  -  —  (sUr) 


1  dug  dUz  d  , 


with  the  following  boundary  condition  at  the  outer  radius: 


(6.4.15) 

(6.4.16) 


br{so,T)  =  1 


So  =  1  jr{so,T)  =  0 


(6.4.17) 


Ur{so,T)  =  0 


Solving  Eqs.(6.4.14)-(6.4.16)  using  Eqs. (6.4.17)  we  obtain: 


bris,T)  =  - 
s 


jr{s,T)  =  0 


Ur{s,  r)  =  0 


(6.4.18) 

(6.4.19) 

(6.4.20) 


The  axial  component  of  the  dimensionless  velocity,  Uz,  can  be  studied  taking 
into  account  the  Eq. (6.4.5)  along  2;: 


duz 

_  A2 

ri  d 

1 

(  duz\ 

dr 

Re 

s  ds 

+ 


Ha^ 

Re 


jebr 


(6.4.21) 


Using  the  hypothesis  of  axial-symmetry  it  is  possible  to  say  that  the  component 
along  6  of  the  current  density  jg  is  null  both  because  of  the  radial  distribution 
of  the  magnetic  held  and  because  of  the  uniformity  of  the  electric  held  applied. 
This  result  leads  the  Eq.  (6.4.21)  to  be  a  homogeneous  PDE  with  homogeneous 
boundary  condition  due  to  the  no-slip  condition  applied  at  the  boundaries  to 
the  velocity  held.  Furthermore  and  without  lack  of  generality  we  can  impose 
that  the  huid  is  at  rest  at  the  initial  time  r  =  0.  This  problem  gives  a  null 
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solution  over  the  time  for  the  Uz'- 


M^(s,r)  =  0  (6.4.22) 

Writing  the  Eq. (6.4.3)  along  the  three  directions  of  cylindrical  reference  frame 
and  being  the  dimensionless  velocity  purely  azimuthal  and  the  magnetic  flux 
density  purely  radial  we  have: 

Cr  =  jr  —  Uebz  +  UzhO  =  0 

=  je  —  UzK  +  Urbz  =  0  (6.4.23) 

Cz  =  iz  -  Urbe  +  uebr  =  jz  +  uebr 

The  Eq.  (6.4.4)  describing  the  irrotationality  of  the  applied  electric  held  is  solved 
taking  into  account  the  hypothesis  that  the  component  of  the  electric  held 
along  the  direction  is  uniform.  In  other  words  if  we  consider  each  armature 
equipotential  we  can  see  that  the  electric  held  along  the  2:  direction  does  not 
vary  along  the  radial  and  azimuthal  directions. 

Finally  we  obtain  the  complete  set  of  equation  for  magneto-huid  problems 
under  the  hypothesis  of  Low-Magnetic  Reynolds  reduced  to  the  the  following 
system  of  three  equations: 


br  = 


1 

s 


(6.4.24) 


jz  G-z  Ugbj- 


(6.4.25) 


due 


Re  ds 


s  ds 


{sue) 


Ha‘^  ,  , 


(6.4.26) 


in  which  Eq.(6.4.24)  is  already  solved  and  just  Eqs.(6.4.25)-(6.4.26)  are  coupled. 
Starting  from  the  hypothesis  of  Low-Magnetic  Reynolds  the  same  equations  as 
in  [25]  has  been  found  and  we  can  use  the  same  solution  proposed  in  the 
paper. 
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Writing  the  Eq. (6.4.26)  with  the  same  symbology  as  in  Reference  [25]  and 
introducing  the  Eq. (6.4.25)  we  have: 


due 


Cl 


d'^ue  1  due 
ds^  s  ds 


u?ue 


+  C2- 


(6.4.27) 


where: 


Cl  =  X^/Re  C2  =  Ha^/Re 

A  =  L/R  n^  =  l  +  Ha^/X^ 


(6.4.28) 


In  [25]  two  solutions  are  given:  one  steady-state  solution  and  one  time-dependent 
solution.  Looking  at  the  purpose  to  obtain  a  lumped  parameter  model  we  show 
more  interest  in  the  stationary  solution  in  which  the  influence  of  the  electric 
current  and  the  shape  of  velocity  prohle  can  be  easily  uncoupled  instead  of 
considering  the  transient  solution  that  gives  a  time  dependent  solution  but  by 
which  is  not  possible  to  derive  a  lumped  parameter  model. 

The  steady-state  solutions  for  the  velocity  and  the  angular  velocity  are: 


Ue  (r)  =  T 

L>0 


Air-A2r^-A3- 


CU.  r  = 


Bn 


-  A 


r.n+1 


(6.4.29) 

(6.4.30) 


where: 


Ai  — 


1 

R 


A2 


1  -  ^  (1  -  R 

(1  -  /32^)  R”  ^  ~  1  - 


(6.4.31) 


and 


/3  =  1-X 


(6.4.32) 
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Figure  6.2:  Analytical  solutions  for  dimensionless  linear  and  angular  veloci¬ 
ties  at  different  Hartmann  Number  Ha. 


r/re 


Figure  6.3:  Magnetic  Field  Density  B  normalized  with  respect  of  the  Magnetic 
Field  Density  Bq  at  the  outer  radius. 

The  shape  of  the  dimensional  magnetic  field  B  can  be  obtained  from  the 
Eqs.(6.4.9)-(6.4.24): 

B  (r)  =  (6.4.33) 

r 

where  Bq  is  the  value  of  the  magnetic  held  at  r  =  r^.  The  solution  is  expressed 
in  dimensional  quantities.  The  coefficients  n  and  (5  in  Eq.  (6.4.29)  are  function 
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of  the  physical  property  of  the  conductive  fluid  and  of  the  dimension  along  the 
direction  r  of  the  cavity  containing  the  fluid.  The  solution  can  be  considered 
expressed  in  the  reference  frame  of  the  spacecraft  with  the  meaning  of  relative 
velocity  between  the  liquid  and  the  case.  With  respect  of  [25]  this  solution  has 
been  derived  explicitly  under  the  hypothesis  of  Low-i?em  and  this  can  help  to 
understand  the  range  of  validity  of  the  solution. 

6.4.2  Bi-dimensional  solution  extended  to  the  cross-section 

The  first  attempt  is  to  hud  an  asymptotic  bi-dimensional  analytical  solution 
for  the  complete  set  of  the  MHD  equations  in  cylindrical  coordinates  with 
radial  magnetic  held  and  axial  electric  current.  Hunt  and  Baylis  [12]  already 
succeeded  in  hnding  this  analytical  solution  but  in  the  case  of  axial  magnetic 
held  and  radial  electric  current.  The  same  procedures  can  be  followed  to  try 
to  hud  a  solution  for  the  case  in  study. 

The  complete  set  of  equations  needs  to  be  used: 


^  T.  dB 

(6.4.34) 

V  ■  B  =  0 

(6.4.35) 

V  X  H  =  47rJ 

(6.4.36) 

V- J  =  0 

(6.4.37) 

=  a(E  + V  X  B) 

(6.4.38) 

-I-  z/pAV  -I-  J  X  B  =  0 

(6.4.39) 

V- V  =  0 

(6.4.40) 

Expanding  the  Eq. (6.4.38)  and  considering  that: 


E  =  ErV  +  EzZ 


d4>  ^  d4>  ^ 


(6.4.41) 
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we  have: 


_  <90 

■Jr  - 

or 


(6.4.42) 


Jz  =  -a 


+  Ue  9  X 


(6.4.43) 


Applying  to  Eq.  (6.4.42)  and  Eq. (6.4.43)  respectively  the  derivative  with  respect 
of  and  r  we  obtain: 


dJr  d^^cj) 

dz  drdz 


(6.4.44) 


dr 


q2^  d(Ue9x  b) 
drdz  dr 


(6.4.45) 


Taking  now  into  account  the  Eq.  (6.4.36)  we  have: 

1  1  dHe 

An  dz 


_  ^ld{rHe) 
^2  -  . 


47r  r  dr 


Proceeding  as  before  we  obtain: 


d.L 


dz 


1  d^He 
An  dz'^ 


dr  47r 


IdHe^d^He  Hg 


r  dr 


dr '2 


It  is  possible  then  to  compare  the  Eq. (6.4.44)  with  Eq. (6.4.48): 

1  d‘^Hg 


dJr  d'^(f) 

dz  drdz 


47r  dz"^ 


and  the  Eq. (6.4.45)  with  Eq. (6.4.49): 


(6.4.46) 

(6.4.47) 


(6.4.48) 

(6.4.49) 


(6.4.50) 
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dJ^ 

dr 


=  a 


q2.  d(Ue9xB 
+  - 


drdz 


dr 


1 

47r 


IdHe  d'^He 

‘T‘ 


He 

rr2 


(6.4.51) 


Subtracting  the  Eq.  (6.4.50)  to  the  Eq.  (6.4.51)  we  have: 


IdH, 


r  dr 


^e  d^He  d‘^He 


e 


dz"^ 


TT  diUeOxB^ 


r 


dr 


■■  0  (6.4.52) 


Considering  a  perfect  radial  applied  magnetic  held  with  equation: 


B  =  fin- 


(6.4.53) 


The  last  term  of  Eq.  (6.4.52)  can  be  written  as: 


diUedxB 


dr 


T"  '^2 


(6.4.54) 


We  hnally  obtain: 


d^He  ,  d^He  ,  IdHe  He  ,  ,  ^  fldUe  Ue  ,  ^ 

+  ~  ~  +  47r(T5ori  (  -  —  )  =  0  (6.4.55) 


dr'i  dz"^  r  dr 


■^2 


The  Navier  stokes  equation  Eq.  (6.4.39)  along  the  6  direction  is: 

1  d  (  dUe\  ,  d^Ue  Ue\  ^ 

h  I  I  r—  I  +  - ^  I  +  J^B  =  0 


r  dr  \  dr  dz"^  r^ 


(6.4.56) 


With  the  aid  of  the  Eqs.(6.4.47),(6.4.53)  the  last  term  become: 


J,B  =  Bor,  (  +  \He 

ly  (P/'y 


(6.4.57) 
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and  then: 


dr'i  dz'^  r  dr  n  dr  j 


(6.4.58) 


What  can  be  seen  is  that  Eqs. (6.4.56)  and  (6.4.58)  have  similar  shape.  The 
hrst  fonr  terms  of  each  eqnation  have  the  same  operators  respectively  applied 
to  Ue  and  Hq.  The  last  terms  are  different:  they  have  different  coefficient  and 
different  signes  inside  the  brackets.  The  traditional  technique  is  to  add  the  two 
equations  after  having  applied  a  non-dimensional  substitution  and  then  solve 
the  equation  in  the  variable  resulting  in  the  sum  between  Ue  and  Hq.  In  this 
case  it  is  not  possible  because  of  the  shape  of  the  last  term.  We  then  try  to 
hnd  an  approximated  solution  that  will  be  compared  with  the  results  obtained 
from  the  FDHM. 

6.4.3  Cartesian  bi-dimensional  solution 

The  axially-symmetric  MHD  problem  extended  to  a  torus  with  square  section 
has  different  solution  compared  with  the  solution  that  can  be  obtained  from 
the  same  problem  in  the  cartesian  form  for  a  straight  pipe. 


Analytical  Solution  2D  —  Straight  pipe  FDHM  —  Angular  Velocity 


Figure  6.4:  Comparison  between  MHD  velocity  field  in  a  straight  pipe  under 
trasverse  magnetic  field  and  MHD  angular  velocity  for  an  axially- 
symmetric  torus  under  a  radial  magnetic  field. 
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There  are  anyway  some  qualitative  analogies.  As  can  be  seen  in  the  results  of 
the  FDHM  and  in  Figure  (6.4),  the  angular  velocity  of  the  conductive  liquid 
shows  a  rigid  body-like  behaviour  of  the  liquid  with  constant  angular  velocity 
in  core  zone  of  the  flux.  Qualitatively  this  distribution  can  be  compared  with 
the  linear  velocity  prohle  solution  of  the  MHD  problem  in  cartesian  coordinates 
for  a  straight  pipe  with  rectangular  cross  section. 

The  basic  equations  to  solve  the  problem  are,  in  dimensionless  form  [13]: 

3u 

V^fe  +  i/a— =  0  (6.4.59) 

dy 


-I-  Ha 


db 

dy 


0 


(6.4.60) 


The  two  equations  can  be  solved  together  with  the  Elsasser  variables: 


A  =  u  +  b,A'  =  u  —  b 


(6.4.61) 


Adding  and  subtracting  Eq. (6.4.60)  from  the  Eq. (6.4.59)  we  obtain: 

d^A  d^A 


dy‘^ 

d^A 

dy"^ 


+ 


dz^ 

d^A' 

dz‘^ 


TT  1 

-|-  Ha— —  —  — 1 
dy 

(6.4.62) 

Ha  =  1 

dy 

(6.4.63) 

The  boundary  conditions  to  be  applied  are  no  slip  condition  at  all  walls,  n  =  0, 
and  considering  the  boundaries  electrically  insulating,  the  induced  magnetic 
held  vanishes  at  the  boundaries,  b  =  0.  If  the  boundaries  is  called  F,  it  means 
that: 


A  =  a!  =  Q  on  F 

Moreover  it  can  be  said  by  symmetry  that: 

A'  (y)=A{-y)  =  Q 


(6.4.64) 


(6.4.65) 
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so  that  it  is  sufficient  to  determine  A  for  the  solution  of  the  problem.  The 
solution  is  obtained  by  expanding  the  unknown  A  into  a  Eourier  series  as: 

OO 

A  (y,  z)  =  A'  {-y,  z)  =  ^  a*  (y)  cos  {\iz) 

i=l,3,5 

(6.4.66) 

where  we  satisfy  the  boundary  conditions  aX  z  =  id  with: 

in 

/A  7  —  I 

2d 

(6.4.67) 

Using  the  Eq. (6.4.61),  for  the  dimensionless  linear  velocity  u 
sionless  induced  magnetic  held  b,  we  have: 

and  the  dimen- 

OO 

u  (y,  z)  =  Ui  (y)  cos  (Xiz) 

2=1, 3, 5 

(6.4.68) 

OO 

b  (d,  z)  =  'Yh  (d)  cos  {Xiz) 

(6.4.69) 

i=l,3,5 


Inserting  Eq.  (6.4.66)  into  Eq.  (6.4.62)  and  using  orthogonality  of  trigonometric 
functions  yields: 


d‘^ai  2  ,  u-  ^sin(Aj2:) 


-ki 


(6.4.70) 


whose  solution  is: 


cii  {y) 


A? 


+  Ci^  cosh  exp  -pi^y  +  exp  -pi^y 


(6.4.71) 


where: 


p,,,,  =  l(HaT  +  (6.4.72) 

In  order  to  obtain  the  coefficient  of  the  Eqs.(6.4.67),(6.4.68)  we  can  use  the 
symmetry  of  the  problem  and  Eq. (6.4.61): 
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Ui  iy) 


ttj  jy)  +  aj  j-y)  ^ 

2 

k- 

^  +  Cii  cosh  {pi^y)  -  Ci^  cosh  {-pi^y) 


(6.4.73) 


,  /  N  ^iiy)-  tti  {-y) 

bi  [y)  = - ^ - = 

=  -Ci^  cosh  {pi^y)  -  Ci^  cosh  {-pi^y) 


(6.4.74) 


which  correspond  to  the  Fourier  modes  for  velocity  and  magnetic  held.  Both 
vanish  at  the  Hartmann  walls,  Ui  (1)  =  bi  (1)  =  0,  so  that  the  coefficients  are 
determined.  We  express  the  Fourier  modes  for  velocity  and  induced  held  in 
Eqs.(6.4.67),(6.4.68)  as: 


Ui  {y) 


^  L  _  fiiyY 
AH  /.(I). 


h  iy) 


h  \ 9i  jy) ' 
A?  U(i). 


(6.4.75) 

(6.4.76) 


where: 


fi  {y)  =  ai^  cosh  {pi^y)  -  at,  cosh  (pi^y)  (6.4.77) 

9i  (y)  =  sinh  (pi^y)  -  ttii  sinh  (pi^y)  (6.4.78) 


and: 


Oil  2  =  sinh  (6.4.79) 

The  solution  u{y,  z)  found  will  now  be  taken  as  it  was  the  angular  velocity 
prohle  Uz  instead  of  the  linear  velocity  prohle.  To  do  this,  the  dimensionless  y 
axis  will  be  renamed  r.  We  then  have: 
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uJz  (r,  z) 


E 

i=l,3,5 


hi 

a! 


h 


2r 


re  +  Vi 


m 


cos 


(6.4.80) 


6.5  Conclusions 

In  this  chapter  a  Lumped  Parameter  Model,  LPM,  has  been  proposed.  The 
numerical  values  of  the  coefficients  of  the  LPM  can  be  computed  analytically 
on  the  basis  of  two  different  analytical  solutions:  a  mono-dimensional  analyt¬ 
ical  solution  and  a  bi-dimensional  analytical  solution.  The  expression  of  each 
coefficient  is  given  in  closed  form  function  of  the  geometry  of  the  torus,  of  the 
physical  properties  of  the  fluid  and  function  of  the  intensity  of  the  magnetic 
held.  It  is  not  known  a  priori  the  precision  of  the  coefficients  computed  with 
both  the  solutions.  In  the  next  chapter  the  LPM  is  then  compared  with  the 
results  obtained  with  the  FDHM. 
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In  this  chapter,  a  comparison  between  the  FDHM  and  the  LPM  is  shown  in 
order  to  understand  the  precision  of  the  estimation  of  the  LPM  coefficients 
based  on  both  the  analytical  solutions. 
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7.1  Simulations 

In  order  to  test  the  accuracy  of  the  lumped  parameter  model  several  simulation 
have  been  done  with  the  FDHM.  The  aim  is  to  compare  the  value  of  the 
coefficients  obtained  through  the  LPM  with  the  ones  obtained  numerically  with 
the  FDHM.  The  analysis  will  be  done  with  all  the  coefficients  and  extended  to 
other  two  quantities:  the  torque  provided  by  the  device  and  the  characteristic 
time  of  the  torque.  The  torque  and  the  characteristic  time  will  be  studied  for 
both  the  voltage  and  current  drive:  the  device  will  be  supplied  respectively 
with  a  voltage  or  current  trend  linear  over  the  time. 

The  LPM  will  be  tested  with  the  monodimensional  analytical  solution  and  the 
bidimensional  cartesian  solution  obtained  in  the  previous  chapter. 

The  simulations  have  been  computed  varying  the  intensity  of  the  magnetic  held 
Bq,  the  electric  conductivity  of  the  liquid  a,  the  inner  radius  r*  and  the  height 
h  of  the  the  torus  containing  the  liquid.  Other  parameter  have  been  considered 
hxed,  such  as  the  outer  radius  r^,  the  viscosity  of  the  liquid  fx  and  its  density 
p.  The  relative  errors  between  the  LPM  and  the  FDHM  are  shown  in  order 
to  highlight  the  accuracy  of  the  solution  computed  by  the  LPM.  The  following 
table  summarizes  the  grid  of  values  of  the  simulations. 


[kg/m^\ 

1000 

1000 

1 

[Pa  s] 

10-3 

10-3 

1 

[S/m] 

10® 

10 

6 

[T] 

0.8 

0.2 

4 

[mm] 

49 

45 

5 

[mm] 

50 

50 

1 

[mm] 

10 

30 

3 

Table  7.1:  Values  grid  for  the  LPM  simulation 


The  results  presented  in  the  next  sections  refers  only  to  liquid  with  electric 
conductivity  a  =  10  S/m  and  a  =  10®  S/m.  For  all  the  intermediate  values 
of  electric  conductivity  the  results  are  given  in  the  Appendices:  a  =  10^  S/m 
in  Appendix  F,  ci  =  10®  S/m  in  Appendix  G,  cr  =  10^  S/m  in  Appendix  H, 
a  =  10®  S/m  in  Appendix  1. 
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7.2  Coefficient  of  the  viscous  shear  moment 


The  coefficient  of  the  viscous  shear  moment  for  the  LPM  based  on  the  Monodi¬ 
mensional  analytical  solution  can  be  expressed  as: 


K, 


jur^h  dooz  (r) 


OJ-ji 


dr 


+ 


2n fir^h  duz  (r) 


(jJti 


dr 


(7.2.1) 


The  derivatives  present  in  Eq. (7.2.1)  can  be  computed  analytically. 

In  the  bidimensional  LPM  all  the  boundaries  need  to  be  taken  into  account. 
The  coefficient  of  viscous  shear  moment  is  computes  as: 


K, 


2n jarf  duz  (r,  z) 


u. 


w  Jo 


dr 


dz  -\- 


2n nrl  duz  (r,  z) 


r=ri 


U. 


w  Jo 


dr 


dz 


2np,  3  dojz  (r,  z) 


(jJii 


dz 


dr  + 


271  fi  r‘'  3  duz  (r,  z) 


2=0 


(jJii 


dz 


(7.2.2) 


dr 


z=h 


The  Kyis  obtained  with  the  Eq.  (7.2.1)  and  Eq.  (7.5.3)  need  to  be  compared  with 
the  numerical  Kyis  obtained  with  the  FDHM.  The  numerical  Kyis  is  obtained 
evaluating  the  velocity  profile  and  its  derivative  with  respect  of  r  and  at  the 
boundaries  and  then  multiplying  each  values  for  the  elementary  surface  it  is 
referred  to: 


K 

^'-VtSpoHM 


2/i7r 


UJy 

jJTT 


{^1+1/2  -  zi-1/2) 


1=1 

K 


duz  (r*,  zi)  duz  (re,  zi) 
n - - hre - 


dr 


dr 


i^k+1/2  -  rl-1/2) 


k=l 


duz  (rfc,  0)  ^  duz  (rfc,  h) 


dr 


dr 


(7.2.3) 


where  N  e  M  are  the  number  of  intervals  of  the  FDHM  computational  grid 
along  2:  and  r. 
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7.3  Coefficient  of  the  Lorentz  Force  moment 

For  both  the  mono-dimensional  and  bi-dimensional  LPM,  the  coefficient  of  the 
Lorentz  Force  moment  Kj  can  be  expressed  as: 

Ki  =  BoTih  (7.3.1) 


The  same  coefficient  can  be  computed  for  the  FDHM  data.  The  numerical  Kj 
is  calculated  as: 


K 


IpDHM 


Ef=l  {^1+1/2  -  Zl-1/2)  (rEl/2  -  ^Li/2 


(7.3.2) 


For  the  current  drive  /  is  the  electric  current  fed  to  the  system.  In  the  voltage 
drive  the  electric  current  is  not  know  a  priori,  so  it  is  taken  as  the  average 
current  flowing  through  the  device.  This  current  is  computed  as: 


I  = 


K 


1=1  k=l 


(7.3.3) 


So  j  comes  directly  from  the  nodal  analysis  in  the  FDHM  and  it  is  the  av¬ 
erage  of  the  electric  currents  in  the  branches. 


7.4  Coefficient  of  the  moment  of  inertia 


The  coefficient  represents  the  difference  between  the  moment  of  inertia  of 
a  rigid  flywheel  with  the  same  dimensions  of  the  torus  containing  the  liquid 
and  the  liquid  spinning  in  the  torus.  In  the  LPM,  with  the  mono-dimensional 
solution,  we  have: 


c 


UJ 


2'Kph  r^Uz  (r)  dr 
Unjlz 


(7.4.1) 


In  the  bidimensional  LPM,  the  coefficient  is  computed  as: 


^UJ 


to -in  I z 


(7.4.2) 
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In  the  FDHM  the  c^j  is  computed  as: 


T^pl^k=il^i=i^ 


Zkd 


Cu,  = 


(^fc+1/2  -  4+1/2)  4l+l/2  -  Zl-1/2) 


2CU'3/i-/'r 


(7.4.3) 


7.5  Coefficient  of  Induced  Voltage 


The  coefficient  of  the  induced  voltage  calculated  in  the  LPM  with  the  monodi¬ 
mensional  solution  is: 


Ky  = 


27TBoh  (r)  dr 

oj^S 


(7.5.1) 


where  S  is  the  surface  perpendicular  to  In  the  bidimensional  LPM,  the 
coefficient  Ky  is  computed  as: 


Ky  = 


2ttBo  jl"  ruz  (r,  z)  drdz 
OJm  S 


The  coefficient  Ky  for  the  FDHM  can  be  computed  as: 

Bon  Ezii  Ef=i  {zi+1/2  -  ^z-i/2) 


K 


Vfdhm 


(7.5.2) 


(7.5.3) 


7.6  Current  Drive 

In  the  case  of  electric  current  drive  we  will  solve  the  dynamic  part  of  the  LPM 
for  a  generic  law  for  the  electric  current: 

u  +  —u  =  —I  (7.6.1) 

Z  ^  Z 

with  general  initial  condition: 

u  (to)  =  oJw  (7.6.2) 
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The  solution  for  u  is: 


_  4- 

u{t)  =  e  P 


K, 


+  /  —  /  (t)  e  dt 


'0 


(7.6.3) 


If  we  now  apply  a  linear  law  to  the  electric  current: 


/  (t)  =  at 


(7.6.4) 


the  solution  (7.6.3)  becomes: 


,  ,  Kj  iLKja 
oa  (t)  =  -rr^at  -  ^ 


K 


^^vis  L 


1  -  e" 


(7.6.5) 


and  the  torque  provided  is: 


T  =  I' u  =  /: 


l-e 


(7.6.6) 


The  (7.6)  shows  that  the  reaction  of  the  system  has  a  characteristic  time: 


ri  = 


K, 


(7.6.7) 


and  the  torque  provided  becomes  constant  for  a  linear  electric  current: 

Ki 


Too  =  lim  T  (t)  =  I'z  a 


t^OQ 


K, 


(7.6.8) 


The  angular  momentum  generated  by  the  device  is  computed  as: 


T  (t)  =  {t) 


(7.6.9) 


where  ca  (t)  is  computed  by  the  Eq. (7.6.5)  for  a  given  time  t.  The  power 
consumption  of  the  device  for  a  linear  current  law  is: 


P  (t)  =  (KyU)  (t)  +  ReqCd)  dt 


(7.6.10) 
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The  quantities  T  and  P  are  intended  as  specific  power  and  specific  angular 
momentum  for  an  electric  current  of  lA.  They  are  computed  as: 


T 


IzCui^  (tEOs) 
0-tEOS 


(7.6.11) 


P 


KyUJ  {IeOs)  +  PeqOtEOS 
atEOS 


(7.6.12) 


The  specific  power  and  the  specific  angular  momentum  in  the  FDHM  are 
computed  as: 

i'FDHM  =  — (7.6.13) 

J^EOS 

Pfdhm  =  (7.6.14) 

^EOS 

where  Ieos,  T eos  and  Peos  are  the  total  electric  current,  the  power  consump¬ 
tion  and  the  angular  momentum  computed  numerically  with  the  FDHM  data 
and  the  relative  to  the  end  of  the  simulation  [EOS).  The  simulations  have 
been  conducted  with  a  current  slope  a  =  lA/s. 


7.7  Voltage  Drive 


In  the  case  the  device  is  voltage  driven  we  need  to  consider  both  the  equations 
of  the  LPM.  Putting  them  together  through  the  electric  current  we  obtain: 


6j 


1 


KiKy 

Req 


Kj 


u  = 


P  P 


with  general  initial  condition: 


(7.7.1) 


U  (to)  =  UJw 


(7.7.2) 


The  solution  for  oo  is: 


u(t) 


-J-(k  ■  I 

=  p  C  Req  P 


+ 


Kj 


TD  T! 

JR'CqJ-z 


V(t) 


dt 


(7.7.3) 
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Considering  now  a  linear  law  for  the  voltage: 

V  (t)  =  bt 


(7.7.4) 


the  solution  (7.7.3)  becomes: 

Ki 


uj{t)  = 


RC I RCy  ~\~  -b^vis^eq 
I'.KiReqb 


bt  + 


{KjKv  +  KyisReqY 


1-e 


KjKy  +  K^^gReq  j 


(7.7.5) 


and  the  torque  provided  is: 


T  =  I'm  =  r 


Kj 


^  KjKv  +  KyisReq 


1-e 


I'.  Rea 


(7.7.6) 


The  (7.7.6)  shows  that  the  reaction  of  the  system  has  a  characteristic  time: 

(7.7.7) 


rv  = 


r/  p 


KjKv  +  KyisRe 


L/jvy  I 

and  the  torque  provided  becomes  constant  for  a  linear  voltage: 


Too  =  lim  T  it) 

t^OO 


Kb 


Ki 

KjKy  T  RyigR^q 


(7.7.8) 


The  torque  generated  in  the  FDHM  is  read  as  the  hnal  value  obtained  during 
the  simulation. 

The  angular  momentum  generated  by  the  device  is  computed  as: 


T  it)  =  IzCujU  (t) 


(7.7.9) 


where  u  (t)  is  computed  by  the  Eq.(7.7.5)  for  a  given  time  t.  The  power 
consumption  of  the  device  for  a  linear  current  law  is: 


p  ^  {bt-Kyu  {t))bt 
Req 


(7.7.10) 
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The  quantities  T  and  P  are  intended  as  specific  power  and  specific  angular 
momentum  for  an  electric  current  of  lA.  They  are  computed  as: 


T 


P 


(tEOs) 


eq 


btEOS  ~  KyOJ  (tEOs)) 


(7.7.11) 


P  = _ (7.7.12) 

ReqbtEOS  ~  ReqKv^  (tEOs) 

The  specific  power  and  the  specific  angular  momentum  in  the  FDHM  are 
computed  as: 

^FDHM  =  ~r~~  (7.7.13) 

J^EOS 

Pfdhm  =  ~j2~  (7.7.14) 

^EOS 

where  Ieosi  T eos  and  Pegs  are  the  total  electric  current,  the  power  consump¬ 
tion  and  the  angular  momentum  computed  numerically  with  the  FDHM  data 
and  the  relative  to  the  end  of  the  simulation  [EOS).  The  simulations  have 
been  conducted  with  a  voltage  slope  b  =  lQ~^V/s. 
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7.8  Results:  Mono  dimensional  LPM  -  a  =  lOS'/m 


cc 
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o  FDHM 
<  LPM -ID 


Figure  7.1:  Kyis 
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Figure  7.2:  Monodimensional  LPM 
Kj  for  (7  =  10  S/m  -  rout  =  50mm 


- Bo  =  0.2  T 

- Bo  =  0.4  T 
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Figure  7.3:  Monodimensional  LPM 
Cuj  for  a  =  10  S/m  -  rout  =  50mm 
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- Bo  =  0.4  T 

—Bo  =  0.6  T 

- Bo  =  0.8  T 

o  FDHM 
<  LPM -ID 


Figure  7.4:  Monodimensional  LPM 
Ky  for  a  =  10  S/m  -  rout  =  50mm 
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- Bq  =  0.2  T 

- So  =  0.4  T 

— So  =  o.6r 
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Figure  7.5:  Monodimensional  LPM 

Relative  Error  -  Kyis  for  a  =  10  S/m  -  rout  =  50mm 
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Figure  7.6:  Monodimensional  LPM 

Relative  Error  -  Kj  for  a  =  10  S/m  -  rout  =  50mm 


- fl„  =  0.2  T 

- Bq  =  0.4  T 

- Bq  =  0.6  T 

- So  =  0.8  T 


Figure  7.7:  Monodimensional  LPM 

Relative  Error  -  for  a  =  10  S/m  -  Vout  =  50mm 


Figure  7.8:  Monodimensional  LPM 

Relative  Error  -  Ky  for  a  =  10  S/m  -  rout  =  50mm 
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7.8.1  Current  Drive 
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Figure  7.9:  Monodimensional  LPM  -  Current  Drive 
r  for  (7  —  10  5/m  -  rout  =  50mm 


- So  =  0.2  T 

- So  =  0.4  T 

- So  =  0.6  7 

- So  =  0.8  T 

o  FDHM 
<  LPM -ID 


Figure  7.10:  Monodimensional  LPM  -  Current  Drive 
Tc  for  a  =  10  S/m  -  rout  =  50mm 
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Figure  7.11:  Monodimensional  LPM  -  Current  Drive 
P  for  a  =  10  S/m  -  rout  =  50mm 
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Figure  7.12:  Monodimensional  LPM  -  Current  Drive 
r  for  a  =  10  S/m  -  rout  =  50mm 
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- fl„  =  0.2  T 

— So  =  0.4  r 

- So  =  0.6  T 

- So  =  0.8  T 


Figure  7.13:  Monodimensional  LPM  -  Current  Drive 

Relative  Error  -  T  for  a  =  10  S/m  -  Tout  =  50mm 
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Figure  7.14:  Monodimensional  LPM  -  Current  Drive 

Relative  Error  -  Tc  for  a  =  10  S/m  -  rout  =  50mm 


- Bo  =  0.2  T 

- Bo  =  0.4  T 

- Bo  =  0.6  T 

- Bo  =  0.8  T 


Figure  7.15:  Monodimensional  LPM  -  Current  Drive 

Relative  Error  -  P  for  a  =  10  S/m  -  rout  =  50mm 


- Bo  =  0.2  T 

- Bo  =  0.4  T 

- Bo  =  0.6  T 

- Bo  =  0.8  T 


Figure  7.16:  Monodimensional  LPM  -  Current  Drive 

Relative  Error  -  T  for  a  =  10  S/m  -  rout  =  50mm 
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7.8.2  Voltage  Drive 
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Figure  7.17:  Monodimensional  LPM  -  Voltage  Drive 
r  for  a  —  10  S/m  -  Vout  =  50mm 
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Figure  7.18:  Monodimensional  LPM  -  Voltage  Drive 
Tc  for  a  =  10  S/m  -  rout  =  50mm 
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Figure  7.19:  Monodimensional  LPM  -  Voltage  Drive 
P  for  a  =  10  S/m  -  rout  =  50mm 
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Figure  7.20:  Monodimensional  LPM  -  Voltage  Drive 
r  for  a  =  10  S/m  -  Vout  =  50mm 
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Figure  7.21:  Monodimensional  LPM  -  Voltage  Drive 

Relative  Error  -  T  for  a  =  10  S/m  -  Tout  =  50mm 
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Figure  7.22:  Monodimensional  LPM  -  Voltage  Drive 

Relative  Error  -  Tc  for  a  =  10  S/m  -  rout  =  50mm 


- fl„  =  0.2  T 

- Bq  =  0.4  T 

- Bq  =  0.6  T 

- So  =  0.8  T 


Figure  7.23:  Monodimensional  LPM  -  Voltage  Drive 

Relative  Error  -  P  for  a  =  10  S/m  -  rout  =  50mm 


- Bo  =  0.2  T 

- Bo  =  0.4  T 

- Bo  =  0.6  T 

- Bo  =  0.8  T 


Figure  7.24:  Monodimensional  LPM  -  Voltage  Drive 

Relative  Error  -  T  for  a  =  10  S/m  -  rout  =  50mm 
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7.9  Results:  Bidimensional  LPM  -  a  =  lOS/m 
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Figure  7.25:  Bidimensional  LPM 
Kyis  for  (T  =  10  S/m  -  rout  =  50mm 
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Figure  7.26:  Bidimensional  LPM 
Kj  for  (7  =  10  S/m  -  Tout  =  50mm 
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Figure  7.27:  Bidimensional  LPM 
Cuj  for  a  =  10  S/m  -  rout  =  50mm 
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Figure  7.28:  Bidimensional  LPM 
Ky  for  a  =  10  S/m  -  rout  =  50mm 
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Figure  7.29:  Bidimensional  LPM 

Relative  Error  -  Kyis  for  a  =  10  S/m  -  rout  =  50mm 
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Figure  7.30:  Bidimensional  LPM 

Relative  Error  -  Kj  for  a  =  10  S/m  -  rout  =  50mm 


- fl„  =  0.2  T 

- Bq  =  0.4  T 

- Bq  =  0.6  T 

- So  =  0.8  T 


Figure  7.31:  Bidimensional  LPM 

Relative  Error  -  for  a  =  10  S/m  -  Vgut  =  50mm 


- Bo  =  0.2  T 

- Bo  =  0.4  T 

- Bo  =  0.6  T 

- Bo  =  0.8  T 


Figure  7.32:  Bidimensional  LPM 

Relative  Error  -  Ky  for  cr  =  10  S/m  -  rout  =  50mm 
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7.9.1  Current  Drive 
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Figure  7.33:  Bidimensional  LPM  -  Current  Drive 
r  for  a  —  10  S/m  -  rout  =  50mm 
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Figure  7.34:  Bidimensional  LPM  -  Current  Drive 
Tc  for  a  =  10  S/m  -  rout  =  50mm 
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Figure  7.35:  Bidimensional  LPM  -  Current  Drive 
P  for  a  =  10  S/m  -  rout  =  50mm 


- So  =  0.2  T 

- So  =  0.4  T 
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Figure  7.36:  Bidimensional  LPM  -  Current  Drive 
r  for  a  =  10  S/m  -  rout  =  50mm 
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Figure  7.37:  Bidimensional  LPM  -  Current  Drive 

Relative  Error  -  T  for  a  =  10  S/m  -  Tout  =  50mm 
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Figure  7.38:  Bidimensional  LPM  -  Current  Drive 

Relative  Error  -  Tc  for  a  =  10  S/m  -  rout  =  50mm 


- fl„  =  0.2  T 

- Bq  =  0.4  T 

- Bq  =  0.6  T 

- So  =  0.8  T 


Figure  7.39:  Bidimensional  LPM  -  Current  Drive 

Relative  Error  -  P  for  a  =  10  S/m  -  rout  =  50mm 
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Figure  7.40:  Bidimensional  LPM  -  Current  Drive 

Relative  Error  -  T  for  a  =  10  S/m  -  Vout  =  50mm 
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7.9.2  Voltage  Drive 
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7.41:  Bidimensional  LPM  -  Voltage  Drive 
r  for  a  —  10  S/m  -  rout  =  50mm 
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7.42:  Bidimensional  LPM  -  Voltage  Drive 
Tc  for  a  =  10  S/m  -  Tout  =  50mTO 
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7.43:  Bidimensional  LPM  -  Voltage  Drive 
P  for  a  =  10  S/m  -  rout  =  50mm 
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7.44:  Bidimensional  LPM  -  Voltage  Drive 
r  for  a  =  10  S/m  -  rout  =  50mm 
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Figure  7.45:  Bidimensional  LPM  -  Voltage  Drive 

Relative  Error  -  T  for  cr  =  10  S/m  -  rout  =  50mm 
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Figure  7.46:  Bidimensional  LPM  -  Voltage  Drive 

Relative  Error  -  Tc  for  a  =  10  S/m  -  rout  =  50mm 


- fl„  =  0.2  T 

- Bq  =  0.4  T 

- Bq  =  0.6  T 

- So  =  0.8  T 


Figure  7.47:  Bidimensional  LPM  -  Voltage  Drive 

Relative  Error  -  P  for  a  =  10  S/m  -  Vout  =  50mm 


h  =  30mm 


- So  =  0.2  T 
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- So  =  0.6  T 

- So  =  0.8  T 


Figure  7.48:  Bidimensional  LPM  -  Voltage  Drive 

Relative  Error  -  T  for  a  =  10  S/m  -  rout  =  50mm 
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7.10  Results:  Monodimensional  LPM  -  a  =  10® /S'/m 
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Figure  7.49:  K^is 
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Figure  7.50:  Monodimensional  LPM 
Kj  for  a  —  10®  Sjm  -  Tout  =  50mm 
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Figure  7.51:  Monodimensional  LPM 
Clu  for  a  =  10®  5/m  -  rout  =  50mm 
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Figure  7.52:  Monodimensional  LPM 
Kv  for  a  =  10®  S/m  -  Tout  =  50mm 
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Figure  7.53:  Monodimensional  LPM 

Relative  Error  -  Kyis  for  a  =  10®  S/m  -  rout  =  50mm 


Figure  7.54:  Monodimensional  LPM 

Relative  Error  -  Kj  for  a  =  10®  S/m  -  rout  =  50mm 


Figure  7.55:  Monodimensional  LPM 

Relative  Error  -  c^y  for  a  =  10®  S/m  -  rout  =  50mm 
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Figure  7.56:  Monodimensional  LPM 

Relative  Error  -  Ky  for  a  =  10®  S/m  -  rout  =  50mm 
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7.10.1  Current  Drive 
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Figure  7.57:  Monodimensional  LPM  -  Current  Drive 
r  for  a  —  10^  S/m  -  rout  =  50mm 
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Figure  7.58:  Monodimensional  LPM  -  Current  Drive 
Tc  for  a  =  10®  S/m  -  rout  =  50mm 
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Figure  7.59:  Monodimensional  LPM  -  Current  Drive 
P  for  a  =  10®  S/m  -  Tout  =  50mm 
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Figure  7.60:  Monodimensional  LPM  -  Current  Drive 
r  for  a  =  10®  S/m  -  rout  =  50mm 
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Figure  7.61:  Monodimensional  LPM  -  Current  Drive 

Relative  Error  -  T  for  a  =  10®  S/m  -  rout  =  50mm 
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Figure  7.62:  Monodimensional  LPM  -  Current  Drive 

Relative  Error  -  Tc  for  a  =  10®  S/m  -  rout  =  50mm 
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Figure  7.63:  Monodimensional  LPM  -  Current  Drive 

Relative  Error  -  P  for  a  =  10®  S/m  -  rout  =  50mm 
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Figure  7.64:  Monodimensional  LPM  -  Current  Drive 

Relative  Error  -  F  for  a  =  10®  S/m  -  rout  =  50mm 
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7.10.2  Voltage  Drive 
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Figure  7.65:  Monodimensional  LPM  -  Voltage  Drive 
r  for  a  =  10®  S/m  -  Vout  =  50mm 


h  =  IO777777 


h  =  2O777777 


1  2  3  4  5 

L  —  [mm] 


2  3  4  5 

L  —  [mml 


h  =  3O777777 


- So  =  0.2  T 

- So  =  0.4  T 

- So  =  0.6r 

- So  =  0.8  T 

o  FDHM 
<  LPM -ID 


Figure  7.66:  Monodimensional  LPM  -  Voltage  Drive 
Tc  for  a  =  10®  S/m  -  rout  =  50mm 


X  10"^  ^  “  IO777777 


L  —  [mml 


X  10^  ^  “  2O777777 


L  —  [mml 


X  10*^  h  =  30777 777 


L  —  [mml 


- So  =  0.2  T 

- So  =  0.4  T 

- So  =  0.6  T 

- So  =  0.8  T 

o  FDHM 
<  LPM -ID 


Figure  7.67:  Monodimensional  LPM  -  Voltage  Drive 
P  for  a  =  10^  S/m  -  rout  =  50mm 


- So  =  0.2  T 

- So  =  0.4  T 

- So  =  0.6  T 

- So  =  0.8  T 

o  FDHM 
<  LPM -ID 


Figure  7.68:  Monodimensional  LPM  -  Voltage  Drive 
r  for  a  =  10®  S/m  -  rout  =  50mm 
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- fl„  =  0.2  T 

— So  =  0.4  r 

- So  =  0.6  T 

- So  =  0.8  T 


Figure  7.69:  Monodimensional  LPM  -  Voltage  Drive 

Relative  Error  -  T  for  a  =  10®  S/m  -  rout  =  50mm 
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Figure  7.70:  Monodimensional  LPM  -  Voltage  Drive 

Relative  Error  -  Tc  for  a  =  10®  S/m  -  Vout  =  50mm 
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Figure  7.71:  Monodimensional  LPM  -  Voltage  Drive 

Relative  Error  -  P  for  a  =  10^  S(m  -  Vout  =  50mm 


- So  =  0.2  T 

- So  =  0.4  T 

- So  =  0.6  T 

- So  =  0.8  T 


Figure  7.72:  Monodimensional  LPM  -  Voltage  Drive 

Relative  Error  -  T  for  a  =  10^  S/m  -  rout  =  50mm 
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7.11  Results:  Bidimensional  LPM  -  a  = 


Figure  7.73:  Bidimensional  LPM 
Kyis  for  a  =  10®  S/m  -  Tout  =  50mm 
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Figure  7.74:  Bidimensional  LPM 
Kj  for  a  =  10®  S/m  -  rout  =  50mm 
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Figure  7.75:  Bidimensional  LPM 
Cuj  for  a  =  10®  S/m  -  rout  =  50mm 

X 10“^  ^  “  10mm  ^  jq-4  h  =  20mm  ^  jq-4  h  =  30mm 


- Bo  =  0.2  T 
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Figure  7.76:  Bidimensional  LPM 
Kv  for  a  =  10®  S/m  -  Tout  =  50mm 
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- fl„  =  0.2  T 

— So  =  0.4  r 

- So  =  0.6  T 

- So  =  0.8  T 


Figure  7.77:  Bidimensional  LPM 

Relative  Error  -  Kyis  for  a  =  10®  S/m  -  rout  =  50mm 


Figure  7.78:  Bidimensional  LPM 

Relative  Error  -  Kj  for  a  =  10®  S/m  -  rout  =  50mm 


Figure  7.79:  Bidimensional  LPM 

Relative  Error  -  for  a  =  10®  S/m  -  rout  =  50mm 
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Figure  7.80:  Bidimensional  LPM 

Relative  Error  -  Ky  for  a  =  10®  S/m  -  rout  =  50mm 
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7.11.1  Current  Drive 
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Figure  7.81:  Bidimensional  LPM  -  Current  Drive 
r  for  a  =  10®  S/m  -  Vout  =  50mm 


Figure  7.82:  Bidimensional  LPM  -  Current  Drive 
Tc  for  a  =  10®  S/m  -  r^ut  =  50mm 
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Figure  7.83:  Bidimensional  LPM  -  Current  Drive 
P  for  a  =  10®  S/m  -  rout  =  50mm 
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Figure  7.84:  Bidimensional  LPM  -  Current  Drive 
r  for  a  =  10®  S/m  -  rout  =  50mm 
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So  =  0.6  T 
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Figure  7.85:  Bidimensional  LPM  -  Current  Drive 

Relative  Error  -  T  for  a  =  10®  S/m  -  rout  =  50mm 
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Figure  7.86:  Bidimensional  LPM  -  Current  Drive 

Relative  Error  -  Tc  for  a  =  10®  S/m  -  rout  =  50mm 


Figure  7.87:  Bidimensional  LPM  -  Current  Drive 

Relative  Error  -  P  for  a  =  10®  S/m  -  rout  =  50mm 
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Figure  7.88:  Bidimensional  LPM  -  Current  Drive 

Relative  Error  -  F  for  a  =  10®  S/m  -  rout  =  50mm 
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7.11.2  Voltage  Drive 
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Figure  7.89:  Bidimensional  LPM  -  Voltage  Drive 
r  for  a  =  10®  S/m  -  Vout  =  50mm 
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Figure  7.90:  Bidimensional  LPM  -  Voltage  Drive 
Tc  for  a  =  10®  S/m  -  rout  =  50mm 
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Figure  7.91:  Bidimensional  LPM  -  Voltage  Drive 
P  for  a  =  10®  S/m  -  rout  =  50mm 
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Figure  7.92:  Bidimensional  LPM  -  Voltage  Drive 
r  for  a  =  10®  S/m  -  Vout  =  50mm 
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Figure  7.93:  Bidimensional  LPM  -  Voltage  Drive 

Relative  Error  -  T  for  a  =  10®  S/m  -  rout  =  50mm 
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Figure  7.94:  Bidimensional  LPM  -  Voltage  Drive 

Relative  Error  -  Tc  for  a  =  10®  S/m  -  rout  =  50mm 
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Figure  7.95:  Bidimensional  LPM  -  Voltage  Drive 

Relative  Error  -  P  for  a  =  10®  S/m  -  rout  =  50mm 
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Figure  7.96:  Bidimensional  LPM  -  Voltage  Drive 

Relative  Error  -  V  for  a  =  10®  S/m  -  rout  =  50mm 
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7.12  Discussion 

Figures  (7.1)  and  (7.5)  show  the  inability  of  the  mono-dimensional  analytical 
solution  to  give  good  results  in  the  estimation  of  K^is  if  the  electric  conductiv¬ 
ity  of  the  liquid  is  low.  The  reason  resides  on  the  fact  that  the  Shercliff  and 
Hartmann  layers  are  of  the  same  magnitude  if  the  electrical  conductivity  of  the 
liquid  is  low.  So,  as  the  mono-dimensional  solution  neglects  the  viscous  shear 
on  the  armatures  of  the  device,  the  computed  K^s  is  underestimated. 

On  the  other  hand,  in  Figures  (7.49)  and  (7.53)  the  mono-dimensional  solution 
shows  its  ability  to  compute  the  Kyis  for  high  conductive  liquid;  the  reason  of 
the  improvement  is  in  the  fact  that  the  Hartmann  layers  are  a  bit  less  than  2 
order  of  magnitude  thinner  than  the  Shercliff  layers,  so  the  neglection  of  the 
influence  of  the  Shercliff  layers  does  not  affect  the  estimation  of  K^ig. 

For  both  high  and  low  electrically  conductive  liquids  the  error  in  the  estimation 
of  Kyis  decreases  for  increasing  height  and  decreasing  radial  dimension  of  the 
torus  and  for  decreasing  magnetic  held. 

The  coefficient  K^s  computed  on  the  basis  of  the  bi-dimensional  solution  shows 
better  accuracy  for  low  electrically  conductive  liquids  than  for  high  electrically 
conductive  liquids.  For  low  electrically  conductive  liquids,  the  accuracy  of  the 
Kyis  based  on  the  bi-dimensional  solution  is  furthermore  more  accurate  than 
the  one  based  on  the  mono-dimensional  solution;  the  opposite  trend  can  be 
read  for  high  electrically  conductive  liquids.  The  error  in  the  estimation  of 
Kyis  increases  for  decreasing  height  and  decreasing  radial  dimension  of  the 
torus  and  for  decreasing  magnetic  held. 

The  computed  coefficient  of  Lorentz  force  moment,  Kj,  is  directly  inhuenced 
neither  by  the  analytical  solution  used  nor  by  the  electric  conductivity  of  the 
liquid.  Figures  (7.2)-(7.54)  show  actually  the  precision  of  the  FDHM  model, 
being  the  Kj  given  in  closed  form. 

The  errors  regarding  the  coefficient  Kj  are  not  directly  ahected  by  the  partic¬ 
ular  analytical  solution  used. 

The  neglection  of  the  Sherclih  layers  does  not  ahect  the  estimation  of  for 
both  high  and  low  electrically  conductive  liquids.  The  reason  resides  in  the 
presence  of  the  not  derived  value  of  the  angular  velocity  in  the  formula  for  the 
computation  of  c^j  and  it  seems  to  be  quite  accurate  for  both  high  and  low 
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electrically  conductive  liquids.  The  coefficient  is  anyway  overestimated  due  to 
the  neglection  of  the  Shercliff  layers.  The  error  in  the  estimation  of  increases 
for  decreasing  height  and  increasing  radial  dimension  of  the  torus  and  for  de¬ 
creasing  magnetic  held. 

Figures  (7.27)-(7.79)  show  that  the  estimation  of  based  on  the  bi-dimensional 
solution  is  more  precise  for  high  electrically  conductive  liquid  than  for  low  elec¬ 
trically  conductive  liquid,  so  following  the  same  trend  of  the  based  on  the 
mono-dimensional  solution.  The  error  in  this  case  is  for  both  high  and  low 
electrically  conductive  liquids  higher  than  the  estimated  with  the  mono¬ 
dimensional  solution. 

As  for  the  computation  of  the  c^;,  here  the  neglection  of  the  Shercliff  layers 
does  not  affect  the  estimation  of  Ky  for  both  high  and  low  electrically  conduc¬ 
tive  liquids.  The  coefficient  is  anyway  overestimated  due  to  the  neglection  of 
the  Shercliff  layers.  The  error  in  the  estimation  of  Ky  increase  for  decreasing 
height  and  increasing  radial  dimension  of  the  torus  and  for  decreasing  magnetic 
held.  The  error  in  the  estimation  of  Ky  increases  for  decreasing  height  and 
increasing  radial  dimension  of  the  torus  and  for  decreasing  magnetic  held. 

The  Ky  based  on  the  bi-dimensional  solution  has  accuracy  comparable  with 
the  one  computed  with  the  mono-dimensional  solution  with  same  trend  with 
respect  of  the  height  and  radial  dimension  of  the  torus  and  the  intensity  of  the 
magnetic  held. 

The  precision  torque  T  computed  with  both  the  monodimensional  and  bidimen- 
sional  solutions  is  mainly  ahected  by  the  estimation  error  of  the  Kvis  being  the 
Cui  and  Kj  quite  accurate.  The  monodimensional  solution  is  able  to  computed 
with  sufficient  accuracy  the  torque  provided  by  the  device  in  the  case  a  high 
conductive  liquid  is  used,  the  bidimensional  solution  guarantees  good  perfor¬ 
mances  for  low  conductive  liquids.  The  comparison  between  the  predictions  in 
the  case  of  current  and  voltage  drive  shows  that  The  torque  T  for  Voltage  Drive 
is  more  accurate  than  the  one  computed  for  the  Current  Drive.  The  reason 
resides  in  the  weight  the  K^s  has  in  in  (7.7.8).  The  error  of  the  K^s  is  reduced 
by  the  low  value  of  the  electrical  resistance  of  the  fluid  and  so,  the  estimated 
torque  appears  to  be  close  to  the  one  derived  from  the  FDHM  simulations. 
About  the  characteristic  time  it  can  be  seen  that  the  error  between  the  char¬ 
acteristic  time  computed  by  the  FDHM  and  the  LPM  is  not  the  result  of  the 
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propagation  error  of  each  coefficient  used  to  compute  Tc-  A  further  analysis 
of  the  time  evolution  of  the  coefficients  and  K^is  made  with  the  results  ob¬ 
tained  with  the  FDHM,  shows  how  these  coefficients,  considered  constant  in 
the  LPM,  are  function  of  the  velocity,  Figures(7.97)  (7.98)  (7.99). 


FDHM  -  c^(t) 


t-[s] 


Figure  7.97:  Evolution  of 


FDHM  - 


Figure  7.98:  Evolution  of  Kyis 


This  analysis  highlights  how  the  device  is  affected  by  non-linear  effect  in  the 
neighbourhood  of  null  angular  velocity.  This,  in  some  way,  resize  the  hypo¬ 
thetical  advantage  not  to  have  jump  of  torque  around  the  zero  crossing  [3]. 
Furthermore,  it  is  evident  from  the  previous  results  that  the  characteristic  time 
is  smaller  in  case  of  Voltage  drive  because  of  the  role  played  by  the  electrical 
resistance,  supposed  to  be  little  for  this  kind  of  device. 
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FDHM  -  uia{t) 


Figure  7.99:  Evolution  of  ujjju 


The  error  on  the  characteristic  time  influences  the  estimated  power  con¬ 
sumption  and  the  estimated  angular  momentum. 


123 


DISTRIBUTION  A.  Approved  for  public  release:  distribution  unlimited. 


Chapter 


8 


Performances  analysis. 


In  this  chapter,  a  performance  analysis  has  been  conducted  by  means  of  the 
LPM.  The  possibility  to  use  liquid  metal  or  an  ideal  liquid  to  be  designed 
as  been  evaluated. 


8.1  Liquid  Metal .  125 

8.2  Ideal  Liquid  .  127 

8.3  Comparison  with  commercial  reaction  wheels  designed  for  Cubesatsl2^ 

8.4  Conclusions .  131 
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8.1  Liquid  Metal 

Considering  mercury  as  conductive  fluid  with  the  following  physical  properties: 


aHg  =  1-04  ■  10®  S/m 

PHg  =  1.36  ■  10^  Kg/m^  (S-l-l) 

PHg  =  1.50  ■  10-3  Pa/s 


the  coefficients  of  the  lumped  parameter  model  together  with  the  torque  and 
the  power  consumption  are  presented  in  Figures(8.1)-(8.4).  All  the  values 
plotted  in  the  following  figure  are  calculated  for  a  torus  with  rectangular  cross 
section  whose  height  is  equal  io  h  =  1cm.  The  symbols  on  the  y  axis  of 
Fig. (8.3)  and  Fig. (8.4)  have  the  following  meaning: 


T  = 


f 

ah 


(8.1.2) 


P  = 


P 

hJ^ 


1 

h 


KyKi 
K  ■ 


(8.1.3) 


where  P  is  the  power  consumption  of  the  actuator  for  an  electric  current  with 
value  of  1  A.  All  the  following  values  are  obtained  by  means  of  the  LPM  based 
on  the  mono- dimensional  solution  with  the  value  of  K^is  computed  by  the 
bi-dimensional  solution. 


-4  7*e  =  25mm 


-4  Te  =  50mm 
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4  6  8 

L  —  [mm] 


4  6  8 
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Figure  8.1:  Kg  function  of  L  for  Bq  =  0.2  1.0  and  3  different  outer  radii 

re  =  25,50,  lOOmm.  Liquid:  Mercury  (Hg). 
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Figure  8.2:  Kyis  function  of  L  for  Bq  =  0.2  1.0  and  3  different  outer  radii 

Te  =  25,  50,  lOOmm.  Liquid:  Mercury  (Hg). 


-4  Te  =  25mm 
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- B„=1T 


Figure  8.3:  Torque  function  of  L  for  Bq  —  0.2  -^1.0  and  3  different  outer  radii 
Te  =  25,50,100mm.  Liquid:  Mercury  (Hg). 


- B„=0.2T 

- B„=0.4T 

- B„=0.6T 

— 

- B„=1T 


Figure  8.4:  Power  consumption  function  of  L  for  Bq  =  0.2  1.0  and  3  differ¬ 

ent  outer  radii  re  =  25,  50, 100mm.  Liquid:  Mercury  (Hg). 
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=  25mm  =  50mm 


Tc  =  100mm 


Figure  8.5:  Characteristic  time  function  of  L  for  cr  =  10  10®  S/m  and 

3  different  outer  radii  r^  =  25, 50,  lOOmm.  Bq  =  1.0  T. 
p  =  1000  Kg/m^ . 


In  Figure(8.2)  it  is  possible  to  see  that  the  influence  of  the  viscosity  is  not 
neglectable  and  it  influences  the  overall  performance  of  the  device.  The  high 
value  of  the  K^is  is  proportional  to  the  electrical  conductivity  of  the  fluid  due 
to  the  coupling  between  the  induced  voltage  and  the  velocity  profile. 

The  high  values  of  compared  to  the  values  of  Kj  lead  to  a  low  level  of 
torque  provided  by  the  actuator.  Furthermore  the  torque  is  limited  in  time  by 
the  total  current  that  can  be  fed  by  the  power  subsystem. 


8.2  Ideal  Liquid 


In  order  to  understand  the  influence  of  the  electric  conductivity  on  the  per¬ 
formances  of  the  device,  the  values  of  K^jis  have  been  computed  for  different 
electrical  conductivities  of  the  liquid  together  with  the  power  consumption  P 
and  the  generated  torque  T.  The  values  are  calculated  for  a  liquid  with  den¬ 
sity  p  =  1000  Kg/mf’  and  viscosity  p  =  10“^  Pa  s;  the  magnetic  held  is  set  to 
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Figure  8.6:  Kyis  function  of  L  for  ct  =  10  10®  S/m  and  3  different  outer 

radii  r^  =  25,  50, 100mm.  Bq  =  1.0  T.  p  =  1000  Kg/m^ . 
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Figure  8.7:  Torque  function  of  L  for  a  =  10  10®  S/m  and  3  different  outer 

radii  r^  =  25,  50, 100mm.  Bq  =  1.0  T.  p  =  1000  Kg/m? . 
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Figure  8.8:  Power  consumption  function  of  L  for  a  =  10  10®  S/m  and 

3  different  outer  radii  r^  =  25, 50, 100mm.  Bq  =  1.0  T. 
p  =  1000  Kg/m^ . 
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Figure  8.9:  Characteristic  time  function  of  L  for  ct  =  10  10®  S/m  and 

3  different  outer  radii  rg  =  25, 50, 100mm.  Bo  =  1.0T. 
p  =  1000  Kg/m^ . 


From  Figures(8.6)  and  (8.7)  it  is  possible  to  see  the  improving  in  the  perfor¬ 
mances  for  low  condncting  flnid.  Dne  to  the  low  conpling  between  the  velocity 
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and  the  magnetic  field,  the  generated  torque  for  a  low  conductive  fluid  be¬ 
come  bigger  of  about  two  orders  of  magnitude  with  respect  to  an  high  conduc¬ 
tive  fluid.  As  expected,  the  power  consumption  has  an  opposite  trend  but  its 
maximum  value  is  anyway  around  IW  if  an  electric  current  of  lA  is  applied, 
Figure(8.8). 

For  high  conductive  fluids  the  estimation  of  the  parameters  is  expected  to  be 
quite  accurate  due  to  the  fact  that  the  length  of  the  boundary  layers  on  the 
upper  and  bottom  walls  scale  as  6sh  ~  l/\/ Ha  while  the  boundary  layers  on 
the  lateral  surfaces  scale  as  Sna  ~  1/Ha.  For  high  Ha  these  two  characteristic 
length  can  differ  of  one  order  of  magnitude,  then  the  error  due  to  the  neglection 
of  the  Shercliff  layers,  next  to  the  armatures,  is  low. 

On  the  other  side  the  prediction  of  and  the  available  torque  could  be  re¬ 
garded  as  not  accurate  for  very  low  conductive  fluid  due  to  the  reason  that  the 
Shercliff  layers  on  the  upper  and  bottom  surface  of  the  torus  are  in  this  case  of 
the  same  order  of  magnitude  as  the  Hartmann  layers  on  the  lateral  walls  of  the 
torus.  Anyway  the  density  of  an  hypothetical  low  conductive  fluid  is  expected 
to  be  low  if  compared  with  the  density  of  a  liquid  metal.  This  implies  that 
in  order  to  have  the  same  moment  of  inertia  for  low  conductive  fluid,  we  need 
to  increase  the  height  of  the  cross  section  if  the  working  fluid  has  low  electric 
conductivity.  This  implies  again  that  the  viscous  shear  on  the  top  and  bottom 
boundaries  could  be  neglectable  with  respect  of  the  viscous  shear  acting  on  the 
lateral  cylindrical  surfaces  of  the  device. 

8.3  Comparison  with  commercial  reaction  wheels 
designed  for  Cubesats 

In  this  section,  in  order  exploits  the  performances  of  the  proposed  Liquid 
Wheel,  a  comparison  with  three  commercial  reaction  wheels  designed  for  Cube¬ 
sats  will  be  presented.  The  sizing  of  the  Liquid  Wheel  has  been  made  with 
the  purpose  to  have  a  final  mass  of  conductive  liquid  comparable  with  the 
overall  mass  of  the  commercial  classical  reaction  wheels  taken  as  example. 
Furthermore  the  dimensions  of  the  simulated  Liquid  Wheel  have  been  made 
suitable  for  the  Cubesat  Standard.  The  three  reaction  wheels  are  available 
at  www.isispace.nl  and  represent  three  different  classes  of  reaction  wheels  ac- 
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cording  to  the  torque  provided.  The  selected  working  fluid  for  the  simulated 
Liquid  Wheel  is  Mercury,  even  if  the  advantage  to  use  an  ad-hoc  conductive 
liquid  has  been  commented  and  proved  in  the  dedicated  section  of  this  chapter. 
The  comparison  is  shown  in  the  table  8.1. 


CubeWheel  S 

CubeWheel  M 

CubeWheel  L  I 

Weight 

55  g 

130  g 

200  g 

Lenght 

28  mm 

46  mm 

57  mm  | 

Width 

31  mm 

46  mm 

57  mm 

Height 

26  mm 

31.5  mm 

31.5  mm  | 

Torque 

0.23  mNm 

1.0  mNm 

2.3  mNm 

Angular  Momentum 

1.7  mNms 

10  mNms 

30  mNms  | 

Power 

<  600  mW 

<  1000  mW 

<  2200  mW 

MHD  Wheel  S 

MHD  Wheel  M 

MHD  Wheel  L 

Weight  (Liquid) 

64  g 

124  g 

181  g 

Outer  Radius 

45  mm 

45  mm 

45  mm  | 

Inner  Radius 

43  mm 

41  mm 

39  mm 

Height 

1  10  mm 

10  mm 

10  mm  1 

Torque(@lA/s) 

0.18  mNm 

0.36  mNm 

0.53  mNm 

Angular  Momentum((§lA) 

1  0.18  mNms 

0.36  mNms 

0.53  mNms  | 

Power((§lA) 

0.55  mW 

0.54  mW 

0.52  mW 

Table  8.1:  Comparison  between  elassical  Reaetion  Wheels  for  Cubesats  and 
Liquid  Wheels 


The  values  listed  in  table  8.1  show  that  the  proposed  actuator,  even  provid¬ 
ing  a  torque  slightly  smaller  than  the  classical  reaction  wheels,  is  characterized 
by  a  power  consumption  that  is  three  orders  of  magnitude  less  than  the  classi¬ 
cal  reaction  wheels.  The  power  consumption  proposed  do  not  take  into  account 
the  efficiency  of  the  power  board  that  have  to  drive  the  actuator,  but  shows 
only  the  power  consumption  of  the  Liquid  Wheel  across  its  electric  armatures. 
The  angular  momentum  provided  by  the  Liquid  Wheel  is  computed  for  a  cur¬ 
rent  of  lA.  In  general  the  possibility  to  increase  the  current  supplied  to  the 
actuator,  leads  to  a  proportional  increase  of  the  provided  angular  momentum. 
But  this  aspect  is  constrained  by  the  capability  of  the  power  board. 

What  is  evident  in  the  numbers  proposed  is  that  the  Liquid  Wheel  can  provide 
comparable  performances  of  small  reaction  wheels,  similar  to  the  ones  in  table 
8.1.  As  the  required  torque  or  angular  momentum  increase  the  performance  of 
the  classical  reaction  become  hard  to  meet  by  a  Liquid  Wheel. 

Even  if  in  paragraph  7.12  it  has  been  shown  that  the  behaviour  of  the  Liquid 
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wheel  is  characterized  by  non-linearities  in  the  vicinity  of  the  null  angular  veloc¬ 
ity,  the  Liquid  Wheel  in  this  zone  is  not  characterized  by  discontinue  torque, 
as  for  the  classical  reaction  wheels.  This  aspect,  together  to  the  possibility 
to  provide  a  clean  torque,  not  influenced  by  ripple  or  mechanical  noise,  still 
represent  positive  aspects  of  the  proposed  actuator. 

8.4  Conclusions 

In  this  chapter  the  possibility  to  use  an  ideal  electrically  conductive  liquid 
with  low  electric  conductivity  is  evaluated.  The  results  show  that  this  can  be 
a  good  solution  in  order  to  improve  the  performances  of  the  device  in  terms 
of  provided  torque.  Unfortunately  such  kind  of  liquid  is  not  yet  available.  It 
should  be  the  a  non  aqueous  solution  of  electrolytes  in  order  to  avoid  the  the 
electrolysis  phenomena  and  so  the  generation  of  hydrogen  {H2)  and  oxygen 
(O). 
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Conclusions 


A  novel  actuator  for  spacecraft  attitude  control  with  liquid  flywheel  has  been 
studied.  Two  different  conhgurations  of  the  device  with  didtributed  magnetic 
held  have  been  studied  on  the  base  of  the  optimization  of  dimensionless  mo¬ 
ment  of  inertia  and  the  minimization  of  the  viscous  shear. 

In  order  to  estimate  the  performances  of  this  actuator,  a  2-dimensional  Finite 
Differences  Hybrid  Model  (FDHM)  has  been  developed  on  the  base  of  the  MHD 
set  equation  under  the  hypothesis  of  low  Magnetic  Reynolds.  The  model  solves 
numerically  the  time  dependent  axially  symmetric  problem  of  a  conductive  liq¬ 
uid  rotating  in  a  torus  with  rectangular  cross  section  due  to  the  interaction 
of  a  radial  magnetic  held  and  an  axial  electric  held.  The  electric  side  of  the 
problem  has  been  solved  by  means  of  the  node  method  applied  to  a  network  of 
electric  resistances  and  voltage  generators  representing  the  back  electromotive 
voltage  induced  by  the  spinning  liquid  through  the  magnetic  held.  The  huid- 
dynamics  side  of  the  problem  has  been  solved  using  a  Crank-Nicolson  method 
over  a  non  uniform  and  collocated  grid.  The  grid  generator  has  been  written 
to  be  sensitive  to  the  Hartmann  number  of  the  problem. 

The  FDHM  has  shown  great  accuracy  over  a  wide  range  of  geometrical  and 
physical  quantities  and  the  analysis  of  the  results  has  helped  a  deeper  under¬ 
standing  of  the  device  and  its  performances  related  to  the  dimension  of  the 
torus,  the  intensity  of  the  magnetic  held,  the  conductivity  of  the  working  huid 
and  the  diherent  drive  chosen. 

A  Lumped  Parameter  Model  (LPM)  has  been  then  derived  using  two  station¬ 
ary  mono  and  bi-dimensional  analytical  solution  of  the  MHD  set  of  equations 
under  the  hypothesis  of  low  Magnetic  Reynolds.  The  LPM  has  been  compared 
with  the  results  obtained  with  the  FDHM,  showing  good  agreement  in  the  es¬ 
timation  of  the  coefficients  of  the  model. 

The  possibility  to  use  an  ideal  huid  diherent  from  liquid  metal  has  been  evalu¬ 
ated  by  means  of  the  LPM.  Due  to  the  low  coupling  between  the  velocity  and 
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the  magnetic  field,  the  generated  torque  of  a  device  using  a  low  conductive  fluid 
is  about  two  orders  of  magnitude  bigger  than  a  device  using  a  high  conductive 
fluid.  As  expected,  the  power  consumption  has  an  opposite  trend. 
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A 


The  Crank-Nicolson  Method 


The  Crank-Nicolson  method  is  a  well-known  hnite  difference  method  for 
the  numerical  integration  of  the  heat  equation  and  closely  related  partial  dif¬ 
ferential  equations. 

It  is  an  implicit  second-order  method  in  time,  it  can  be  written  as  an  implicit 
Runge-Kutta  method,  and  it  is  numerically  stable.  The  method  was  developed 
by  John  Crank  and  Phyllis  Nicolson  in  the  mid  20th  century. 

Considerind  a  ID  partial  differential  equation: 


du  f  du  d‘^u\ 


(A.0.1) 


the  time  derivative  on  the  left  hand  side  can  be  replaced  with  a  forward  finite 
difference: 


du  _  -  u] 

dt  At 

j,n 


(A.0.2) 


The  Crank-Nicolson  method  appears  to  be  a  combination  of  the  forwards  and 
backwards  Euler  method  at  two  different  steps  in  time,  n  an  n  -|-  1. 


Ml 


At 


1 

2 


d‘^u\  /  du  d‘^u\ 


(A.0.3) 
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j-l,n+l  j,n+l  ]+l,ml 


Figure  A.l:  Crank- Nicols  on  stencil  for  a  ID  problem 


As  the  quantities  defined  at  the  time  step  n+1  are  unknown,  in  order  to  update 
the  solution  at  the  time  n  +  1,  as  algebraic  system  must  be  solved: 

^n+i  (A.0.4) 

The  shape  of  the  matrix  A  depends  on  the  studied  partial  differential  equation 
and  on  the  chosen  finite  difference  scheme  used  to  approximate  the  spatial  part 
of  the  partial  differential  equation.  The  vector  B  is  not  null  if  the  source  term 
of  the  differential  equation  is  not  equal  to  zero. 
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Results:  FDHM  for  a  =  10^  S/m 


B.l  Current  Drive .  137 

B.2  Voltage  Drive  .  139 
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B.l  Current  Drive 
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Figure  B.l:  Current  Drive  -  Torque  for  a  =  10^  S/m  -  FDHM 
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Figure  B.2:  Current  Drive  -  Characteristic  time  for  a  =  10^  S/m  -  FDHM 
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Figure  B.3:  Current  Drive  -  Specific  Power  for  a  =  10^  S/m  -  FDHM 
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Figure  B.4:  Current  Drive  -  Specific  Angular  Momentum  for  a  =  10^  S/m  - 
FDHM 
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Figure  B.5:  Current  Drive  -  Relative  Dynamic  Error  for  a  =  10^ 
FDHM 
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Figure  B.6:  Current  Drive  -  Standard  deviation  of  the  electric  current  for 
a  =  10^  S/m  -  FDHM 


138 


DISTRIBUTION  A.  Approved  for  public  release:  distribution  unlimited. 


App.  B  Results:  FDHM  for  a  =  10^  S/m 


B.2  Voltage  Drive 
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Figure  B.7:  Voltage  Drive  -  Torque  for  a  =  10^  S/m  -  FDHM 
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Figure  B.8:  Voltage  Drive  -  Characteristic  time  for  a  =  10^  S/m  -  FDHM 
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Figure  B.9:  Voltage  Drive  -  Specific  Power  for  a  =  10^  S/m  -  FDHM 
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App.  B  Results:  FDHM  for  a  =  10^  S/m 


a  =  10^5/m  —  50mm 


X 10”^  ^  “  10mm 


L  —  [mm] 


^  j^-4  h  =  20mm 


L  —  [mm] 


X 10“^  ^  “  30mm 


L  —  [mm] 


- So  =  0.2T 


- So  =  0.4T 


- So  =  0.6T 


- So  =  0.8T 


Figure  B.IO:  Voltage  Drive  -  Specific  Angular  Momentum  for  a  =  10^  S/m  - 
FDHM 
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Figure  B.ll:  Voltage  Drive  -  Relative  Dynamic  Error  for  a  =  10^  S/m  - 
FDHM 
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Figure  B.12:  Voltage  Drive  -  Standard  deviation  of  the  electric  current  for 
a  =  10^  S/m  -  FDHM 


140 


DISTRIBUTION  A.  Approved  for  public  release:  distribution  unlimited. 


Appendix 


C 


Results:  FDHM  for  a  =  10^  S/m 


C.  1  Current  Drive .  142 
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App.  C  Results:  FDHM  for  a  =  10^  S/m 
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App.  C  Results:  FDHM  for  a  =  10^  S/m 
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145 


DISTRIBUTION  A.  Approved  for  public  release:  distribution  unlimited. 


Appendix 


D 


Results:  FDHM  for  cr  =  10"^  S/m 


D.  1  Current  Drive .  147 

D.2  Voltage  Drive  .  149 


146 


DISTRIBUTION  A.  Approved  for  public  release:  distribution  unlimited. 


App.  D  Results:  FDHM  for  a  =  10'^  S/m 
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Figure  D.l:  Current  Drive  -  Torque  for  a  =  10^  S/m  -  FDHM 
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Figure  D.3:  Current  Drive  -  Specific  Power  for  a  =  10^  S/m  -  FDHM 
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App.  D  Results:  FDHM  for  a  =  10'^  S/m 
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App.  D  Results:  FDHM  for  a  =  10'^  S/m 
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Figure  D.9:  Voltage  Drive  -  Specific  Power  for  a  =  10^  S/m  -  FDHM 
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App.  E  Results:  FDHM  for  a  =  10®  S/m 
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App.  E  Results:  FDHM  for  a  =  10®  S/m 
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Figure  E.5:  Current  Drive  -  Relative  Dynamic  Error  for  a  =  10®  S/m  - 
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App.  E  Results:  FDHM  for  a  =  10®  S/m 
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App.  E  Results:  FDHM  for  a  =  10®  S/m 
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App.  F  Results:  Comparison  between  FDHM  and  LPM:  a  =  10^  S/m 
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App.  F  Results:  Comparison  between  FDHM  and  LPM:  a  =  10^  S/m 
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App.  F  Results:  Comparison  between  FDHM  and  LPM:  a  =  10^  S/m 
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X 10'"^  ^  “  10mm 


L  —  \mm\ 
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X 10^^  =  30mm 


L  —  \mm] 
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Figure  F.12:  Monodimensional  LPM  -  Current  Drive 
r  for  a  =  10^  S/m  -  Vout  =  50mm 
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Errrd  -  f  ,  6  ,  ?  ,  6  ,  Ewrei  -  Torque 


App.  F  Results:  Comparison  between  FDHM  and  LPM:  a  =  10^  S/m 


- Bq  =  0.2  T 

- So  =  0.4  T 

- So  =  0.6  T 

- So  =  0.8  T 


Figure  F.13:  Monodimensional  LPM  -  Current  Drive 

Relative  Error  -  F  for  a  =  10^  S/m  -  rout  =  50mm 


- fl„  =  0.2  T 

- Bq  =  0.4  T 

- flo  =  0.6T 

- So  =  0.8  T 


Figure  F.14:  Monodimensional  LPM  -  Current  Drive 

Relative  Error  -  To  for  a  =  10^  S/m  -  rout  =  50mm 


h  =  10777,771 


X  10"^  ^  “  20777771 


L  —  [mm] 


X  10’^  ft-  =  30777777 


- So  =  0.2  T 

- So  =  0.4  T 

- So  =  0.6  T 

- So  =  0.8  T 


Figure  F.15:  Monodimensional  LPM  -  Current  Drive 

Relative  Error  -  P  for  a  =  10^  S/m  -  rout  =  50mm 


- fl„  =  0.2  T 

- Bq  =  0.4  T 

- Bq  =  0.6  T 

- Bq  =  0.8  T 


Figure  F.16:  Monodimensional  LPM  -  Current  Drive 

Relative  Error  -  F  for  a  =  10^  S/m  -  rout  =  50mm 
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App.  F  Results:  Comparison  between  FDHM  and  LPM:  a  =  10^  S/m 


F.1.2  Voltage  Drive 


- _B„  =  0.2  T 

- B„  =  0.4  T 

- Bo  =  0.6  T 

- Bo  =  0.8  T 

o  FDHM 
<  LPM -ID 


Figure  F.17:  Monodimensional  LPM  -  Voltage  Drive 
r  for  a  —  10^  Sjm  -  rout  =  50mm 


- Bo  =  0.2  T 

- Bo  =  0.4  T 

- Bo  =  0.6r 

- Bo  =  0.8  T 

o  FDHM 
<  LPM -ID 


Figure  F.18:  Monodimensional  LPM  -  Voltage  Drive 
Tc  for  a  =  10^  Sjm  -  rout  =  50mm 


- Bo  =  0.2  T 

- Bo  =  0.4  T 

—Bo  =  0.6  T 

- Bo  =  0.8  T 

o  FDHM 
<  LPM -ID 


Figure  F.19:  Monodimensional  LPM  -  Voltage  Drive 
P  for  a  =  10^  Sjm  -  rout  =  50mm 


- Bo  =  0.2  T 

- Bo  =  0.4  T 

- Bo  =  0.6  T 

- Bo  =  0.8  T 

o  FDHM 
<  LPM -ID 


Figure  F.20:  Monodimensional  LPM  -  Voltage  Drive 
r  for  a  =  10^  S/m  -  Tout  =  50mm 
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App.  F  Results:  Comparison  between  FDHM  and  LPM:  a  =  10^  S/m 


h  =  10mm 


h  =  20mm 


h  =  30mm 


- Bq  =  0.2  T 

- So  =  0.4  T 

- So  =  0.6  T 

- So  =  0.8  T 


Figure  F.21:  Monodimensional  LPM  -  Voltage  Drive 

Relative  Error  -  F  for  a  =  10^  S/m  -  Tout  =  50mm 


- fl„  =  0.2  T 

- Bq  =  0.4  T 

- flo  =  0.6T 

- So  =  0.8  T 


Figure  F.22:  Monodimensional  LPM  -  Voltage  Drive 

Relative  Error  -  for  a  =  10^  S/m  -  Tout  =  50mm 


h  =  10mm 


L  —  [mm] 


X 10"^  ^  “  20mm 


L  —  [mm] 


X  10’^  ft.  =  30mm 


- So  =  0.2  T 

- So  =  0.4  T 

- So  =  0.6  T 

- So  =  0.8  T 


Figure  F.23:  Monodimensional  LPM  -  Voltage  Drive 

Relative  Error  -  P  for  a  =  10^  S/m  -  rout  =  50mm 


- fl„  =  0.2  T 

- Bq  =  0.4  T 

- Bq  =  0.6  T 

- Bq  =  0.8  T 


Figure  F.24:  Monodimensional  LPM  -  Voltage  Drive 

Relative  Error  -  F  for  a  =  10^  S/m  -  rout  =  50mm 
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—  [dimensionless] 


App.  F  Results:  Comparison  between  FDHM  and  LPM:  a  =  10^  S/m 


F.2  Bidimensional  LPM  -  a  =  10‘^S/m 


- Bo  =  0.2T 

- Bo  =  0.4  T 

—Bo  =  0.6  T 

- Bo  =  0.8  T 
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<  LPM -2D 


Figure  F.25:  Bidimensional  LPM 
Kvis  for  a  =  10^  S/m  -  Tout  =  50mm 

X  10'"^  ^  “  lOmm  h  =  2Umm  ^  jq-4  h  =  30mm 


- Bo  =  0.2  T 

- Bo  =  0.4  T 

- Bo  =  0.6  T 

- Bo  =  0.8  T 
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Figure  F.26:  Bidimensional  LPM 
Kj  for  a  =  10^  Sjm  -  Tout  =  50mm 


- Bo  =  0.2T 

- Bo  =  0.4  T 

- Bo  =  0.6  T 

- Bo  =  0.8  T 

o  FDHM 
4  LPM  -  2D 
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Figure  F.27:  Bidimensional  LPM 
Cui  for  <j  =  10^  S/m  -  rout  =  50mm 

^  /i  =  2Umm 


^  4  *■ 
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2  3  4 

L  —  \mm\ 


2  3  4 

L  —  \mm] 


- Bo  =  0.2  T 

- Bo  =  0.4  T 

- — -Bo  =  0.6r 

- Bo  =  0.8  T 

o  FDHM 
4  LPM -2D 


Figure  F.28:  Bidimensional  LPM 
Kv  for  a  =  10^  S/m  -  Tout  =  50mm 
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App.  F  Results:  Comparison  between  FDHM  and  LPM:  a  =  10^  S/m 


h  =  IQmm 


h  =  20mm 


h  =  30mm 


- Bq  =  0.2  T 

- So  =  0.4  T 

- So  =  0.6  T 

- So  =  0.8  T 


Figure  F.29:  Bidimensional  LPM 

Relative  Error  -  Kyis  for  a  =  10^  S/m  -  rout  =  SOmm 


Figure  F.30:  Bidimensional  LPM 

Relative  Error  -  Kj  for  a  =  10^  S/m  -  rout  =  SOmm 


- Bo  =  0.2  T 

- Bo  =  0.4  T 

- Bo  =  0.6  T 

- Bo  =  0.8  T 


Figure  F.31:  Bidimensional  LPM 

Relative  Error  -  c^y  for  a  =  10^  S/m  -  rout  =  50mm 


X 10'^  ^  “  107rim 
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X  lO  "^  =  30mm 


- So  =  0.2  T 

- So  =  0.4  T 

- So  =  0.6  T 

- So  =  0.8  T 


Figure  F.32:  Bidimensional  LPM 

Relative  Error  -  Ky  for  a  =  10^  S/m  -  rout  =  50mm 
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App.  F  Results:  Comparison  between  FDHM  and  LPM:  a  =  10^  S/m 


F.2.1  Current  Drive 


Figure  F.33:  Bidimensional  LPM  -  Current  Drive 
r  for  a  —  10^  Sjm  -  Vout  =  50mm 


- So  =  0.2  T 

- So  =  0.4  T 

— So  =  o.6r 

- So  =  0.8  T 
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<  LPM  -  2D 


Figure  F.34:  Bidimensional  LPM  -  Current  Drive 
Tj,  for  a  =  10^  S/m  -  Vout  =  50mm 
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L  —  \mm] 


h  =  2U7nm 


L  —  Imm] 


h  =  30mm 


L  —  \mm] 
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Figure  F.35:  Bidimensional  LPM  -  Current  Drive 
P  for  a  =  10^  S/m  -  rout  =  50mm 


X 10'"^  ^  “  10mm 


L  —  \mm\ 
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L  —  Imm] 
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<  LPM  -  2D 


Figure  F.36:  Bidimensional  LPM  -  Current  Drive 
r  for  a  =  10^  S/m  -  rout  =  50mm 
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App.  F  Results:  Comparison  between  FDHM  and  LPM:  a  =  10^  S/m 


h  =  IQmm 


h  =  20mm 


h  =  30mm 


- Bq  =  0.2  T 

- So  =  0.4  T 

- So  =  0.6  T 

- So  =  0.8  T 


Figure  F.37:  Bidimensional  LPM  -  Current  Drive 

Relative  Error  -  F  for  a  =  10^  S/m  -  rout  =  50mm 


h  =  10mm 


L  —  [mm] 


- So  =  0.2  T 

- So  =  0.4  T 

- — -So  =  o.6r 
- So  =  0.8  T 


Figure  F.38:  Bidimensional  LPM  -  Current  Drive 

Relative  Error  -  To  for  a  =  10^  S/m  -  Tout  =  50mm 


- fl„  =  0.2  T 

- Bq  =  0.4  T 

- Bq  =  0.6  T 

- So  =  0.8  T 


Figure  F.39:  Bidimensional  LPM  -  Current  Drive 

Relative  Error  -  P  for  a  =  10^  S/m  -  rout  =  50mm 


- Bo  =  0.2  T 

- Bo  =  0.4  T 

- Bo  =  0.6  T 

- Bo  =  0.8  T 


Figure  F.40:  Bidimensional  LPM  -  Current  Drive 

Relative  Error  -  F  for  a  =  10^  S/m  -  rout  =  50mm 
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App.  F  Results:  Comparison  between  FDHM  and  LPM:  a  =  10^  S/m 


F.2.2  Voltage  Drive 


Figure  F.41:  Bidimensional  LPM  -  Voltage  Drive 
r  for  a  —  10^  S/m  -  Tout  =  50mm 
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Figure  F.42:  Bidimensional  LPM  -  Voltage  Drive 
Tj,  for  a  =  10^  S/m  -  rout  =  50mm 
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Figure  F.43:  Bidimensional  LPM  -  Voltage  Drive 
P  for  a  =  10^  Sjm  -  rout  =  50mm 
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- So  =  0.2  T 

- So  =  0.4  T 

- So  =  0.6r 
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o  FDHM 
<  LPM  -  2D 


- So  =  0.2  T 

- So  =  0.4  T 

—  So  =  0.6T 

- So  =  0.8  T 

o  FDHM 
<  LPM  -  2D 


- So  =  0.2  T 

- So  =  0.4  T 

- So  =  0.6  T 
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Figure  F.44:  Bidimensional  LPM  -  Voltage  Drive 
r  for  a  =  10^  Sjm  -  rout  =  50mm 
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App.  F  Results:  Comparison  between  FDHM  and  LPM:  a  =  10^  S/m 


h  =  10mm 


h  =  20mm 


h  =  30mm 


- Bq  =  0.2  T 

- So  =  0.4  T 

- So  =  0.6r 

- So  =  0.8  T 


Figure  F.45:  Bidimensional  LPM  -  Voltage  Drive 

Relative  Error  -  F  for  a  =  10^  S/m  -  rout  =  50mm 


- Bo  =  0.2  T 

- Bo  =  0.4  T 

- So  =  0.6  1’ 

- Bo  =  0.8  T 


Figure  F.46:  Bidimensional  LPM  -  Voltage  Drive 

Relative  Error  -  Tc  for  a  =  10^  S/m  -  rout  =  50mm 


- Bo  =  0.2  T 

- Bo  =  0.4  T 

- Bo  =  O.ST 

- Bo  =  0.8  T 


Figure  F.47:  Bidimensional  LPM  -  Voltage  Drive 

Relative  Error  -  P  for  a  =  10^  S/m  -  rout  =  50mm 
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Figure  F.48:  Bidimensional  LPM  -  Voltage  Drive 

Relative  Error  -  F  for  a  —  10^  S/m  -  rout  =  50mm 


168 


DISTRIBUTION  A.  Approved  for  public  release:  distribution  unlimited. 


Appendix 


G 


Comparison  between  FDHM  and  LPM: 

a  =  10^  S/m 


G.l  Monodimensional  LPM  -  a  =  10^ S/m  .  170 

G.1.1  Gurrent  Drive .  172 

G.l. 2  Voltage  Drive  .  174 

G.2  Bidimensional  LPM  -  a  =  lO^S'/m .  176 
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App.  G  Results:  Comparison  between  FDHM  and  LPM:  a  =  10^  S/m 


G.l  Mono  dimensional  LPM  -  a  =  10^  S/m 
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Figure  G.l:  Kyis  for  a  =  10^  S/m  -  Tout  =  50mm  -  Monodimensional  LPM 
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Figure  G.2:  Monodimensional  LPM 
Kj  for  a  =  10^  S/m  -  Tout  =  50mm 


- Bo  =  0.2  T 

- Bo  =  0.4  T 

- Bo  =  0.6  T 

- Bo  =  0.8  T 
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Figure  G.3:  Monodimensional  LPM 
Cuj  for  a  =  10^  S/m  -  rout  =  50mm 


- So  =  0.2T 

- Bo  =  0.4  T 
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o  FDHM 
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Figure  G.4:  Monodimensional  LPM 
Ky  for  a  =  10^  S/m  -  rout  =  50mm 
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Err^el  -  c„  EWrel  -  Kj 


App.  G  Results:  Comparison  between  FDHM  and  LPM:  a  =  10^  S/m 
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Figure  G.5:  Monodimensional  LPM 

Relative  Error  -  Kyis  for  a  =  10^  S/m  -  Vout  =  50mm 
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Figure  G.6:  Monodimensional  LPM 

Relative  Error  -  Kj  for  a  =  10^  S/m  -  rout  =  50mm 


Figure  G.7:  Monodimensional  LPM 

Relative  Error  -  c^j  for  a  =  10^  S/m  -  rout  =  50mm 
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Figure  G.8:  Monodimensional  LPM 

Relative  Error  -  Ky  for  a  =  10^  S/m  -  rout  =  50mm 
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App.  G  Results:  Comparison  between  FDHM  and  LPM:  a  =  10^  S/m 


G.1.1  Current  Drive 
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G.9:  Monodimensional  LPM  -  Current  Drive 
r  for  a  =  10^  S/m  -  rout  =  50mm 
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G.ll:  Monodimensional  LPM  -  Current  Drive 
P  for  a  =  10^  S/m  -  rout  =  50mm 
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Figure  G.12:  Monodimensional  LPM  -  Current  Drive 
r  for  a  =  10^  S/m  -  rout  =  50mm 
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Errrei  -  f  Errrei  p  ^  Errrei  Tc  Errrei  -  Torque 


App.  G  Results:  Comparison  between  FDHM  and  LPM:  a  =  10^  S/m 
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Figure  G.13:  Monodimensional  LPM  -  Current  Drive 

Relative  Error  -  F  for  a  =  10^  S/m  -  rout  =  50mm 
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Figure  G.14:  Monodimensional  LPM  -  Current  Drive 

Relative  Error  -  Tc  for  a  =  10^  S/m  -  rout  =  50mm 


- fl„  =  0.2  T 

- Bq  =  0.4  T 

- Bq  =  0.6  T 

- So  =  0.8  T 


Figure  G.15:  Monodimensional  LPM  -  Current  Drive 

Relative  Error  -  P  for  a  =  10^  S/m  -  rout  =  50mm 


- Bo  =  0.2  T 

- Bo  =  0.4  T 

- Bo  =  0.6  T 

- Bo  =  0.8  T 


Figure  G.16:  Monodimensional  LPM  -  Current  Drive 

Relative  Error  -  F  for  a  =  10^  S/m  -  rout  =  50mm 
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App.  G  Results:  Comparison  between  FDHM  and  LPM:  a  =  10^  S/m 


G.1.2  Voltage  Drive 


X  10'^  ^  “  10777,771 


L  —  fmml 


X  10"^  ^  “  20777777 


L  —  \mm\ 


XlO'^  /l  =  30777777 


L  —  \mm\ 


- So  =  0.2  T 

- So  =  0.4  T 

- So  =  0.6r 

- So  =  0.8  T 

o  FDHM 
<  LPM -ID 


Figure  G.17:  Monodimensional  LPM  -  Voltage  Drive 
r  for  a  —  10^  Sjm  -  rout  =  50mm 


- So  =  0.2  T 

- So  =  0.4  T 

- So  =  0.6  T 

- So  =  0.8  T 

o  FDHM 
<  LPM -ID 


Figure  G.18:  Monodimensional  LPM  -  Voltage  Drive 
Tc  for  a  =  10^  Sjm  -  rout  =  50mm 


- So  =  0.2  T 

- So  =  0.4  T 

—  So  =  0.6T 

- So  =  0.8  T 

o  FDHM 
<  LPM -ID 


Figure  G.19:  Monodimensional  LPM  -  Voltage  Drive 
P  for  a  =  10^  S/m  -  rout  =  50mm 


X  10'"^  ^  “  IO777777 


L  —  fmml 


X 10"^  ^  ~  2O777777 


L  —  fmml 


X  h  =  30777 777 


L  —  fmml 


- So  =  0.2  T 

- So  =  0.4  T 

- So  =  0.6  T 

- So  =  0.8  T 

o  FDHM 
<  LPM -ID 


Figure  G.20:  Monodimensional  LPM  -  Voltage  Drive 
r  for  a  =  10^  S/m  -  rout  =  50mm 
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App.  G  Results:  Comparison  between  FDHM  and  LPM:  a  =  10^  S/m 


- Bq  =  0.2  T 

- So  =  0.4  T 

-  So  =  0.6  T 

- So  =  0.8  T 


Figure  G.21:  Monodimensional  LPM  -  Voltage  Drive 

Relative  Error  -  F  for  a  =  10^  S/m  -  rout  =  50mm 


h  =  10mm 


L  —  [mm] 


h  =  20mm 


L  —  [mm] 


h  =  30mm 


- So  =  0.2  T 

- So  =  0.4  T 

- So  =  0.6r 

- So  =  0.8  T 


Figure  G.22:  Monodimensional  LPM  -  Voltage  Drive 

Relative  Error  -  Tc  for  a  =  10^  S/m  -  rout  =  50mm 


- fl„  =  0.2  T 

- Bq  =  0.4  T 

- Bq  =  0.6  T 

- So  =  0.8  T 


Figure  G.23:  Monodimensional  LPM  -  Voltage  Drive 

Relative  Error  -  P  for  a  =  10^  S/m  -  rout  =  50mm 


h  =  30mm 


- So  =  0.2  T 

- So  =  0.4  T 

- So  =  0.6  T 

- So  =  0.8  T 


Figure  G.24:  Monodimensional  LPM  -  Voltage  Drive 

Relative  Error  -  F  for  a  =  10^  S/m  -  rout  =  50mm 
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App.  G  Results:  Comparison  between  FDHM  and  LPM:  a  =  10^  S/m 


G.2  Bidimensional  LPM  -  a  =  lO^S/m 


- Bo  =  0.2T 

- Bo  =  0.4  T 

—Bo  =  0.6  T 

- Bo  =  0.8  T 

o  FDHM 
<  LPM -2D 


X  10'"^  ^  “  10mm 


1  2  3  4  5 

L  —  \mm\ 


Figure  G.25:  Bidimensional  LPM 
Kvis  for  a  =  10^  S/m  -  Tout  =  50mm 


X 10“^  ^  ~  20mm  ^  jq-4  h  =  30mm 


- Bo  =  0.2  T 

- Bo  =  0.4  T 

- Bo  =  0.6  T 

- Bo  =  0.8  T 

o  FDHM 
4  LPM  -  2D 


Figure  G.26:  Bidimensional  LPM 
Kj  for  a  =  10^  S/m  -  Tout  =  50mm 

h  =  107nm  h  =  2Umm  h  =  30mm 


L  —  \mm]  L  —  \mm]  L  —  \mm] 


- Bo  =  0.2T 

- Bo  =  0.4  T 

- Bo  =  0.6  T 

- Bo  =  0.8  T 

o  FDHM 
4  LPM  -  2D 


1  2  3  4  5 

L  —  \mm] 


Figure  G.27:  Bidimensional  LPM 
Cui  for  <j  =  10^  S/m  -  rout  =  50mm 


X  10~^  ^  “  20mm 


1  2  3  4  5 

L  —  \mm] 


- Bo  =  0.2  T 

- Bo  =  0.4  T 

- — -Bo  =  0.6r 

- Bo  =  0.8  T 

o  FDHM 
4  LPM -2D 


Figure  G.28:  Bidimensional  LPM 
Kv  for  a  —  10^  S/m  -  Tout  =  50mm 
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App.  G  Results:  Comparison  between  FDHM  and  LPM:  a  =  10^  S/m 


- Bq  =  0.2  T 

- So  =  0.4  T 

— So  =  o.6r 
- So  =  0.8  T 


Figure  G.29:  Bidimensional  LPM 

Relative  Error  -  Kyis  for  a  =  10^  S/m  -  Vout  =  50mm 


^  jq-13  h  =  10mm 


L  —  [mm] 


^jq-12  h  =  20mm 


L  —  [mm] 


1  2  3  4  5 

L  —  [mm] 


Figure  G.30:  Bidimensional  LPM 

Relative  Error  -  Kj  for  a  =  10^  S/m  -  Vout  =  50mm 


X  10'^  ^  “  10777,772 


L  —  [mm] 


X  10’^  ^  “  20772771  jq-3  h  =  30777772 


- So  =  0.2  T 

- So  =  0.4  T 

- So  =  0.6  T 

- So  =  0.8  T 


Figure  G.31:  Bidimensional  LPM 

Relative  Error  -  for  a  =  10^  S/m  -  rout  =  50mm 


X  10'^  ^  “  10772772 


L  —  [mm] 


X  10"^  ^  ~  20772772 


L  —  [mm] 


X  lO  "^  /l  =  30771772 


- So  =  0.2  T 

- So  =  0.4  T 

- So  =  0.6  T 

- So  =  0.8  T 


Figure  G.32:  Bidimensional  LPM 

Relative  Error  -  Ky  for  a  =  10^  S/m  -  rout  =  50mm 
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App.  G  Results:  Comparison  between  FDHM  and  LPM:  a  =  10^  S/m 


G.2.1  Current  Drive 


- _B„  =  0.2  T 

- B„  =  0.4  T 

- - =  0.6  T 

- Bq  =  0.8  T 

o  FDHM 
<  LPM -ID 


Figure  G.33:  Bidimensional  LPM  -  Current  Drive 
r  for  a  =  10^  Sjm  -  rout  =  50mm 


- So  =  0.2  T 

- So  =  0.4  T 

— So  =  o.6r 

- So  =  0.8  T 

o  FDHM 
<  LPM  -  2D 


Figure  G.34:  Bidimensional  LPM  -  Current  Drive 
Tc  for  a  =  10^  S/m  -  r^ut  =  50mm 


h  =  10mm  fi  =  ‘20mm 


L  —  \mm]  L  —  Imm] 


h  =  30mm 


L  —  \mm] 


- So  =  0.2  T 

- So  =  0.4  T 

—  So  =  0.6T 

- So  =  0.8  T 

o  FDHM 
<  LPM -ID 


Figure  G.35:  Bidimensional  LPM  -  Current  Drive 
P  for  a  =  10^  S/m  -  rout  =  50mm 

X  10'"^  ^  “  10mm  ^  jq-4  fi  =  2Umm  ^  jq-4  h  =  30mm 


L  —  \mm\ 


- So  =  0.2  T 

- So  =  0.4  T 

- So  =  0.6  T 

- So  =  0.8  T 

o  FDHM 
<  LPM  -  2D 


Figure  G.36:  Bidimensional  LPM  -  Current  Drive 
r  for  a  =  10^  S/m  -  rout  =  50mm 
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App.  G  Results:  Comparison  between  FDHM  and  LPM:  a  =  10^  S/m 


h  =  10mm 


h  =  20mm 


h  =  30mm 


L  —  [mm]  L  —  [mm] 


- Bq  =  0.2  T 

- So  =  0.4  T 

— So  =  o.6r 
- So  =  0.8  T 


Figure  G.37:  Bidimensional  LPM  -  Current  Drive 

Relative  Error  -  F  for  a  =  10^  S/m  -  rout  =  50mm 


h  =  10mm 


L  —  [mm] 


- So  =  0.2  T 

- So  =  0.4  T 

- — -So  =  o.6r 
- So  =  0.8  T 


Figure  G.38:  Bidimensional  LPM  -  Current  Drive 

Relative  Error  -  Tc  for  a  =  10^  S/m  -  rout  =  50mm 


- fl„  =  0.2  T 

- Bq  =  0.4  T 

- Bq  =  0.6  T 

- So  =  0.8  T 


Figure  G.39:  Bidimensional  LPM  -  Current  Drive 

Relative  Error  -  P  for  a  =  10^  S/m  -  rout  =  50mm 


- Bo  =  0.2  T 

- Bo  =  0.4  T 

- Bo  =  0.6  T 

- Bo  =  0.8  T 


Figure  G.40:  Bidimensional  LPM  -  Current  Drive 

Relative  Error  -  F  for  a  =  10^  S/m  -  rout  =  50mm 
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App.  G  Results:  Comparison  between  FDHM  and  LPM:  a  =  10^  S/m 


G.2.2  Voltage  Drive 


h  =  10777,771 


n-5  h  =  20777777 


XlO'^  /l  =  30777777 


L  —  \mm] 

Figure  G.41:  Bidimensional  LPM  -  Voltage  Drive 
r  for  a  =  10^  S/m  -  rout  =  50mm 


h  =  IO777777 


h  =  2O777777 


L  —  \mm] 


h  =  3O777777 


2 

_  1.5 
CfL 

I  1 

G 

0.5 

0 


L  —  \mm] 


L  —  Imm] 


Figure  G.42:  Bidimensional  LPM  -  Voltage  Drive 
To  for  a  =  10^  S/m  -  rout  =  50mm 


h  =  IO777777 


h  =  2O777777 


L  —  \mm] 


h  =  3O777777 


L  —  Imm] 

Figure  G.43:  Bidimensional  LPM  -  Voltage  Drive 
P  for  a  =  10^  S/m  -  Tout  =  50mm 


L  —  fmml 


y4  h  =  IO777777 


Y-i  /7  =  2O777777 


^10-4  h  =  30777 777 


CO 

I 


L  —  fmml 


L  —  fmml 


L  —  fmml 


Figure  G.44:  Bidimensional  LPM  -  Voltage  Drive 
r  for  a  =  10^  Sjm  -  rout  =  50mm 


- Bo  =  0.2  T 

- Bo  =  0.4  T 

- Bo  =  0.6  T 

- Bo  =  0.8  T 

o  FDHM 
<  LPM  -  2D 


- Bo  =  0.2  T 

- Bo  =  0.4  T 

- Bo  =  0.6r 

- Bo  =  0.8  T 

o  FDHM 
<  LPM  -  2D 


- Bo  =  0.2  T 

- Bo  =  0.4  T 

—  Bo  =  0.6T 

- Bo  =  0.8  T 

o  FDHM 
<  LPM  -  2D 


- Bo  =  0.2  T 

- Bo  =  0.4  T 

- Bo  =  0.6  T 

- Bo  =  0.8  T 

o  FDHM 
<  LPM  -  2D 
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App.  G  Results:  Comparison  between  FDHM  and  LPM:  a  =  10^  S/m 


- Bq  =  0.2  T 

- So  =  0.4  T 

— So  =  o.6r 
- So  =  0.8  T 


Figure  G.45:  Bidimensional  LPM  -  Voltage  Drive 

Relative  Error  -  T  for  a  =  10^  S/m  -  rout  =  50mm 


- So  =  0.2  T 

- So  =  0.4  T 

— So  =  o.6r 
- So  =  0.8  T 


Figure  G.46:  Bidimensional  LPM  -  Voltage  Drive 

Relative  Error  -  Tc  for  a  =  10^  S/m  -  rout  =  SOmm 


- fl„  =  0.2  T 

- Bq  =  0.4  T 

- Bq  =  O.ST 

- Bq  =  0.8  T 


Figure  G.47:  Bidimensional  LPM  -  Voltage  Drive 

Relative  Error  -  P  for  a  =  10^  S/m  -  rout  =  50mm 


- Bo  =  0.2  T 

- Bo  =  0.4  T 

- So  =  0.6  1’ 

- Bo  =  0.8  T 


Figure  G.48:  Bidimensional  LPM  -  Voltage  Drive 

Relative  Error  -  F  for  a  =  10^  S/m  -  rout  =  50mm 
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Appendix 


H 


Comparison  between  FDHM  and  LPM: 

a  =  10^  S/m 


H.l  Monodimensional  LPM  -  a  =  10"^ S'/m  .  183 

H.1.1  Current  Drive .  185 

H.1.2  Voltage  Drive  .  187 

H.2  Bidimensional  LPM  -  a  =  lO^S'/m .  189 

H.2.1  Current  Drive .  191 

H.2. 2  Voltage  Drive  .  193 
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App.  H  Results:  Comparison  between  FDHM  and  LPM:  a  =  10'^  S/m 


H.l  Mono  dimensional  LPM  -  a  =  lO^S/m 


h  =  lOmm 


h  =  ‘2{)mm 


h  =  30mm 


I 


- Bq  =  0.2  T 

- Bo  =  0.4  T 

- Bo  =  0.6r 

- Bo  =  0.8  T 

o  FDHM 
<  LPM -ID 


L  —  fmm] 


-  [mml 


L  —  fmml 


Figure  H.l:  Kyis  for  a  =  10^  S/m  -  Tout  =  50mm  -  Monodimensional  LPM 


-  lOmm 


h  =  2Umm 


3  4 

•  fmm] 


- So  =  0.2T 

- Bo  =  0.4  T 

- Bo  =  0.6r 

- Bo  =  0.8  T 

o  FDHM 
<  LPM -ID 


Figure  H.2:  Monodimensional  LPM 
Kj  for  a  =  10^  S/m  -  Tout  =  50mm 


Figure  H.3:  Monodimensional  LPM 
Cuj  for  a  =  10^  S/m  -  rout  =  50mm 


- Bo  =  0.2  T 

- Bo  =  0.4  T 

- Bo  =  0.6  T 

- Bo  =  0.8  T 

o  FDHM 
<  LPM  -  ID 


Figure  H.4:  Monodimensional  LPM 
Kv  for  a  =  10^  S/m  -  rout  =  50mm 
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App.  H  Results:  Comparison  between  FDHM  and  LPM:  a  =  10'^  S/m 


h  =  IQmm 


h  =  20mm 


h  =  30mm 


- Bq  =  0.2  T 

- So  =  0.4  T 

- So  =  0.6  T 

- So  =  0.8  T 


Figure  H.5:  Monodimensional  LPM 

Relative  Error  -  Kyis  for  a  =  10^  S/m  -  Tout  =  50mm 


- fl„  =  0.2  T 

- Bq  =  0.4  T 

- flo  =  0.6T 

- So  =  0.8  T 


Figure  H.6:  Monodimensional  LPM 

Relative  Error  -  Kj  for  a  =  10^  S/m  -  Tout  =  50mm 


X 10“^  ^  “  107rim 


X 10'"^  ^  “  20mm 


L  —  [mm] 


jq-4  h  =  30mm 


- So  =  0.2  T 

- So  =  0.4  T 

- So  =  0.6  T 

- So  =  0.8  T 


Figure  H.7:  Monodimensional  LPM 

Relative  Error  -  c^j  for  a  =  10^  S/m  -  rout  =  50mm 


X 10'^  ^  “  107rim 


L  —  [mm] 


X  ^  “  20mm  jq-4  h  =  30mm 


- So  =  0.2  T 

- So  =  0.4  T 

- So  =  0.6  T 

- So  =  0.8  T 


Figure  H.8:  Monodimensional  LPM 

Relative  Error  -  Ky  for  a  =  10"^  S/m  -  rout  =  50mm 


184 


DISTRIBUTION  A.  Approved  for  public  release:  distribution  unlimited. 


App.  H  Results:  Comparison  between  FDHM  and  LPM:  a  =  10'^  S/m 


H.1.1  Current  Drive 


- _B„  =  0.2  T 

- B„  =  0.4  T 

- Bo  =  0.6r 

- Bq  =  0.8  T 

o  FDHM 
<  LPM -ID 


Figure  H.9:  Monodimensional  LPM  -  Current  Drive 
r  for  a  =  10^  Sfm  -  Vout  =  50mm 


- So  =  0.2  T 

- So  =  0.4  T 

— So  =  o.6r 

- So  =  0.8  T 

o  FDHM 
<  LPM -ID 


Figure  H.IO:  Monodimensional  LPM  -  Current  Drive 
Tc  for  a  =  10^  Sfm  -  rout  =  50mm 


- So  =  0.2  T 

- So  =  0.4  T 

—  So  =  0.6T 

- So  =  0.8  T 

o  FDHM 
<  LPM -ID 


Figure  H.ll:  Monodimensional  LPM  -  Current  Drive 
P  for  a  =  10^  S/m  -  rout  =  50mm 


X 10'"^  ^  “  10mm 


L  —  \mm\ 


X 10“^  ^  “  ‘20mm 


L  —  Imm] 


X  h  =  30mm 


L  —  \mm] 


- So  =  0.2  T 

- So  =  0.4  T 

- So  =  0.6  T 

- So  =  0.8  T 

o  FDHM 
<  LPM -ID 


Figure  H.12:  Monodimensional  LPM  -  Current  Drive 
f  for  a  =  10"^  S/m  -  rout  =  SOmm 
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App.  H  Results:  Comparison  between  FDHM  and  LPM:  a  =  10'^  S/m 


h  =  IQmm 


h  =  20mm 


h  =  30mm 


- Bq  =  0.2  T 

- So  =  0.4  T 

- So  =  0.6  T 

- So  =  0.8  T 


Figure  H.13:  Monodimensional  LPM  -  Current  Drive 

Relative  Error  -  T  for  a  =  10^  S/m  -  rout  =  50mm 


h  =  10mm 


L  —  [mm] 


h  =  20mm 


L  —  [mm] 


h  =  30mm 


- So  =  0.2  T 

- So  =  0.4  T 

- — -So  =  o.6r 
- So  =  0.8  T 


Figure  H.14:  Monodimensional  LPM  -  Current  Drive 

Relative  Error  -  Tc  for  a  =  10"^  S/m  -  rout  =  50mm 


Figure  H.15:  Monodimensional  LPM  -  Current  Drive 

Relative  Error  -  P  for  a  =  10^  S/m  -  rout  =  50mm 


- fl„  =  0.2  T 

- Bq  =  0.4  T 

- Bq  =  0.6  T 

- Bq  =  0.8  T 


Figure  H.16:  Monodimensional  LPM  -  Current  Drive 

Relative  Error  -  T  for  a  =  10^  S/m  -  rout  =  50mm 
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App.  H  Results:  Comparison  between  FDHM  and  LPM:  a  =  10'^  S/m 


H.1.2  Voltage  Drive 


- _B„  =  0.2  T 

- B„  =  0.4  T 

- Bo  =  0.6r 

- Bq  =  0.8  T 

o  FDHM 
<  LPM -ID 


Figure  H.17:  Monodimensional  LPM  -  Voltage  Drive 
f  for  a  —  10^  Sjm  -  rout  =  50mm 


- So  =  0.2  T 

- So  =  0.4  T 

— So  =  o.6r 

- So  =  0.8  T 

o  FDHM 
<  LPM -ID 


Figure  H.18:  Monodimensional  LPM  -  Voltage  Drive 
Tc  for  a  =  10^  Sjm  -  rout  =  50mm 


- So  =  0.2  T 

- So  =  0.4  T 

—So  =  0.6  T 

- So  =  0.8  T 

o  FDHM 
<  LPM -ID 


Figure  H.19:  Monodimensional  LPM  -  Voltage  Drive 
P  for  a  =  10^  S/m  -  rout  =  50mm 


- So  =  0.2  T 

- So  =  0.4  T 

- So  =  0.6  T 

- So  =  0.8  T 

o  FDHM 
<  LPM -ID 


Figure  H.20:  Monodimensional  LPM  -  Voltage  Drive 
f  for  a  =  10"^  S/m  -  rout  =  SOmm 
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App.  H  Results:  Comparison  between  FDHM  and  LPM:  a  =  10'^  S/m 


h  =  10mm 


h  =  20mm 


h  =  30mm 


- Bq  =  0.2  T 

- So  =  0.4  T 

- So  =  0.6r 

- So  =  0.8  T 


Figure  H.21:  Monodimensional  LPM  -  Voltage  Drive 

Relative  Error  -  T  for  a  =  10^  S/m  -  Tout  =  50mm 


h  =  10mm 


L  —  [mm] 


h  =  20mm 


L  —  [mm] 


h  =  30mm 


- So  =  0.2  T 

- So  =  0.4  T 

- — -So  =  o.6r 
- So  =  0.8  T 


Figure  H.22:  Monodimensional  LPM  -  Voltage  Drive 

Relative  Error  -  Tc  for  a  =  10"^  S/m  -  rout  =  50mm 


- fl„  =  0.2  T 

- Bq  =  0.4  T 

- Bq  =  0.6  T 

- So  =  0.8  T 


Figure  H.23:  Monodimensional  LPM  -  Voltage  Drive 

Relative  Error  -  P  for  a  =  10^  S/m  -  rout  =  50mm 


h  =  107rim 


L  —  [mm] 


h  =  20mm 


L  —  [mm] 


h  =  30mm 


- So  =  0.2  T 

- So  =  0.4  T 

- So  =  0.6  T 

- So  =  0.8  T 


Figure  H.24:  Monodimensional  LPM  -  Voltage  Drive 

Relative  Error  -  T  for  a  =  10^  S/m  -  rout  =  50mm 
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—  [dimensionless] 


App.  H  Results:  Comparison  between  FDHM  and  LPM:  a  =  10'^  S/m 


H.2  Bidimensional  LPM  -  a  =  lO^S/m 


X 10'^  ^  “  10mm 


X 10'*^  ^  “  2Umm 


X  10  h  =  30mm 
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<  LPM -2D 


Figure  H.25:  Bidimensional  LPM 
Kvis  for  a  =  10^  S/m  -  rout  =  50mm 

X  10'"^  ^  “  10mm  iq-4  h  =  2Umm  ^  iq-4  h  =  30mm 


- Bo  =  0.2  T 

- Bo  =  0.4  T 

- Bo  =  0.6  T 

- Bo  =  0.8  T 
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Figure  H.26:  Bidimensional  LPM 
Kj  for  a  =  10^  Sjm  -  Tout  =  50mm 


- Bo  =  0.2T 

- Bo  =  0.4  T 

- Bo  =  0.6  T 

- Bo  =  0.8  T 

o  FDHM 
4  LPM  -  2D 
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Figure  H.27:  Bidimensional  LPM 
Cuj  for  a  =  10"^  S/m  -  rout  =  50mm 


h  =  2Umm 


1  2  3  4  5 

L  —  \mm] 


- Bo  =  0.2  T 

- Bo  =  0.4  T 
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4  LPM -2D 


Figure  H.28:  Bidimensional  LPM 
Ky  for  a  —  10^  S/m  -  rout  =  50mm 
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App.  H  Results:  Comparison  between  FDHM  and  LPM:  a  =  10“^  S/m 


- Bq  =  0.2  T 

- So  =  0.4  T 

— So  =  o.6r 
- So  =  0.8  T 


Figure  H.29:  Bidimensional  LPM 

Relative  Error  -  Kyis  for  a  =  10^  S/m  -  rout  =  50mm 


- fl„  =  0.2  T 

- Bq  =  0.4  T 

- flo  =  0.6T 

- So  =  0.8  T 


Figure  H.30:  Bidimensional  LPM 

Relative  Error  -  Kj  for  a  =  10"^  S/m  -  rout  =  50mm 


Figure  H.31:  Bidimensional  LPM 

Relative  Error  -  c^j  for  a  =  10^  S/m  -  rout  =  50mm 


- Bo  =  0.2  T 

- Bo  =  0.4  T 

- Bo  =  0.6  T 

- Bo  =  0.8  T 


Figure  H.32:  Bidimensional  LPM 

Relative  Error  -  Ky  for  a  =  10^  S/m  -  rout  =  50mm 
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App.  H  Results:  Comparison  between  FDHM  and  LPM:  a  =  10'^  S/m 


H.2.1  Current  Drive 


- _B„  =  0.2  T 

- B„  =  0.4  T 

- Bo  =  0.6r 

- Bq  =  0.8  T 

o  FDHM 
<  LPM -ID 


Figure  H.33:  Bidimensional  LPM  -  Current  Drive 
r  for  a  =  S/m  -  Vout  =  50mm 
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Figure  H.34:  Bidimensional  LPM  -  Current  Drive 
Tc  for  a  =  10^  Sfm  -  rout  =  50mm 


- So  =  0.2  T 

- So  =  0.4  T 

—So  =  0.6  T 

- So  =  0.8  T 

o  FDHM 
<  LPM -ID 


Figure  H.35:  Bidimensional  LPM  -  Current  Drive 
P  for  (7  =  10^  5/m  -  rout  =  50mm 


- So  =  0.2  T 

- So  =  0.4  T 

- So  =  0.6  T 

- So  =  0.8  T 

o  FDHM 
<  LPM  -  2D 


Figure  H.36:  Bidimensional  LPM  -  Current  Drive 
f  for  a  =  10"^  S/m  -  rout  =  50mm 
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App.  H  Results:  Comparison  between  FDHM  and  LPM:  a  =  10“^  S/m 


- Bq  =  0.2  T 

- So  =  0.4  T 

- So  =  0.6  T 

- So  =  0.8  T 


Figure  H.37:  Bidimensional  LPM  -  Current  Drive 

Relative  Error  -  T  for  a  =  10^  S/m  -  rout  =  50mm 


h  =  lOmm 


L  —  [mm] 


h  =  20m77i 


L  —  [mm] 


h  =  30mm 


- So  =  0.2  T 

- So  =  0.4  T 

- — -So  =  o.6r 
- So  =  0.8  T 


Figure  H.38:  Bidimensional  LPM  -  Current  Drive 

Relative  Error  -  Tc  for  a  =  10"^  S/m  -  rout  =  50mm 


- Bo  =  0.2  T 

- Bo  =  0.4  T 

- Bo  =  0.6  T 

- Bo  =  0.8  T 


Figure  H.39:  Bidimensional  LPM  -  Current  Drive 

Relative  Error  -  P  for  a  =  10^  S/m  -  rout  =  50mm 


- Bo  =  0.2  T 

- Bo  =  0.4  T 

- Bo  =  0.6  T 

- Bo  =  0.8  T 


Figure  H.40:  Bidimensional  LPM  -  Current  Drive 

Relative  Error  -  T  for  a  =  10^  S/m  -  rout  =  50mm 
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App.  H  Results:  Comparison  between  FDHM  and  LPM:  a  =  10'^  S/m 


H.2.2  Voltage  Drive 


X  10'^  ^  “  10777,771 


L  —  \mm\ 


- So  =  0.2  T 

- So  =  0.4  T 

- So  =  0.6r 

- So  =  0.8  T 

o  FDHM 
<  LPM  -  2D 


Figure  H.41:  Bidimensional  LPM  -  Voltage  Drive 
f  for  a  =  10^  Sjm  -  rout  =  50mm 


- So  =  0.2  T 

- So  =  0.4  T 

— So  =  o.6r 

- So  =  0.8  T 

o  FDHM 
<  LPM  -  2D 


Figure  H.42:  Bidimensional  LPM  -  Voltage  Drive 
Tc  for  a  =  10^  Sjm  -  rout  =  50mm 


- So  =  0.2  T 

- So  =  0.4  T 

—So  =  0.6  T 

- So  =  0.8  T 

o  FDHM 
<  LPM  -  2D 


Figure  H.43:  Bidimensional  LPM  -  Voltage  Drive 
P  for  (7  =  10^  S/m  -  rout  =  50mm 


- So  =  0.2  T 

- So  =  0.4  T 

- So  =  0.6  T 

- So  =  0.8  T 

o  FDHM 
<  LPM  -  2D 


Figure  H.44:  Bidimensional  LPM  -  Voltage  Drive 
f  for  a  =  10'^  S/m  -  rout  =  50mm 
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Errrd  -  f  Errrei  -  P  Err^d  -  Err^d  -  Torque 


App.  H  Results:  Comparison  between  FDHM  and  LPM:  a  =  10'^  S/m 
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h  =  30mm 
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-  So  =  0.6  T 
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Figure  H.45:  Bidimensional  LPM  -  Voltage  Drive 
Relative  Error  -  T  for  a  =  10^  S/m  - 
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Figure  H.46:  Bidimensional  LPM  -  Voltage  Drive 

Relative  Error  -  Tc  for  a  =  S/m  -  rout  =  SOmm 


- Bo  =  0.2  T 

- Bo  =  0.4  T 

- Bo  =  O.ST 

- Bo  =  0.8  T 


Figure  H.47:  Bidimensional  LPM  -  Voltage  Drive 

Relative  Error  -  P  for  a  =  10^  S/m  -  rout  =  50mm 


h  =  10mm 


L  —  [mm] 


h  =  20mm 


L  —  [mm] 


h  =  30mm 
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- So  =  0.2  T 

- So  =  0.4  T 

— So  =  o.6r 
- So  =  0.8  T 


Figure  H.48:  Bidimensional  LPM  -  Voltage  Drive 

Relative  Error  -  T  for  a  —  10'^  S/m  -  rout  =  50mm 
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Appendix 

I  Comparison  between  FDHM  and  LPM: 

a  =  10^  S/m 


1.1  Monodimensional  LPM  -  a  =  10^ S/m  .  196 

1.1.1  Current  Drive .  198 

1.1.2  Voltage  Drive  . 200 

1.2  Bidimensional  LPM  -  a  =  lO^S/m . 202 

1.2.1  Current  Drive . 204 

1.2.2  Voltage  Drive  . 206 
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App.  I  Results:  Comparison  between  FDHM  and  LPM:  a  =  10®  S/m 


I.l  Mono  dimensional  LPM  -  a  =  lO^S/m 


v4  h  =  lOmm 
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y4  h  =  30mm 
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- Bo  =  0.4  T 

- Bo  =  0.6r 

- Bo  =  0.8  T 
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L  —  fmm] 

Figure  I.l:  Kyis  for  a  =  10^  S/m  -  rout  =  50mm  -  Monodimensional  LPM 
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h  =  2Umm 
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- Bo  =  0.4  T 
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Figure  1.2:  Monodimensional  LPM 
Kj  for  a  =  10®  S/m  -  rout  =  50mm 


h  =  2Umm 


h  =  30mm 


L  —  fmm] 
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- Bo  =  0.2  T 
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Figure  1.3:  Monodimensional  LPM 
Cui  for  a  =  10®  S/m  -  rout  =  50mm 


X 10"^  ^  “  10mm 


L  —  fmm] 


X  ^  “  2Umm  ^  jq-4  h  =  30mm 


- So  =  0.2T 

- Bo  =  0.4  T 
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Figure  1.4:  Monodimensional  LPM 
Kv  for  a  =  10®  S/m  -  rout  =  50mm 
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App.  I  Results:  Comparison  between  FDHM  and  LPM:  a  =  10®  S/m 


- fl„  =  0.2  T 

- So  =  0.4  r 

- So  =  0.6  T 

- So  =  0.8  T 


Figure  1.5:  Monodimensional  LPM 

Relative  Error  -  K^s  for  a  =  10®  S/m  -  Vout  =  50mm 


^  jq-12  h  =  10mm 
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Figure  1.6:  Monodimensional  LPM 

Relative  Error  -  Kj  for  a  =  10®  S/m  -  r^ut  =  50mm 


X  lO  "^  ^  “  10mm 


X 10^  ^  ~  20mm 


L  —  [mm] 


X  10^^  ft.  =  30mm 


Figure  1.7:  Monodimensional  LPM 

Relative  Error  -  for  a  =  10®  S/m  -  rout  =  50mm 


X 10'^  ^  “  107nm 


L  —  [mm] 
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L  —  [mm] 


^  jq-3  ft  =  30mm 


- So  =  0.2  T 

- So  =  0.4  T 

- So  =  0.6  T 

- So  =  0.8  T 


Figure  1.8:  Monodimensional  LPM 

Relative  Error  -  Ky  for  a  =  10®  S/m  -  Tout  =  50mm 
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App.  I  Results:  Comparison  between  FDHM  and  LPM:  a  =  10®  S/m 


1.1.1  Current  Drive 
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Figure  1.9:  Monodimensional  LPM  -  Current  Drive 
r  for  a  —  10^  Sjm  -  rout  =  50mm 


- So  =  0.2  T 
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— So  =  o.6r 
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Figure  1. 10:  Monodimensional  LPM  -  Current  Drive 
Tc  for  a  =  10^  Sjm  -  rout  =  50mm 
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Figure  I.ll:  Monodimensional  LPM  -  Current  Drive 
P  for  a  =  10^  Sjm  -  rout  =  50mm 


- So  =  0.2  T 

- So  =  0.4  T 

- So  =  0.6  T 

- So  =  0.8  T 

o  FDHM 
<  LPM -ID 


Figure  1.12:  Monodimensional  LPM  -  Current  Drive 
f  for  a  =  10^  Sjm  -  rout  =  50mm 
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Erryci  -  r  Ervya  -  P  Err^d  -  Err  yd  -  Torque 


App.  I  Results:  Comparison  between  FDHM  and  LPM:  a  =  10®  S/m 
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Figure  1.13:  Monodimensional  LPM  -  Current  Drive 

Relative  Error  -  F  for  a  =  10®  S/m  -  rout  =  50mm 


h  =  10mm 
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- — -So  =  o.6r 
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Figure  1.14:  Monodimensional  LPM  -  Current  Drive 

Relative  Error  -  for  a  =  10®  S/m  -  Vout  =  50mm 


- fl„  =  0.2  T 

- Bq  =  0.4  T 

- Bq  =  0.6  T 

- So  =  0.8  T 


Figure  1.15:  Monodimensional  LPM  -  Current  Drive 

Relative  Error  -  P  for  a  =  10®  S/m  -  rout  =  50mm 


h  =  10mm  h  =  20mm 


L  —  [mm]  L  —  [mm] 


h  =  30mm 


- So  =  0.2  T 

- So  =  0.4  T 

- So  =  0.6  T 

- So  =  0.8  T 


Figure  1.16:  Monodimensional  LPM  -  Current  Drive 

Relative  Error  -  F  for  a  =  10®  S/m  -  rout  =  50mm 


199 


DISTRIBUTION  A.  Approved  for  public  release:  distribution  unlimited. 


App.  I  Results:  Comparison  between  FDHM  and  LPM:  a  =  10®  S/m 


1.1.2  Voltage  Drive 
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Figure  1.17:  Monodimensional  LPM  -  Voltage  Drive 
r  for  a  =  10^  S/m  -  rout  =  50mm 


1  2  3  4  5 

L  —  fmml 


h  =  20777777  h  =  30777777 


12345  12345 

L  —  \mm\  L  —  \mm\ 


- So  =  0.2  T 

- So  =  0.4  T 

- So  =  0.6r 

- So  =  0.8  T 

o  FDHM 
<  LPM -ID 


Figure  1.18:  Monodimensional  LPM  -  Voltage  Drive 
Tc  for  a  =  10^  S/m  -  rout  =  50mm 
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- So  =  0.8  T 
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Figure  1.19:  Monodimensional  LPM  -  Voltage  Drive 
P  for  a  =  10^  S/m  -  rout  =  50mm 
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- So  =  0.8  T 
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<  LPM -ID 


Figure  1.20:  Monodimensional  LPM  -  Voltage  Drive 
f  for  a  =  10®  S/m  -  Tout  =  50mm 
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Ervyci  -  r  Ervya  -  P  Err^ei  -  ErVya  -  Torque 


App.  I  Results:  Comparison  between  FDHM  and  LPM:  a  =  10®  S/m 
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Figure  1.21:  Monodimensional  LPM  -  Voltage  Drive 

Relative  Error  -  F  for  a  =  10®  S/m  -  Tout  =  50mm 
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- — -So  =  o.6r 
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Figure  1.22:  Monodimensional  LPM  -  Voltage  Drive 

Relative  Error  -  for  a  =  10®  S/m  -  rout  =  50mm 


h  =  10mm 


h  =  20mm 


h  =  30mm 


L  —  [mm] 


- So  =  0.2  T 

- So  =  0.4  T 

- So  =  0.6  T 

- So  =  0.8  T 


Figure  1.23:  Monodimensional  LPM  -  Voltage  Drive 

Relative  Error  -  P  for  a  =  10®  S/m  -  Vgut  =  50mm 


- fl„  =  0.2  T 

- Bq  =  0.4  T 

- Bq  =  0.6  T 

- Bq  =  0.8  T 


Figure  1.24:  Monodimensional  LPM  -  Voltage  Drive 

Relative  Error  -  F  for  a  =  10®  S/m  -  rout  =  50mm 
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App.  I  Results:  Comparison  between  FDHM  and  LPM:  a  =  10®  S/m 


1.2  Bidimensional  LPM  -  a  =  lO^S/m 


Figure  1.25:  Bidimensional  LPM 
Kyis  for  a  =  10®  S/m  -  Tout  =  50mm 
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Figure  1.26:  Bidimensional  LPM 
Kj  for  a  =  10®  S/m  -  Tout  =  50mm 
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Figure  1.27:  Bidimensional  LPM 
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Figure  1.28:  Bidimensional  LPM 
Kv  for  a  —  10®  S/m  -  Tout  =  50mm 
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Figure  1.29:  Bidimensional  LPM 
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Figure  1.30:  Bidimensional  LPM 

Relative  Error  -  Kj  for  a  =  10®  S/m  -  rout  =  SOmm 
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Figure  1.31:  Bidimensional  LPM 

Relative  Error  -  c^j  for  a  =  10®  S/m  -  rout  =  50mm 
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Figure  1.32:  Bidimensional  LPM 

Relative  Error  -  Ky  for  a  =  10®  S/m  -  rout  =  50mm 
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1.2.1  Current  Drive 
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Figure  1.33:  Bidimensional  LPM  -  Current  Drive 
r  for  a  =  10^  Sjm  -  rout  =  50mm 
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Figure  1.34:  Bidimensional  LPM  -  Current  Drive 
Tc  for  a  =  10^  Sjm  -  rout  =  50mm 
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Figure  1.35:  Bidimensional  LPM  -  Current  Drive 
P  for  a  =  10^  Sjm  -  rout  =  50mm 
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Figure  1.36:  Bidimensional  LPM  -  Current  Drive 
r  for  a  =  10®  S/m  -  Tout  =  50mm 
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Figure  1.37:  Bidimensional  LPM  -  Current  Drive 

Relative  Error  -  F  for  a  =  10®  S/m  -  rout  =  50mm 
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Figure  1.38:  Bidimensional  LPM  -  Current  Drive 

Relative  Error  -  for  a  =  10®  S/m  -  rout  =  50mm 
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Figure  1.39:  Bidimensional  LPM  -  Current  Drive 

Relative  Error  -  P  for  a  =  10®  S/m  -  rout  =  50mm 
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Figure  1.40:  Bidimensional  LPM  -  Current  Drive 

Relative  Error  -  F  for  a  —  10®  S/m  -  rout  =  50mm 
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App.  I  Results:  Comparison  between  FDHM  and  LPM:  a  =  10®  S/m 
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Figure  1.41:  Bidimensional  LPM  -  Voltage  Drive 
t  for  a  =  10®  S/m  -  Tout  =  50mm 
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Figure  1.42:  Bidimensional  LPM  -  Voltage  Drive 
Tc  for  a  =  10^  Sjm  -  ro^t  =  50mm 
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Figure  1.43:  Bidimensional  LPM  -  Voltage  Drive 
P  for  a  =  10^  Sjm  -  rout  =  50mm 
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Figure  1.44:  Bidimensional  LPM  -  Voltage  Drive 
r  for  a  =  10^  Sjm  -  rout  =  50mm 
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App.  I  Results:  Comparison  between  FDHM  and  LPM:  a  =  10®  S/m 
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Figure  1.45:  Bidimensional  LPM  -  Voltage  Drive 

Relative  Error  -  F  for  a  =  10®  S/m  -  Tout  =  50mm 


h  =  10mm 


h  =  20mm 


L  —  [mm] 


h  =  30mm 


- So  =  0.2  T 

- So  =  0.4  T 

- So  =  0.6r 

- So  =  0.8  T 


Figure  1.46:  Bidimensional  LPM  -  Voltage  Drive 

Relative  Error  -  for  a  =  10®  S/m  -  rout  =  50mm 
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Figure  1.47 :  Bidimensional  LPM  -  Voltage  Drive 

Relative  Error  -  P  for  a  =  10®  S/m  -  rout  =  50mm 
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Figure  1.48:  Bidimensional  LPM  -  Voltage  Drive 

Relative  Error  -  F  for  a  —  10®  S/m  -  rout  =  50mm 
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